序列问题

https://zybuluo.com/ysner/note/1174224

题面

现有一个长度为(n)的序列({ai}),找出两个非空的集合(S)、(T)。
这两个集合要满足以下的条件:

• 两个集合中的元素都为整数,且都在([1, n])里,即(S_i,T_iin[1, n])。
• 对于集合(S)中任意一个元素(x),集合(T)中任意一个元素(y),满足(x<y)。
• 对于大小分别为(p,q)的集合(S)与(T),满足
  (a[S_1]igoplus a[S_2]igoplus a[S_3] ... igoplus a[s_p] = a[t_1]&a[t_2]&a[t_3] ...&a[t_q])

询问一共有多少对这样的集合((S,T))。

• (60pts) (nleq100)
• (100pts) (nleq1000,a_i<1024)

解析

(nleq100)

直接设(lt[i][j])表前(i)个数中 已取数 异或和 为(j)的方案数。
(rt[i][j])同理。
然后枚举分割点和相等值统计答案。
注意一下开始分割点要放在取了的数上，否则会计重。

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=1024;
int n,a[N];
ll lt[1200][1200],rt[1200][1200][2],ans;
il ll gi() 
{
  re ll x=0,t=1;
  re char ch=getchar();
  while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
  if(ch=='-') t=-1,ch=getchar();
  while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
  return x*t;
}
il void wri(re int x)
{
  if(x<0) putchar('-'),x=-x;
  if(x>9) wri(x/10);
  putchar(x%10+'0');
}
int main()
{
  freopen("sequence.in","r",stdin);
  freopen("sequence.out","w",stdout);
  n=gi();
  fp(i,1,n) a[i]=gi();
  fp(i,1,n) lt[i][a[i]]=1,rt[i][a[i]][1]=1;
  fp(i,1,n)
    fp(j,0,N)
    {
      lt[i][j^a[i]]+=lt[i-1][j];
      lt[i][j]+=lt[i-1][j];
    }
  fq(i,n,1)
    fp(j,0,N)
    {
      rt[i][j&a[i]][1]+=rt[i+1][j][1]+rt[i+1][j][0];
      rt[i][j][0]+=rt[i+1][j][1]+rt[i+1][j][0];
    }
  fp(i,1,n)
    fp(j,0,N)
    ans+=lt[i-1][j]*rt[i][j][1];
  printf("%lld
",ans);
  fclose(stdin);
  fclose(stdout);
  return 0;
}

复杂度(O(1024*n*2))

(nleq1000)

注意到(lt[i][j])的值可能以(3!)的比例呈现。。。
所以我们需要高精。。。
然而我们发现强行高精乘效率很低。
于是我们要把(dp)方程式转化一下。
设(dp[i][j][0/1])表示第(i-n)个数中，已取数"and和"（(0)）或"xor和"(&)"and和"（(1)）为(j)的方案数。
应选择倒推，这样最后的合法情况中(j)一定为(0)。
则
(dp[i][j& a[i]][0]+=dp[i+1][j][0])
(dp[i][j][0]+=dp[i+1][j][0])
(dp[i][jigoplus a[i]][1]+=dp[i+1][j][1]+dp[i+1][j][0])
(dp[i][j][1]+=dp[i+1][j][1])
(注意短式不能与长式合并，因长式中第二位取值受到限制）
答案在(dp[1][0][1])中。
这样只用高精加即可。
同时，如果把(dp)数组转为高精，就需要滚动以省空间。
滚动前：

  n=gi();
  fp(i,1,n) a[i]=gi();
  fp(i,1,n) dp[i][a[i]][0]=1;
  fq(i,n,1)
    fp(j,0,1023)
    {
      dp[i][j&a[i]][0]+=dp[i+1][j][0];
      dp[i][j^a[i]][1]+=dp[i+1][j][1]+dp[i+1][j][0];
      dp[i][j][0]+=dp[i+1][j][0];
      dp[i][j][1]+=dp[i+1][j][1];
    }
  printf("%lld
",dp[1][0][1]);

滚动后：
注意到(dp)数组进行转移前要先集体赋值，否则可能会加上上两次的方案。并且赋值为(1)改为(++)。

 n=gi();
  fp(i,1,n) a[i]=gi();
  fq(i,n,1)
    {
      re int now=i&1,las=(i&1)^1;
      fp(j,0,N-1) dp[now][j][0]=dp[las][j][0],dp[now][j][1]=dp[las][j][1];
      dp[now][a[i]][0]++;
      fp(j,0,N-1)
      {
        dp[now][j&a[i]][0]+=dp[las][j][0];
        dp[now][j^a[i]][1]+=dp[las][j][1]+dp[las][j][0];
    }
    }
  cout<<dp[1][0][1]<<endl;

更加毒瘤

发现出题人很卡复杂度，高精运算复杂度过高还是会(GG)。
于是我们需要压位。
并且如果高精数位数很大速度也很慢！！！
于是压九位，高精数开五十位即可。
注意一下输出方式。

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=1024;
int n,a[N];
class BigInteger
{
public:
  int sz,num[50];
  BigInteger(){sz=1;memset(num,0,sizeof(num));}
  void print()
  {
    fq(i,sz,1) if(i^sz) printf("%09d",num[i]);else printf("%d",num[i]);
    puts("");
  }
}dp[2][1200][2],ysn;
BigInteger operator + (BigInteger A,BigInteger B)
{
  BigInteger Ans;
  re int s=max(A.sz,B.sz);
  Ans.sz=s;
  fp(i,1,s) Ans.num[i]=A.num[i]+B.num[i];
  fp(i,1,s)
    if(Ans.num[i]>999999999)
      {
	Ans.num[i+1]+=Ans.num[i]/1000000000;
	Ans.num[i]=Ans.num[i]%1000000000;
      }
  if(Ans.num[s+1]) ++Ans.sz;
  return Ans;
}
il ll gi() 
{
  re ll x=0,t=1;
  re char ch=getchar();
  while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
  if(ch=='-') t=-1,ch=getchar();
  while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
  return x*t;
}
il void wri(re int x)
{
  if(x<0) putchar('-'),x=-x;
  if(x>9) wri(x/10);
  putchar(x%10+'0');
}
int main()
{
  freopen("sequence.in","r",stdin);
  freopen("sequence.out","w",stdout);
  n=gi();ysn.num[1]=1;ysn.sz=1;
  fp(i,1,n) a[i]=gi();
  fq(i,n,1)
    {
      re int now=i&1,las=(i&1)^1;
      fp(j,0,N-1) dp[now][j][0]=dp[las][j][0],dp[now][j][1]=dp[las][j][1];
      dp[now][a[i]][0]=dp[now][a[i]][0]+ysn;
      fp(j,0,N-1)
      {
        dp[now][j&a[i]][0]=dp[now][j&a[i]][0]+dp[las][j][0];
        dp[now][j^a[i]][1]=dp[now][j^a[i]][1]+dp[las][j][1]+dp[las][j][0];
      }
    }
  dp[1][0][1].print();
  fclose(stdin);
  fclose(stdout);
  return 0;
  }