superoj887 set

题目描述：

输入格式

第一行一个数：N 。
接下来一行 N 个数，其中第 i 数 ai 表示初始时刻第 i 堆石子中的石子数。

输出格式

输出一行一个数，即答案。

样例数据 1

输入　 [复制]

2 
1 2

输出

4

备注

【样例说明】
四种状态分别是：
{1} , {}
{2} , {}
{1,2} , {}
{2} , {1}

对于 100 分的数据： 0≤N,ai≤1000 。

分析：

        因为这道题 A赢B的局面==B赢A的局面 所以

 ans=（3^n-dp[相等情况]）/2 mod p

f[石子数][异或差]

f[i][j^a[i]]=(f[i][j^a[i]]+f[i-1][j])%mm;给A

f[i][j^a[i]]=(f[i][j^a[i]]+f[i-1][j])%mm;给B

f[i][j]=(f[i][j]+f[i-1][j])%mm;不给

f[n][0]为dp[相等情况]

最后mod p 用逆元做

#include<ctime>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
long long a[5001];
long long n;
long long f[4001][4096];
long long mm=1998585857;
long long read()
{
	long long k=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0') {if(c=='-') f=-1; c=getchar();}
	while(c>='0'&&c<='9') {k=k*10+(c-'0'); c=getchar();}
	return k*f;
}

long long ksm(long long a,long long b,long long c)
{
	long long ans=1;
	
	while(b>0)
	 {
	 	if(b%2==1) ans=(ans*a)%c;
	 	b=b/2;
	 	a=(a*a)%c;
	 }
	return ans;
}

int main()
{
	//..freopen("set.in","r",stdin);
	//freopen("set.out","w",stdout);
	long long i,j,k;
	n=read();
	for(i=1;i<=n;i++)
	  {
	  	 a[i]=read();
	  }
  
	f[0][0]=1;
	for(i=1;i<=n;i++)
	  {
	  	for(j=0;j<=4095;j++)
	  	 { 
	  	   f[i][j^a[i]]=(f[i][j^a[i]]+f[i-1][j])%mm;
	  	   f[i][j^a[i]]=(f[i][j^a[i]]+f[i-1][j])%mm;
	  	   
	  	   f[i][j]=(f[i][j]+f[i-1][j])%mm;
	     }
	  }
	long long aa=ksm(3,n,mm);
	aa=aa-f[n][0];
	if(aa<0) aa+=mm;
	long long bb=ksm(2,mm-2,mm);
	
	cout<<(aa*bb)%mm;
	return 0;
}