CF1364E X-OR

题目链接：CF1364E

有一个固定的长度为 n 的排列 P，其值域为 [0,n-1]，你可以进行不超过 4269 次询问，之后你需要输出这个排列 P​。

你可以按照 ? a b 的格式进行询问，之后你会得到 P_a 与 P_b 的按位或。

保证 3\le n\le 2048，0\le P_i\le n-1。

Tutorial

首先发现如果找到 0 的位置，其他的全部问一遍就都出来了，因此转化为如何找 0。

不妨先随机两个位置 x,y，再枚举其他所有 i：

• P_x P_y > P_y P_iP_x 必然非 0，因此用 i 代替 x 继续循环。
• P_xP_y = P_yP_i：P_y 必然非 0，因此用 i 代替 y 继续循环。
• P_xP_y <P_yP_iP_i 必然非 0，因此不换。

那么最后肯定剩下 x,y，其中一个必为 0。

于是我们可以再随机一个数 i：

• P_x P_i > P_yP_iP_x 必然非 0，因此 P_y=0。
• P_xP_i < P_yP_iP_y 必然非 0，因此 P_x=0。
• P_xP_i = P_yP_i：分辨不出来啥再重新随一个:(

---

简单分析一下询问次数：

首先第一部分找到两个中必有 0 的一对数字，扫了一次所有的点，而其中第二种情况需要多询问一次新的 P_xP_i的值（为下次循环做准备），而产生第二种情况的充要条件是 P_x \& P_y = P_x,P_i\& P_y = P_i\sum_{i=1}^n\frac 1{n^3} 2^{2\operatorname{popcount}(i)}\approx 0.0056843422353267669677734375，可以看出这是非常小的。

因此这部分的询问次数大约为 0.0056843422353267669677734375\times 2048+2048\approx 2059。

其次第二部分随到第三种情况的概率：\sum_{i=1}^n \frac 1{n^2}2^{\operatorname{popcount}(i)}\approx 0.04223537445068359375。

不妨设其期望次数为 f，则有 f=0.04223537445068359375\times (f+1) +(1-0.04223537445068359375)\times 1

因此 f=1.0441。

所以总询问次数期望为：2059+1+2048=4108，远远小于 4269。

Solution

#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define W while
#define I inline
#define RI register int
#define LL long long
#define Cn const
#define CI Cn int&
#define gc getchar
#define D isdigit(c=gc())
#define pc(c) putchar((c))
#define min(x,y) ((x)<(y)?(x):(y))
#define max(x,y) ((x)>(y)?(x):(y))
using namespace std;
namespace Debug{
    Tp I void _debug(Cn char* f,Ty t){cerr<<f<<'='<<t<<endl;}
    Ts I void _debug(Cn char* f,Ty x,Ar... y){W(*f!=',') cerr<<*f++;cerr<<'='<<x<<",";_debug(f+1,y...);}
    Tp ostream& operator<<(ostream& os,Cn vector<Ty>& V){os<<"[";for(Cn auto& vv:V) os<<vv<<",";os<<"]";return os;}
    #define gdb(...) _debug(#__VA_ARGS__,__VA_ARGS__)
}using namespace Debug;
namespace FastIO{
    Tp I void read(Ty& x){char c;int f=1;x=0;W(!D) f=c^'-'?1:-1;W(x=(x<<3)+(x<<1)+(c&15),D);x*=f;}
    Ts I void read(Ty& x,Ar&... y){read(x),read(y...);}
    Tp I void write(Ty x){x<0&&(pc('-'),x=-x,0),x<10?(pc(x+'0'),0):(write(x/10),pc(x%10+'0'),0);}
    Tp I void writeln(Cn Ty& x){write(x),pc('\n');}
}using namespace FastIO;
Cn int N=2e5+10;
int n,Ans[N];
I int Q(CI x,CI y){
    RI t;cout<<"? "<<x<<" "<<y<<endl;
    cin>>t;return t;
}
int main(){
    mt19937 mt(time(nullptr));
    RI i,x,y,v,t,q1,q2,p0;read(n);x=mt()%n+1,y=mt()%n+1;W(x==y) y=mt()%n+1;
    for(v=Q(x,y),i=1;i<=n;i++) if(i^x&&i^y){
        if(v>(t=Q(y,i))) x=i,v=t;
        else if(v==t) y=i,v=Q(x,y);
    }do{
    t=mt()%n+1;W(t==xt==y) t=mt()%n+1;
    if((q1=Q(x,t))>(q2=Q(y,t))) p0=y;else if(q1<q2) p0=x;}while(!p0);
    for(i=1;i<=n;i++) Ans[i]=i==p0?0:Q(p0,i);
    for(cout<<"! ",i=1;i<=n;i++) cout<<Ans[i]<<" ";cout<<endl;
    return 0;
}

原文地址：https://cloud.tencent.com/developer/article/2111506