POJ 3071 Football 概率DP

#include<cstdio>
#include<cstring>
const int LMT=7;
double dp[LMT+3][1<<LMT],mat[1<<LMT][1<<LMT];
int n;
void init(void)
{
	int i;
	memset(dp,0,sizeof(dp));
	for(i=0;i<(1<<n);i++)dp
[i]=1.0;
}
int main(void)
{
	int i,j,t,lim,pos,d;
	double ans;
	while(~scanf("%d",&n)&&~n)
	{
		init();
		ans=-1;pos=-1;
		lim=1<<n;
		for(i=0;i<lim;i++)
			for(j=0;j<lim;j++)
				scanf("%lf",&mat[i][j]);
		for(i=n-1;i>=0;i--)
			for(j=0;j<lim;j++)
			{
				d=j/(1<<(n-i-1));
				if(!(d&1))
				{
					for(t=(d+1)*(1<<(n-i-1));t<(d+2)*(1<<(n-i-1));t++)
						dp[i][j]+=dp[i+1][j]*dp[i+1][t]*mat[j][t];
				}
				else
				{
					for(t=(d-1)*(1<<(n-i-1));t<d*(1<<(n-i-1));t++)
					dp[i][j]+=dp[i+1][j]*dp[i+1][t]*mat[j][t];
				}
			}
			for(i=0;i<lim;i++)
				if(ans<dp[0][i])
				{
					ans=dp[0][i];
					pos=i+1;
				}
                printf("%d\n",pos);
	}
	return 0;
}