有一个猜奖者和一个主持人，一共有 \(n\) 扇门，只有一扇门后面有奖，主持人事先知道哪扇门后有奖，而猜奖者不知道。每一轮，猜奖者选择它认为的有奖概率最大（如果有多个最大，随机选一个）的一扇门，主持人从剩下的且门后没有奖的门中随机打开一扇。直到剩两扇门时，猜奖者做出的选择就是他最后的选择。
现在由你来安排主持人每次打开哪一扇门，猜奖者不知道有内幕，他还认为主持人是从可以打开的门中随机一扇打开。你要使猜奖者获奖概率最低，求这个概率。
（Discover Probability，你的快乐老家 ）
Solution
这个题真是有趣，一步一步来
当然还是先用个 EXCRT 的板子把外面的数论模型解掉
选手是怎么计算每扇门的概率的？
设现在还有 \(n\) 扇门，第 \(i\) 扇当前的概率是 \(p[i]\) ，不妨设选手选了 \(x\)，主持人打开了 \(y\)
根据样本空间的对称性，\(p[x]\) 保持不变，而 \(p[y]\) 变化为 \(0\)，于是对任意 \(z\neq x,y\)，有 \(p'[z]=p[z] \frac{1-p[x]}{1-p[x]-p[y]}\)
我们顺便发现，由于每次选手选择的门都是当前概率最小的，并且在一波操作后它的概率没变而其它门的概率变大，所以这扇门永远是概率最小的
问题转化
有 \(n\) 个球站队，队分成两段，靠近队头的那段叫做 head，一开始所有的球都在 head 里，所有球中有一个是 good 球，其它球是 bad 球
每次从 head 随机取一个球 \(x\)，从剩下的所有 bad 球中去掉一个，并且把 \(x\) 放到队尾。被追加到队尾的球不再是 head 段中的球
当只剩下两个球的时候，如果 head 球是 good，那么就 WIN，否则就 LOSE
小范围的 DP
设 \(f[i][j][k]\) 表示还剩下 \(i\) 个球，head 中有 \(j\) 个球，good 是队列中第 \(k\) 个球，从这个状态开始玩的最小概率，暴力转移即可，时间复杂度 \(O(n^3)\)
打表
将不同 \(n\) 时的答案打印出来，很容易发现在 \(n>10\) 的时候，选手就必败了
那还等什么，写个暴搜打表交
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N = 205;
string s[12]={"error","error","0.500000","0.666667","0.625000","0.466667",
"0.416667","0.342857","0.291667","0.253968","0.225000","0.000000"};
// Input: n,ai[],bi[]
// Method: solve()
// Output: (returned)
namespace excrt {
const int maxn=100010;
int n;
int ai[maxn],bi[maxn]; //x=a%b
int mul(int a,int b,int mod){
int res=0;
while(b>0){
if(b&1) res=(res+a)%mod;
a=(a+a)%mod;
b>>=1;
}
return res;
}
int exgcd(int a,int b,int &x,int &y){
if(b==0){x=1;y=0;return a;}
int gcd=exgcd(b,a%b,x,y);
int tp=x;
x=y; y=tp-a/b*y;
return gcd;
}
int solve(){
int x,y,k;
int M=bi[1],ans=ai[1];
for(int i=2;i<=n;i++){
int a=M,b=bi[i],c=(ai[i]-ans%b+b)%b;
int gcd=exgcd(a,b,x,y),bg=b/gcd;
if(c%gcd!=0) return -1;
x=mul(x,c/gcd,bg);
ans+=x*M;
M*=bg;
ans=(ans%M+M)%M;
}
return (ans%M+M)%M;
}
}
signed main(){
int n;
cin>>n;
excrt::n=n;
for(int i=1;i<=n;++i) cin>>excrt::ai[i]>>excrt::bi[i]; //x=a%b
n=excrt::solve();
cout<<s[min(max(0ll,n),11ll)];
}