GFOJ686 序列操作（分段DP，未高精）

数据会举出，还需要有高精度。

#include<cstdio>
#include<cstring>
#include<set>
using namespace std;

const int maxSize=1000;
const int maxValue=1024;
int n,a[maxSize+5],i,j,k,m,ans;
int fAnd[maxSize+5][maxValue],fXor[maxSize+5][maxValue],visited[maxSize+5];
set <int> sAnd,sXor;
set <int>::iterator iter,iter1;

void output(int f[][maxValue])
{
	int i,j;
	
	for(i=1;i<=n;i++)
	{
		printf("%d: ",i);
		for(j=0;j<16;j++)
			printf("%d ",f[i][j]);
		printf("\n");
	}
}

//求出选任意个数xor值的方案数
//f[2]:求出[2..n]任意个数xor值的方案数 
//f[3]:求出[3..n]任意个数xor值的方案数
//f[4]:求出[4..n]任意个数xor值的方案数
//
//f[n]:求出[n..n]任意个数and值的方案数
//所以要倒着计算 
void initAnd()
{
	fAnd[n][a[n]]=1;
	memset(visited,0,sizeof(visited));
	
	sAnd.clear();
	sAnd.insert(a[n]);
	for(i=n-1;i>=2;i--)
	{
		//printf("And %d:",i);
		for(iter=sAnd.begin();iter!=sAnd.end();iter++)
		{
			j=*iter;
			fAnd[i][j]+=fAnd[i+1][j];
			//printf("%d ",j);
			fAnd[i][j & a[i]] += fAnd[i+1][j];
			sAnd.insert(j&a[i]);
		}	
		sAnd.insert(a[i]);
		fAnd[i][a[i]]++;
		//printf("\n");
	}
}


void output1()
{
	freopen("b.txt","w",stdout);
	printf("            Xor  And\n");
	for(j=1;j<=n;j++)
	for(i=0;i<maxValue;i++)	
		if(fXor[j][i]+fAnd[j][i])	printf("%3d %6d %3d %3d\n",j,i,fXor[j][i],fAnd[j][i]);
}


//求出选任意个数xor值的方案数 
void initXor()
{
	int i,j;
	
	fXor[1][a[1]]=1;
	sXor.clear();
	sXor.insert(a[1]);
	for(i=2;i<=n-1;i++)
	{
		for(iter=sXor.begin();iter!=sXor.end();iter++)
		{
			j=*iter;
			//fXor[i][j]+=fXor[i-1][j];
			
			fXor[i][j ^ a[i]] += fXor[i-1][j];
			sXor.insert(j^a[i]);
		}	
		sXor.insert(a[i]);
		fXor[i][a[i]]++;
	}	
}


void solve()
{
	ans=0;
	for(i=1;i<=n-1;i++)		//分成两部分：[1..i],[i+1..n] 
		for(j=0;j<1024;j++)
			{
				ans = fXor[i][j]*fAnd[i+1][j];
			}
	printf("%d\n",ans);
}

int main()
{
	freopen("a.txt","r",stdin) ;
	
	scanf("%d",&n);
	for(i=1;i<=n;i++) scanf("%d",&a[i]);
	
	initAnd();initXor();
	//printf("Xor:\n");output(fXor);printf("\n\nAnd:\n");output(fAnd);
	output1();
	solve();
	return 0;
}

以下增加高精度，未完成。

#include<cstdio>
#include<cstring>
#include<set>
#include<vector>
using namespace std;

const int maxSize=1000;
const int maxValue=1024;
const int mo=10000000;
 
int n,a[maxSize+5],i,j,k,m;

vector<int> fAnd[maxSize+5][maxValue],fXor[2][maxValue]; 

vector<int> ans;
set <int> sAnd,sXor;
set <int>::iterator iter,iter1;

void output(int f[][maxValue])
{
	int i,j;
	
	for(i=1;i<=n;i++)
	{
		printf("%d: ",i);
		for(j=0;j<16;j++)
			printf("%d ",f[i][j]);
		printf("\n");
	}
}

//高精度加法 p1=p1+p2 
//注意长度可能不一样 
void plusBig(vector<int> &p1,vector<int> p2)
{
	int size=max(p1.size(),p2.size());
	int x,i,c=0;
	
	while(p1.size()<size)		p1.push_back(0);
	while(p2.size()<size) p2.push_back(0);
	for(i=0;i<size;i++)
	{
		x=p1[i]+p2[i]+c;
		p1[i]=x%mo;
		c=x/mo;
	}
	
}

//求出选任意个数xor值的方案数
//f[2]:求出[2..n]任意个数xor值的方案数 
//f[3]:求出[3..n]任意个数xor值的方案数
//f[4]:求出[4..n]任意个数xor值的方案数
//
//f[n]:求出[n..n]任意个数and值的方案数
//所以要倒着计算 
void initAnd()
{
	vector<int> p;
	
	p.push_back(1);
	fAnd[n][a[n]].push_back(1);
	
	sAnd.clear();
	sAnd.insert(a[n]);
	for(i=n-1;i>=2;i--)
	{
		//printf("And %d:",i);
		for(iter=sAnd.begin();iter!=sAnd.end();iter++)
		{
			j=*iter;
			plusBig(fAnd[i][j],fAnd[i+1][j]);	//			fAnd[i][j]+=fAnd[i+1][j];
			
			plusBig(fAnd[i][j&a[i]],fAnd[i+1][j]);			//fAnd[i][j & a[i]] += fAnd[i+1][j];
			sAnd.insert(j&a[i]);
		}	
		sAnd.insert(a[i]);
		plusBig(fAnd[i][a[i]],p);
		//printf("\n");
	}
}


void output1()
{
//	freopen("b.txt","w",stdout);
//	printf("            Xor  And\n");
//	for(j=1;j<=n;j++)
//	for(i=0;i<maxValue;i++)	
//		if(fXor[j][i]+fAnd[j][i])	printf("%3d %6d %3d %3d\n",j,i,fXor[j][i],fAnd[j][i]);
}


//求出选任意个数xor值的方案数 
void initXor()
{
	int i,j,now=0;
	vector<int> p;
	
	p.push_back(1);
	
	fXor[0][a[1]].push_back(1);
	sXor.clear();
	sXor.insert(a[1]);
	
	plusBig(ans,fAnd[2][a[1]]);
	
	for(i=2;i<=n-1;i++)
	{
		now=1-now;
		for(iter=sXor.begin();iter!=sXor.end();iter++)
		{
			j=*iter;
			//fXor[i][j]+=fXor[i-1][j];
			
			plusBig(fXor[now][j ^ a[i]],fXor[1-now][j]);//fXor[now][j ^ a[i]] += fXor[1-now][j];
			sXor.insert(j^a[i]);
		}	
		sXor.insert(a[i]);
		plusBig(fXor[now][a[i]],p); //fXor[now][a[i]]++;
		
		
	}	
}



int main()
{
	freopen("a.txt","r",stdin) ;
	
	scanf("%d",&n);
	for(i=1;i<=n;i++) scanf("%d",&a[i]);
	
	initAnd();initXor();
	//printf("Xor:\n");output(fXor);printf("\n\nAnd:\n");output(fAnd);
	//output1();
	printf("%d\n",ans.size());
	//for(i=ans.size()-1;i;i--)		printf("%08d",ans[i]);
	
	
	return 0;
}