题目链接:http://codeforces.com/problemset/problem/901/B
多项式为n次,gcd的次数又必须是n从中可以发现gcd每次只能将1阶,所以构造方程f(x)=q(x)*x+r(x)。可以看出这是个递推,直接从后面开始往前面递推,最后的gcd状态b一定为0,所以最后状态表示为gcd(a,0),又因为方程系数的绝对值小于等于1所以a为1,于是列出末状态dp[1][1]=1,dp[2][1]=1,dp[2][2]=0。dp表示的是多项式A递推到n+1那么dp[n+1]表示A的各项系数dp[n]则为B的各项系数,递推的方程为:dp[i][j]=j==1?dp[i-1][j]:(dp[i-1][j]+dp[i-2][j-2])%2在这里可以直接mod2(因为此式子相当于f(x)=q(x)*x+r(x)两边系数mod2等式不变)
代码:
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
const int maxn=200;
int dp[maxn][maxn]={};
int main()
{
int n;
while(scanf("%d",&n)==1){
memset(dp,0,sizeof(dp));
dp[1][1]=1;
dp[2][1]=1;
for(int i=3;i<=n+1;i++)
for(int j=1;j<=i;j++){
if(j==1) dp[i][j]=dp[i-1][j];
else dp[i][j]=(dp[i-1][j]+dp[i-2][j-2])%2;
}
printf("%d\n",n);
for(int i=n+1;i>1;i--) printf("%d ",dp[n+1][i]);
printf("%d\n",dp[n+1][1]);
printf("%d\n",n-1);
for(int i=n;i>1;i--) printf("%d ",dp[n][i]);
printf("%d\n",dp[n][1]);
}
return 0;
}