NOJI 148 fibonacci数列(二) 矩阵乘法二分幂

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
struct app
{
int p[30][30];
};
app mult(app a, app b)
{
app c;
int i , j, k;
for( i =0; i<2; i++)
for(j =0; j<2; j++)
{
c.p[i][j]=0;
for(k =0; k<2; k++)
{
c.p[i][j]=(c.p[i][j]+a.p[i][k]*b.p[k][j])%10000;
}
}
return c;
}

app pow(app a, int k)
{
app b;
memset(b.p, 0, sizeof(b.p));
for(int i=0;i<2; i++)
{
b.p[i][i]=1;
}
while(k)
{
if(k&1)
b=mult(b,a);
a=mult(a,a);
k>>=1;
}
return b;
}
int main()
{
int n;
while(~scanf("%d", &n)&&n!=-1)
{
app a;
a.p={{1,1},{1,0}};
app d=pow(a, n);
printf("%d\n",d.p[0][1]);
}
return 0;
}