• 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
    ",cas,1ll*f[n][1]*inv%mod);
    	}
    	return 0;
    }
    
  • 相关阅读:
    Django之数据库表的创建和ORM相关操作
    Django后续和Ajax初识
    阿里云Maven中央仓库配置
    java/javascript 时间操作工具类
    原生javascript实现文件异步上传
    MySQL中的存储函数和存储过程的简单示例
    java同步锁的正确使用
    浅谈javascript的面向对象思想
    java与javascript对cookie操作的工具类
    json字符串与json对象的相互转换
  • 原文地址:https://www.cnblogs.com/SovietPower/p/8664745.html
Copyright © 2020-2023  润新知