Codeforces Beta Round #16 (Div. 2 Only), problem: (E) Fish 概率DP

Even

于 2013-04-14 21:06:08 发布

阅读量1k

点赞数

分类专栏： ACM

本文链接：https://blog.youkuaiyun.com/cqlf__/article/details/8801222

版权

ACM 专栏收录该内容

13 篇文章

订阅专栏

本文介绍了一个使用C++实现的动态规划算法案例，该案例通过预处理位运算来计算特定状态转移的概率问题。文章详细展示了如何初始化数据结构、进行状态转移计算以及输出最终结果的过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<cstdio>
#include<cstring>
const int LMT=(1<<18)+10;
double dp[LMT],a[19][19],__one[LMT];
void get_one(void)
{
	int i,x;
	for(i=1;i<LMT;i++)
	{
		x=i;
		while(x)
		{
			if(x&1)__one[i]++;
			x>>=1;
		}
	}
}
int main(void)
{
	int n,lim,i,ii,jj;
	get_one();
	while(~scanf("%d",&n))
	{
		memset(dp,0,sizeof(dp));
	     lim=(1<<n)-1;
	     dp[lim]=1.0;
		 for(ii=0;ii<n;ii++)
			 for(jj=0;jj<n;jj++)
				 scanf("%lf",&a[ii][jj]);
		 for(i=lim-1;i>0;i--)
			 if(__one[i]>0)
			 {
				 for(jj=0;jj<n;++jj)
					 if(!(i&(1<<jj)))
						 for(ii=0;ii<n;++ii)
							 if(i&(1<<ii))
								 dp[i]+=a[ii][jj]*dp[i|(1<<jj)];
								 dp[i]/=(__one[i]+1)*__one[i]*0.5;
			 }
			 for(i=0;i<n;i++)
				 printf("%.6lf ",dp[1<<i]);
			 printf("\n");
	}
	return 0;
}