http://poj.org/problem?id=3070
#include <iostream>
#include <string.h>
#include <stdlib.h>
#include <cstdio>
#include <algorithm>
#define mod 10000
using namespace std;
struct m
{
int a[3][3];
} init,res;
int n;
m Mult(m x,m y)
{
m tmp;
for(int i=0; i<2; i++)
{
for(int j=0; j<2; j++)
{
tmp.a[i][j]=0;
for(int k=0; k<2; k++)
{
tmp.a[i][j]=(tmp.a[i][j]+x.a[i][k]*y.a[k][j])%mod;
}
}
}
return tmp;
}
m Pow(m x,int n)
{
m tmp;
for(int i=0; i<2; i++)//将矩阵tmp初始化成单位矩阵[1,0]
{
for(int j=0; j<2; j++)// [0,1]
{
tmp.a[i][j]=(i==j);
}
}
while(n)
{
if(n&1) tmp=Mult(tmp,x);
n>>=1;
x=Mult(x,x);
}
return tmp;
}
int main()
{
while(scanf("%d",&n)!=EOF)
{
if(n==-1) break;
init.a[0][0]=1;
init.a[0][1]=1;
init.a[1][0]=1;
init.a[1][1]=0;
res=Pow(init,n);
cout<<res.a[1][0]<<endl;
}
return 0;
}
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.
- 25.
- 26.
- 27.
- 28.
- 29.
- 30.
- 31.
- 32.
- 33.
- 34.
- 35.
- 36.
- 37.
- 38.
- 39.
- 40.
- 41.
- 42.
- 43.
- 44.
- 45.
- 46.
- 47.
- 48.
- 49.
- 50.
- 51.
- 52.
- 53.
- 54.
- 55.
- 56.
- 57.
- 58.
- 59.
- 60.