HDU.5819.Knights(概率DP)

题目链接

参考一下这的.

\(Description\)
数轴上有n个骑士，分别位于1,2,3,...,n，它们的移动速度相同，初始移动方向已知。当两个骑士相遇时，各有50%的概率获胜，失败的骑士就死了。
当骑士移动到0,n+1位置时方向会反转。求最右边的骑士最后存活的概率。

\(Solution\)
题目中速度、时间都是无关紧要的。第n个骑士若要赢，一定是打败了左边所有存活的骑士。
那么设\(f[i][j]\)表示前i个骑士中有j个骑士存活的概率。
考虑怎么求。第i个骑士如果向左，那它应把前i-1个骑士中k-j个向右的骑士都打败(然后go_die->1/2)，才能剩下j个.
即\(f[i][j]=\sum_{k=j}^{i-1}f[i-1][k]*(\frac{1}{2})^{k-j+1}\)
把\(k=j\)时的分离出来，可以化简为\[f[i][j]=\frac{1}{2}*(f[i-1][j]+f[i][j+1])\]
如果骑士i向右，那么\[f[i][j]=f[i-1][j-1]\]
但是如果只有一个骑士向右，还可能是它打败了左边所有骑士，到0点折回的。
所以给每个\(f[i][1]+=\sum_{j=1}^{i-1}f[i-1][j]*(\frac{1}{2})^j\)
即\[(f[i][1]+=\frac{1}{2}*(f[i-1][1]+f[i][2]))=f[i-1][1]+f[i][2]\]
最终答案\(Ans=\sum_{i=1}^{n-1}f[n-1][i]*(\frac{1}{2})^i\)，即\[Ans=\frac{f[n][1]}{2}\]

注意骑士n是要强制向左走的（本来就是）。

#include <cstdio>
#include <cctype>
#include <cstring>
#define gc() getchar()
#define inv (500000004)
#define mod (1000000007)//除以2全都化为乘逆元(题目要求)，这样也能直接用整数 
const int N=1005;

int n,f[N][N];

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}

int main()
{
    int T=read();
    for(int cas=1; cas<=T; ++cas)
    {
        memset(f,0,sizeof f);
        n=read(), read();
        f[1][1]=1;
        for(int i=2; i<=n; ++i)
        {
            if(read()&&i!=n)
                for(int j=1; j<=i; ++j)
                    f[i][j]=f[i-1][j-1];
            else{
                for(int j=i-1; j>1; --j)
                    f[i][j]=1ll*inv*(f[i-1][j]+f[i][j+1])%mod;
                f[i][1]=(f[i-1][1]+f[i][2])%mod;
            }
        }
        printf("Case #%d: %I64d\n",cas,1ll*f[n][1]*inv%mod);
    }
    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/8664745.html