1 0101111 10 100000001
1111000 001000010
lcy | We have carefully selected several similar problems for you: 1757 1588 2604 2256 2294
/*
能够把转换关系找出来的话，这题就是个比较裸的快速幂了
用a[i]*代表上一层的a[i],通过推样例会发现如下规律：
a[i]=(a[i]*+a[i-1]*)%2,其中a[1]=(a[1]*a[n]*)%2
第一次敲快速幂，存做模板
*/
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
#include <queue>
#include <map>
#include <stack>
#include <set>
#include <cmath>
#include <vector>
#include<cstdlib>
#define INF 1<<30
#pragma comment (linker,"/STACK:102400000,102400000")
using namespace std;
#define maxn 150
#define moo 2
char s[maxn];
int n,m;
struct Mat
{
int a[maxn][maxn];
};
Mat res1;
Mat muti(Mat a,Mat b)
{
Mat res;
memset(res.a,0,sizeof(res.a));
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
{
for(int k=0;k<n;k++)
res.a[i][j]=(res.a[i][j] + a.a[i][k]*b.a[k][j])%moo;
}
return res;
}
Mat pow(Mat res,int m)
{
Mat ss;
memset(ss.a,0,sizeof(ss.a));
for(int i=0;i<n;i++)
ss.a[i][i]=1;
while(m)
{
if(m%2)
ss=muti(ss,res);
res=muti(res,res);
m/=2;
}
return ss;
}
int main()
{
int ans[maxn];
Mat oo;
while(~scanf("%d",&m))
{
scanf("%s",s);
n=strlen(s);
for(int i=0;i<n;i++)
ans[i]=s[i]-'0';
int ant[maxn];
memset(ant,0,sizeof(ant));
memset(res1.a,0,sizeof(res1.a));
res1.a[0][0]=res1.a[n-1][0]=res1.a[n-1][n-1]=1;
for(int i=0;i<n-1;i++)
res1.a[i][i]=res1.a[i][i+1]=1;
oo=pow(res1,m);
for(int j=0;j<n;j++)
for(int i=0;i<n;i++)
ant[j]=(ant[j]+ans[i]*oo.a[i][j])%moo;
for(int i=0;i<n;i++)
cout<<ant[i];
cout<<endl;
}
return 0;
}
版权声明：本文为博主原创文章，欢迎指教~
原文地址：http://blog.csdn.net/u013712847/article/details/48005245
