矩阵快速幂
struct Matrix
{
LL m[66][66];
};
Matrix operator * (Matrix a,Matrix b)
{
Matrix c;
for(int i=0; i<M; i++)
for(int j=0; j<M; j++)
c.m[i][j]= 0;
for(int k=0; k<M; k++)
for(int i=0; i<M; i++)
{
if(a.m[i][k] <= 0) continue;
for(int j=0; j<M; j++)
{
c.m[i][j]=(c.m[i][j]+a.m[i][k]*b.m[k][j]+MOD)%MOD;
}
}
return c;
}
Matrix operator ^ (Matrix a,LL b)
{
Matrix c;
for(int i=0; i<n; i++)
for(int j=0; j<n; j++)
c.m[i][j] = c.m[i+n][j+n] = ( i == j );
while(b)
{
if(b&1) c= c * a ;
b >>= 1;
a = a * a ;
}
return c;
}
Matrix operator + (Matrix a,Matrix b)
{
for(int i=0; i<M; i++)
for(int j=0; j<M; j++)
a.m[i][j]=(a.m[i][j]+b.m[i][j]+MOD)%MOD;
return a;
}
博弈SG函数
(一) 记忆化搜索dfs版
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int SG[10100],n,m,s[102],k;
int dfs(int x)
{
if(SG[x]!=-1) return SG[x];
bool vis[110];
memset(vis,0,sizeof(vis));
for(int i=0;i<k;i++)
{
if(x>=s[i])
{
dfs(x-s[i]);
vis[SG[x-s[i]]]=1;
}
}
int e;
for(int i=0;;i++)
if(!vis[i])
{
e=i;
break;
}
return SG[x]=e;
}
int main()
{
int cas,i;
while(scanf("%d",&k)!=EOF)
{
if(!k) break;
memset(SG,-1,sizeof(SG));
for(i=0;i<k;i++) scanf("%d",&s[i]);
sort(s,s+k);
scanf("%d",&cas);
while(cas--)
{
int t,sum=0;
scanf("%d",&t);
while(t--)
{
int num;
scanf("%d",&num);
sum^=dfs(num);
}
if(sum==0) printf("L");
else printf("W");
}
printf("\n");
}
return 0;
}
(二) 预处理打表版
#include<stdio.h>
#include<string.h>
#include <string>
#include <iostream>
using namespace std;
const int N = 10008;
int s[108],t;
int sg[N];
bool has[N];
void sg_solve(int *s,int t,int N)
{
int i,j;
memset(sg,0,sizeof(sg));
for(i=1;i<=N;i++)
{
memset(has,0,sizeof(has));
for(j=0;j<t;j++)
if(i - s[j] >= 0)
has[sg[i-s[j]]] = 1;
for(j=0;j<=N;j++)
if(!has[j])
break;
sg[i] = j;
}
return ;
}
int main()
{
int i,j,n,m,h;
while(scanf("%d",&t),t)
{
string ans="";
for(i=0;i<t;i++)
scanf("%d",&s[i]);
sg_solve(s,t,N);
scanf("%d",&n);
for(i=0;i<n;i++)
{
scanf("%d",&m);
int res = 0;
for(j=0;j<m;j++)
{
scanf("%d",&h);
res ^= sg[h];
}
ans+=res?'W':'L';
}
cout<<endl;
}
return 0;
}