bzoj2326【HNOI2011】数学作业

来源：互联网发布：php 负载均衡解决方案编辑：程序博客网时间：2024/04/29 04:16

2326: [HNOI2011]数学作业

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 1648 Solved: 972
[Submit][Status][Discuss]

Description

（f[n]）（10^k 1 1）（f[n-1]）

（ n ）=（0 1 1）（ n-1 ）

（ 1 ）（0 0 1）（ 1 ）

然后分段矩阵乘法。这道题调了整整一晚上，忘记每一步运算都取模，这道题很容易超过long long。

#include<iostream>#include<cstdio>#include<cmath>#include<cstdlib>#include<cstring>#include<algorithm>#define F(i,j,n) for(int i=j;i<=n;i++)#define D(i,j,n) for(int i=j;i>=n;i--)#define ll long longusing namespace std;ll n,m;struct matrix{ll f[4][4];matrix(){memset(f,0,sizeof(f));}friend matrix operator *(const matrix &a,const matrix &b){matrix c;F(i,1,3) F(j,1,3) F(k,1,3) c.f[i][j]=(c.f[i][j]%m+a.f[i][k]%m*b.f[k][j]%m)%m;return c;}}ans,tmp;inline void calc(ll t,ll x){memset(tmp.f,0,sizeof(tmp.f));tmp.f[1][1]=t%m;tmp.f[1][2]=tmp.f[1][3]=tmp.f[2][2]=tmp.f[2][3]=tmp.f[3][3]=1;for(ll y=x-t/10+1;y;y>>=1,tmp=tmp*tmp) if (y&1) ans=tmp*ans;}int main(){scanf("%lld%lld",&n,&m);memset(ans.f,0,sizeof(ans.f));F(i,1,3) ans.f[i][i]=1;ll t=10;while (n>=t){calc(t,t-1);t*=10;}calc(t,n);printf("%lld\n",ans.f[1][3]);return 0;}

1 0