码迷,mamicode.com
首页 > 其他好文 > 详细

POJ 3070

时间:2014-08-30 21:41:50      阅读:268      评论:0      收藏:0      [点我收藏+]

标签:blog   os   io   for   div   log   sp   amp   c   

水题,结果MOD 10000

#include <iostream>
#include <cstdio>
using namespace std;

int n=2,M=10000;
struct Matrax {
    int m[35][35];
};
 
Matrax a,per;
 
void initial(){
    int i,j;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            per.m[i][j]=(i==j);
        }
    }
    a.m[0][0]=a.m[0][1]=a.m[1][0]=1;
    a.m[1][1]=0;
}
 
Matrax multi(Matrax a,Matrax b){
    Matrax c;
    int k,i,j;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            c.m[i][j]=0;
            for(k=0;k<n;k++){
                c.m[i][j]=(c.m[i][j]+a.m[i][k]*b.m[k][j])%M;
            }
            c.m[i][j]%=M;
        }
    }
    return c;
}
 
Matrax Power(int k){
    Matrax c,p,ans=per;
    p=a;
    while(k){
        if(k&1){
            ans=multi(ans,p);
            k--;
        }
        else {
            k/=2;
            p=multi(p,p);
        }
    }
    return ans;
}

int main(){
	int k;
	initial();
	while(scanf("%d",&k)!=EOF){
		if(k==-1) break;
		Matrax ans=Power(k);
		printf("%d\n",ans.m[0][1]);
	}
	return 0;
}

  

POJ 3070

标签:blog   os   io   for   div   log   sp   amp   c   

原文地址:http://www.cnblogs.com/jie-dcai/p/3946854.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!