2021.1.4 省选训练 CF113D Museum

暴力卡常

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79

//可以加火车头（亲测用处不大） #include<iostream> #include<cstdio> #include<algorithm> #include<cstring> #define N 23//开大会TLE（？） #define M 550 using namespace std; inline int read() {     int x=0,f=1;char ch=getchar();     while(ch<'0' || ch>'9'){if(ch=='-') f=-1;ch=getchar();}     while(ch>='0' && ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}     return x*f; } int n,m; double a[N],b[N],out[N],f[2][N][N],ans[N]; double E=1e-15; int head[N],nxt[M],v[M],cnt; inline void add(int x,int y) {     v[++cnt]=y;     nxt[cnt]=head[x];     head[x]=cnt; } inline void DP(int s1,int s2,int x,int y) {     double tmp[2][2];//减少乘法运算     tmp[0][0]=b[s1]*b[s2]*f[x][s1][s2];//期望的运算方法（错）     tmp[0][1]=b[s1]*a[s2]*f[x][s1][s2];     tmp[1][0]=a[s1]*b[s2]*f[x][s1][s2];     tmp[1][1]=a[s1]*a[s2]*f[x][s1][s2];     for(int i=head[s1];i;i=nxt[i])         for(int j=head[s2];j;j=nxt[j])             f[y][v[i]][v[j]]+=tmp[0][0];     for(int i=head[s1];i;i=nxt[i]) f[y][v[i]][s2]+=tmp[0][1];     for(int j=head[s2];j;j=nxt[j]) f[y][s1][v[j]]+=tmp[1][0];     f[y][s1][s2]+=tmp[1][1]; } int main() {     int s1,s2;     n=read();m=read();s1=read();s2=read();     if(n==1){printf("1\n");return 0;}//注意特判（错）     if(s1==s2)     {         for(int i=1;i<=n;i++)             if(i==s1) printf("1 ");             else printf("0 ");         return 0;     }     int x,y;     for(int i=1;i<=m;i++)     {         x=read();y=read();         add(x,y),add(y,x);         out[x]+=1.0;         out[y]+=1.0;     }     for(int i=1;i<=n;i++) scanf("%lf",&a[i]),b[i]=(1.0-a[i])/out[i];     x=0;y=1;     f[x][s1][s2]=1.0;     int t=150000;     while(t--)//x->y  推表式     {         for(int i=1;i<=n;i++)             for(int j=1;j<=n;j++)             {                 if(f[x][i][j]<=E) continue;//忽略过小的精度误差（卡精度）                 if(i==j) continue;                 DP(i,j,x,y);             }         for(int k=1;k<=n;k++) ans[k]+=f[y][k][k];         x^=1;y^=1;         memset(f[y],0,sizeof(f[y]));//注意清空（错）     }     for(int i=1;i<=n;i++) printf("%.12lf ",ans[i]);//需要手动调整输出精度     return 0; }

　　

高斯消元（待填）