superoj887 set

題目描述:

輸入格式

第一行一個數:N 。
接下來一行 N 個數,其中第 i 數 ai 表示初始時刻第 i 堆石子中的石子數。

輸出格式

輸出一行一個數,即答案。

樣例數據 1

輸入  [複製]


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;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章