//九度OJ 1081 清华09机试题之递推数列
//http://ac.jobdu.com/problem.php?pid=1081
//利用二分求幂法和蛮生疏的函数递归调用编写的矩阵的大数幂。
//相信N年后的我再回来看看这稚气未脱的刚入道时自己写的代码,会心一笑吧哈哈
#include <stdio.h>
#define N 10000
typedef struct M{
long long h[2][2];
}M;
M c,e;//e在main函数里面被设置成了单位矩阵。
long long a,b,s;
M copy(M a,M b)
{
a=b;
return a;
}
M mul(M a,M b)
{
M tempm;
for(int i=0;i<2;i++)
for(int j=0;j<2;j++)
{
tempm.h[i][j]=0;
for(int k=0;k<2;k++)tempm.h[i][j]+=(a.h[i][k])*(b.h[k][j]);
tempm.h[i][j]%=N;
}
return tempm;
}
M fang(int k,M linshi)
{
if(k==0)return e;
if(k==1)return c;
linshi=fang(k>>1,linshi);
linshi=mul(linshi,linshi);
if(k&1)linshi=mul(linshi,c);
return linshi;
}
int main()
{
int k;
c.h[1][0]=1;c.h[1][1]=0;
e.h[0][0]=e.h[1][1]=1;
e.h[1][0]=e.h[0][1]=0;
while(~scanf("%lld %lld %d %d %d",&a,&b,&c.h[0][0],&c.h[0][1],&k))
{
M res=c;
if(k==0){printf("%lld\n",a%N);continue;}
if(k==1){printf("%lld\n",b%N);continue;}
res=fang(k-1,res);
printf("%lld\n",(res.h[0][0]*b+res.h[0][1]*a)%N);
}
return 0;
}
九度OJ 1081 清华09机试题之递推数列
最新推荐文章于 2017-01-15 15:33:58 发布