hdu 6710 2019 CCPC 网络赛 09题

题意：就是给你一个二分图，左右两边两点之间存在边的概率是0.5，问两点的期望距离是多少？

思路：组合数学和dp。

 

设dp[i][j][k][l]

第一个状态是当前的步长

第二个状态是左边还有多少个点没链接

第三个状态是右边还有多少个点没链接

第四个状态是当前可用来链接下一批的点是多少个

 

举一个转移的例子，当我们从第一个点出发可以链接到右边的m个点（不包含终点），则第二个状态相应减少m个点，第四状态则为m。最后只要统计偶数步时的概率乘上相应步数，就可以得到答案。

 

一批点链接到达一批点，则使用组合数学进行处理。

PS：中午在D教学楼困死了，比赛时题目都没看清。
 

#include<stdio.h>
#include<algorithm>
using namespace std;

#define ll long long

int n,m,mod;
int dp[62][32][32][32];
int c[32][32];
int p[32][32];
int t[32];
int q[32][32];
int qp[32][32][32][32];


ll qpow(ll x,ll y){
    ll e=1;
    while(y){
        if(y&1)e=e*x%mod;
        x=x*x%mod;
        y>>=1;
    }
    return e;
}

void init(){

    int lim=max(n,m)+1;
    for(int i=0;i<=lim;i++){
        c[i][0]=1;
        for(int j=1;j<=i;j++){
            c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
        }
    }

    t[0]=1;
    t[1]=qpow(2,mod-2);
    for(int i=2;i<=lim;i++)
        t[i]=1LL*t[i-1]*t[1]%mod;

    for(int i=1;i<=lim;i++){
        q[i][0]=1;
        q[i][1]=t[i];
        for(int j=2;j<=lim;j++){
            q[i][j]=(1LL*q[i][j-1]*t[i])%mod;
        }
    }

    for(int i=1;i<=lim;i++){
        int e=0;
        for(int j=1;j<=i;j++){
            e=(e+1LL*c[i][j]*t[i])%mod;
        }
        p[i][0]=1;
        for(int j=1;j<=lim;j++){
            p[i][j]=(1LL*p[i][j-1]*e)%mod;
        }
    }

    for(int i=0;i<=lim;i++){
        for(int j=0;j<=lim;j++){
            for(int k=0;k<=lim;k++){
                for(int l=0;l<=lim;l++){
                    qp[i][j][k][l]=(1LL*p[i][j]*q[k][l])%mod;
                }
            }
        }
    }
}

int main(){

    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d %d %d",&n,&m,&mod);

        init();

        n--;
        m--;

        for(int i=0;i<=n+m+1;i++)
            for(int j=0;j<=n+1;j++)
                for(int k=0;k<=m+1;k++)
                    for(int l=0;l<=30;l++)
                        dp[i][j][k][l]=0;

        dp[0][n][m][1]=1;
        for(int i=1;i<=(n+m+1);i++){
            for(int j=0;j<=n;j++){
                for(int k=0;k<=m;k++){
                    if(i&1){
                        for(int l=1;l<=k;l++){
                            int u=k-l;
                            int v=max(n-j,1);
                            for(int o=1;o<=v;o++){
                                dp[i][j][u][l]=(dp[i][j][u][l]+1LL*dp[i-1][j][k][o]*c[k][l]%mod*qp[o][l][o][1+u])%mod;
                            }
                        }
                    }
                    else{
                        for(int l=1;l<=j;l++){
                            int u=j-l;
                            int v=m-k;
                            for(int o=1;o<=v;o++){
                                dp[i][u][k][l]=(dp[i][u][k][l]+1LL*dp[i-1][j][k][o]*c[j][l]%mod*qp[o][l][o][u])%mod;
                            }
                        }

                    }

                }
            }
        }

        ll ans=0;
        for(int i=0;i<=n+m+1;i+=2){
            for(int j=0;j<=n;j++){
                for(int k=0;k<=m;k++){
                    for(int l=1;l<=n+1;l++){
                        ans=(ans+1LL*(i+1)*dp[i][j][k][l]%mod*p[l][1])%mod;
                    }
                }
            }
        }

        printf("%lld\n",(ans+mod)%mod);
    }
    return 0;
}