[loj502]ZQC的截图

给每一个人一个随机数$R_{i}$，将一个消息中所有人的的$R_{i}$在三进制下相加（多次出现需要多个$R_{i}$），最终之和若为0，即判定答案为-1，若为某个$R_{i}$或$R_{i}+R_{i}$（三进制下），则为$i$，否则为-2

显然这一做法是随机的，但其每一次失败都意味者结果为0，但存在若干个$R_{i}$出现次数不为3的倍数，且其之和恰为0（关于结果是$R_{i}$或$R_{i}+R_{i}$可以看作减去$R_{i}$或$R_{i}+R_{i}$）

注意到$R_{i}$是随机的，因此若干个$R_{i}$之和也是随机的，为0的概率恰为$\frac{1}{V}$（其中$V$为值域，应为3的幂次）

但事实上，这里还有两个问题：

1.如何快速判定是否为$R_{i}$或$R_{i}+R_{i}$，由于有$o(n)$个值，且空间仅有128MB，哈希范围仅能为$10^{7}$，因此需要在哈希的基础上，每一个位置再开一个vector（此时哈希范围仅能为$10^{6}$）

2.三进制下暴力计算是$\log_{3}V$的，选择$V=3^{35}$，并预处理出$3^{7}$以内任意两数加法，那么单次加法的复杂度降为$o(5)$（这些都只是常数优化，实际理论复杂度还是$o(\log_{3}V)$）

#include<bits/stdc++.h>
using namespace std;
#define N 1000005
#define K 2187
#define mod 1000003
#define ll long long
vector<int>mat[mod];
int n,m,x,y,ans,init[K][K];
ll R[N],dep[N<<1];
int read(){
    int x=0,flag=0;
    char c=getchar();
    while ((c<'0')||(c>'9')){
        if (c=='-')flag=1;
        c=getchar();
    }
    while ((c>='0')&&(c<='9')){
        x=x*10+c-'0';
        c=getchar();
    }
    if (flag)x=-x;
    return x; 
}
void write(int x){
    if (x<0){
        putchar('-');
        x=-x;
    }
    int num[11];
    num[0]=0;
    while (x){
        num[++num[0]]=x%10;
        x/=10;
    }
    if (!num[0])putchar('0');
    for(int i=num[0];i;i--)putchar(num[i]+'0');
    putchar('\n');
}
int add_low(int x,int y){
    int s=1,ans=0;
    for(int i=0;i<7;i++){
        ans=ans+(x+y)%3*s;
        s*=3,x/=3,y/=3;
    }
    return ans;
}
ll add_high(ll x,ll y){
    ll s=1,ans=0;
    for(int i=0;i<5;i++){
        ans=ans+init[x%K][y%K]*s;
        s*=K,x/=K,y/=K;
    }
    return ans;
}
int main(){
    srand(time(0));
    for(int i=0;i<K;i++)
        for(int j=0;j<K;j++)init[i][j]=add_low(i,j);
    n=read(),m=read();
    for(int i=1;i<=n;i++){
        for(int j=0;j<35;j++)R[i]=R[i]*3+rand()%3;
        mat[R[i]%mod].push_back(i);
        mat[add_high(R[i],R[i])%mod].push_back(i);
    }
    for(int i=1;i<=m;i++){
        x=read(),y=read();
        x^=ans,y^=ans;
        dep[i]=add_high(dep[y],R[x]);
        if (!dep[i])ans=-1;
        else{
            ans=-2;
            for(int j=0;j<mat[dep[i]%mod].size();j++){
                ll x=R[mat[dep[i]%mod][j]];
                if ((x==dep[i])||(add_high(x,x)==dep[i])){
                    ans=mat[dep[i]%mod][j];
                    break;
                }
            }
        }
        write(ans);
    }
    return 0;
}