省队集训Round3 DAY6

T1

题解

这道题应该是可以通过组合数直接计算的，但是我不会数学证明，所以就用了一种简单粗暴的方式。
$ans=\sum_{2|i} C(n,i)*C(n,n-i)$
$ans=\sum_{2|i} C(n,i)^2$
$10^18$肯定不能直接枚举，如果要计算的话也需要用到lucas定理。观察发现这题的模数比较小，所以从模数入手。
考虑Lucas定理。
$n=a[0]*p^0+a[1]*p^1+a[2]*p^2+...$
$m=b[0]*p^0+b[1]*p^1+b[2]*p^2+...$
对n,m进行p进制分解，那$C(n,m) \% p=C(a[0],b[0])*C(a[1],b[1])*....\% p$
如果某一位$a[i]<b[i]$，那么答案在模p意义下一定是0，也就是没有贡献，所以我们现在考虑计算出所有有贡献的答案。
确定了a之后，我们就可以通过枚举计算每一位的选择。因为i必须是偶数，所以我们可以进行数位dp.
$f[i][0/1]$表示确定了第i位后，m是偶数/奇数的总贡献。
转移过程直接看代码吧。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define p 1000003
#define LL long long 
using namespace std;
LL jc[2000000],inv[2000000],f[10][10],mi[10];
LL n;
int a[10],cnt;
LL quickpow(LL num,int x)
{
    LL ans=1,base=num%p;
    while (x) {
        if (x&1) ans=ans*base%p;
        x>>=1;
        base=base*base%p;
    }
    return ans;
}
LL calc(int n,int m)
{
    return jc[n]*inv[m]%p*inv[n-m]%p;
}
LL pow(LL x){
    return x*x%p;
}
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    cin>>n;
    jc[0]=1;
    for (int i=1;i<=p;i++) jc[i]=(LL)jc[i-1]*i%p;
    for (int i=0;i<=p;i++) inv[i]=quickpow(jc[i],p-2);
    LL x=n; 
    while (x) {
        a[++cnt]=(LL)x%p;
        x/=p;
    }
    mi[0]=1;
    for (int i=1;i<=cnt;i++) mi[i]=mi[i-1]*mi[i-1]%p;
    for (int i=0;i<=a[1];i++)
     (f[1][i*mi[0]%2]+=pow(calc(a[1],i)))%=p;
    for (int i=2;i<=cnt;i++)
     for (int j=0;j<=a[i];j++){
        LL t=j*mi[i-1]%2;
        (f[i][t]+=f[i-1][0]*pow(calc(a[i],j))%p)%=p;
        (f[i][(t+1)%2]+=f[i-1][1]*pow(calc(a[i],j))%p)%=p;
     }
    cout<<f[cnt][0]<<endl;
}

T3

题解

可持久化trie
首先我们考虑如果只有异或操作该怎么做。因为是对全局取xor，所以我们可以把要异或的值记录下来，每次查询的时候考虑上异或值即可，不会对trie的形态造成影响。
然后考虑and，对于每一位分开考虑，如果是and 1的话没有影响。关键就是and 0,会使所有数当前位的取值变成0。
or的话，如果是or 1,会使所有数当前位变成1。我们不妨对于每一位开一个标记数组，表示当前位是否全部相同以及值是多少。如果某一位在操作的过程中由值不全相同变成相同，那么我们就重构可持久化trie，如果某一位全部相同，我们就把所有数同一放到左儿子，方便查询。
因为每一位只有一次由不同变成相同的机会，所以重构的次数不会太多，最多只有30次。
注意异或的时候别忘了修改每一位的标记。

代码

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
#define N 100003
using namespace std;
int xornum,n,m,a[N],ch[N*20][2],size[N*20],pd[N],opt[N],root[N],sz;
void insert(int i,int x)
{
    int pre=root[i-1];
    root[i]=++sz;
    int now=root[i];
    for (int i=30;i>=0;i--){
        int t=(x>>i)&1;
        if (pd[i]) t=0;
        size[now]=size[pre]+1;
        ch[now][t^1]=ch[pre][t^1];
        ch[now][t]=++sz;
        now=ch[now][t]; pre=ch[pre][t];
    }
    size[now]=size[pre]+1;
}
int query(int i,int j,int x)
{
    int ans=0;
    int R=i; int L=j;
    for (int t=30;t>=0;t--) {
        if (pd[t]) ans+=(1<<t)*opt[t],R=ch[R][0],L=ch[L][0];
        else {
            int son=(xornum>>t)&1;
            int k=size[ch[R][son]]-size[ch[L][son]];
            if (k>=x) {
                ans+=(1<<t)*0;
                R=ch[R][son],L=ch[L][son];
            }
            else {
                x-=k;
                ans+=(1<<t)*1;
                R=ch[R][son^1],L=ch[L][son^1];
            }
        }
    } 
    return ans;
}
int main()
{
    freopen("c.in","r",stdin);
    freopen("c.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]);
    insert(1,0);
    for (int i=1;i<=n;i++) 
     insert(i+1,a[i]);
    for (int i=1;i<=m;i++) {
        char s[10]; scanf("%s",s+1);
        int l,r,x;
        if (s[1]=='X') {
            scanf("%d",&x);
            xornum^=x;
            for (int j=0;j<=30;j++)
             if (pd[j]&&((x>>j)&1)) opt[j]^=1;
        }
        if (s[1]=='O'){
            scanf("%d",&x); bool flag=false;
            for (int j=0;j<=30;j++)
             if ((x>>j)&1) {
                if (pd[j]) opt[j]=1;
                else{
                    flag=true;
                    pd[j]=opt[j]=1;
                }
             }
             if (flag) {
                sz=0;
                memset(size,0,sizeof(size));
                memset(ch,0,sizeof(ch));
                memset(root,0,sizeof(root));
                insert(1,0);
                for (int j=1;j<=n;j++) insert(j+1,a[j]);
             }
        }
        if (s[2]=='n'){
            scanf("%d",&x); bool flag=false;
            for (int j=0;j<=30;j++)
             if (((x>>j)&1)==0){
                if (pd[j]) opt[j]=0;
                else{
                    flag=true;
                    pd[j]=1;opt[j]=0;
                }
             }
            if (flag){
                sz=0;
                memset(size,0,sizeof(size));
                memset(ch,0,sizeof(ch));
                memset(root,0,sizeof(root));
                insert(1,0);
                for (int j=1;j<=n;j++) insert(j+1,a[j]);
            }
        }
        if (s[2]=='s') {
            scanf("%d%d%d",&l,&r,&x);
            printf("%d\n",query(root[r+1],root[l],x));
        }
    }
}