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]任意个数and值的方案数
//所以要倒着计算
void initAnd()
{
fAnd
[a
]=1;
memset(visited,0,sizeof(visited));

sAnd.clear();
sAnd.insert(a
);
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]任意个数and值的方案数
//所以要倒着计算
void initAnd()
{
vector<int> p;

p.push_back(1);
fAnd
[a
].push_back(1);

sAnd.clear();
sAnd.insert(a
);
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;
}