BZOJ 2326 数学作业(分段矩阵快速幂)
生活随笔
收集整理的這篇文章主要介紹了
BZOJ 2326 数学作业(分段矩阵快速幂)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
實際上,對于位數相同的連續段,可以用矩陣快速冪求出最后的ans,那么題目中一共只有18個連續段。
分段矩陣快速冪即可。
?
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<queue> #include<cmath> #define ll long long using namespace std; const int maxn=5; ll n,m; struct Mat { ll mat[maxn][maxn]; Mat() { memset(mat,0,sizeof(mat)); } Mat operator *(const Mat &b) { Mat c; for(int k=1;k<=3;k++) { for(int i=1;i<=3;i++) { for(int j=1;j<=3;j++) { c.mat[i][j]+=mat[i][k]*b.mat[k][j]; c.mat[i][j]%=m; } } } return c; } }; Mat ans; Mat cal(ll k,ll t) { Mat res; Mat b; ll y=t-k/10+1; b.mat[1][1]=k%m; b.mat[1][2]=b.mat[2][2]=b.mat[2][3]=b.mat[3][3]=b.mat[1][3]=1; for(int i=1;i<=3;i++)res.mat[i][i]=1; while(y) { if(y&1)res=res*b; y>>=1; b=b*b; } return res; } int main() { ans.mat[3][1]=1; scanf("%lld%lld",&n,&m); ll t=10; for(int i=1;i<=18;i++) { ans=cal(t,t-1)*ans; t*=10; if(t>n)break; } ans=cal(t,n)*ans; printf("%lld\n",ans.mat[1][1]); return 0; } View Code?
轉載于:https://www.cnblogs.com/lishiyao/p/6739149.html
總結
以上是生活随笔為你收集整理的BZOJ 2326 数学作业(分段矩阵快速幂)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: eclipse修改代码后都需要clean
- 下一篇: 字符串在内存中的存储——C语言进阶