CodeChef Sereja and Game [DP 概率 博弈论]

https://www.codechef.com/problems/SEAGM

题意：

n个数（可能存在相同的数），双方轮流取数。如果在一方选取之后，所有
已选取数字的GCD变为1，则此方输。
问：
1 若双方均采取最优策略，先手是否必胜？
2 若双方随机取数，先手获胜的概率为多少？
$n,ai \le 100$

---

状态比较难想，核心是找到一个划分阶段的顺序：根据$GCD$划分阶段

$GCD$是只会减小不会增加的

课件上的状态是$f[i][j]$表示当前$GCD$为$i$，没选的$i$的倍数有$j$个，感觉有点奇怪...

看了一下官方题解，意识到只要记录$j$为当前已经选的有$j$个就好了，已经选的一定是$i$的倍数，这样就和其他的状态比较像了

转移还是比较好想的

$1.\ f[i][j] \rightarrow f[i][j+1]\ :\ j<mult[i]$

$2.\ f[i][j] \rightarrow f[gcd(i,k)][j+1]\ :\ 1 \le gcd(i,k) \le i$

记忆化搜索倒推就行了

PS：给$gcd$加上记忆化之后$0s$就跑过去了....

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=;
const double eps=1e-;
inline int read(){
    char c=getchar();int x=,f=;
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x*f;
}
int n,a[N];
int g[N][N];
int gcd(int a,int b){return g[a][b] ? g[a][b] : g[a][b]=(b==?a:gcd(b,a%b));}
//int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
int f[N][N];
double p[N][N];
bool dfsWin(int u,int c){//printf("dfsWin %d %d\n",u,c);
    int &re=f[u][c];
    if(c==n) re=;//has chosen all
    if(u==) re=;//win
    if(re!=-) return re;
 
    re=;
    int mult=;
    for(int i=;i<=n;i++) if(gcd(u,a[i])==u) mult++;
    if(c<mult&&!dfsWin(u,c+)) re=;
    else{
        for(int i=;i<=n;i++)
            if(gcd(u,a[i])>&&gcd(u,a[i])!=u)
                if(!dfsWin(gcd(u,a[i]),c+)) {re=;break;}
    }
    return re;
}
double dfsPro(int u,int c){//printf("dfsPro %d %d\n",u,c);
    double &re=p[u][c];
    if(c==n) re=0.0;
    if(u==) re=1.0;
    if(re>-0.9) return re;
 
    re=0.0;
    int mult=;
    for(int i=;i<=n;i++) if(gcd(u,a[i])==u) mult++;
    if(c<mult) re+= (double)(mult-c) / (n-c) * (-dfsPro(u,c+));
    for(int i=;i<=n;i++)
        if(gcd(u,a[i])>&&gcd(u,a[i])!=u)
            re+=(double) / (n-c) *(-dfsPro(gcd(u,a[i]),c+));
    if(abs(re)<eps) re=;
    return re;
}
int main(){
    freopen("in","r",stdin);
    int T=read();
    while(T--){
        n=read(); int g=;
        for(int i=;i<=n;i++) a[i]=read(),g=gcd(a[i],g);
        if(g>){printf("%d %.4lf\n",n&,double(n&));continue;}
 
        for(int i=;i<N;i++) for(int j=;j<N;j++) f[i][j]=-,p[i][j]=-1.0;
        int flag=dfsWin(,);
        printf("%d ",flag);
 
        double prob=dfsPro(,);
        printf("%.4lf\n",prob);
    }
}