結果

問題 No.162 8020運動
ユーザー kuuso1kuuso1
提出日時 2015-03-08 17:58:57
言語 C90
(gcc 11.4.0)
結果
RE  
実行時間 -
コード長 1,814 bytes
コンパイル時間 381 ms
コンパイル使用メモリ 27,816 KB
実行使用メモリ 39,852 KB
最終ジャッジ日時 2023-09-06 21:55:29
合計ジャッジ時間 40,407 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <stdlib.h>

double dp[2][1<<14];

int BitCnt(int s){
	int ret=0;
	int i;
	for(i=0;i<20;i++){
		if( ((s>>i)&1)>0 )ret++;
	}
	return ret;
}

int main(){
	double P0,P1,P2;
	int A;
	scanf("%d\n",&A);
	scanf("%f %f %f\n",&P0,&P1,&P2);
	P0/=100.0;P1/=100.0;P2/=100.0;
	int s,t,i;
	
	double **table;
	table=malloc(sizeof(double *) *(1<<14));
	int Size[1<<14];
	for(s=0;s<(1<<14);s++){
		Size[s]=1<<BitCnt(s);
		table[s]=malloc(sizeof(double)*Size[s]);
	}
	
	
	for(s=0;s<1<<14;s++){
		for(t=0;t<=s;t++){
			if((s|t)!=s){
				t+=(t&-t)-1;
				continue;
			}
			double p=1.0;
			for(i=0;i<14;i++){
				if(((s>>i)&1)==0)continue;
				if(((t>>i)&1)==0){
					if(i==0){
						if(((s>>(i+1))&1)==0){
							p*=P0;
						}else{
							p*=P1;
						}
					}else{
						if(((s>>(i+1))&1)==0 && ((s>>(i-1))&1)==0){
							p*=P0;
						}else if(((s>>(i+1))&1)==1 && ((s>>(i-1))&1)==1){
							p*=P2;
						}else{
							p*=P1;
						}
					}
				}else{
					if(i==0){
						if(((s>>(i+1))&1)==0){
							p*=(1.0-P0);
						}else{
							p*=(1.0-P1);
						}
					}else{
						if(((s>>(i+1))&1)==0 && ((s>>(i-1))&1)==0){
							p*=(1.0-P0);
						}else if(((s>>(i+1))&1)==1 && ((s>>(i-1))&1)==1){
							p*=(1.0-P2);
						}else{
							p*=(1.0-P1);
						}
					}
				}
			}
			table[s][t]=p;
		}
	}
	
		
	int now=0,next=1;
	for(s=0;s<1<<14;s++)dp[now][s]=0.0;
	dp[0][(1<<14)-1]=1.0;

	for(i=0;i<80-A;i++){
		now=i%2;
		next=now^1;
		for(s=0;s<1<<14;s++)dp[next][s]=0.0;
		
		for(s=0;s<(1<<14);s++){
			for(t=0;t<=s;t++){
				if((s|t)!=s){
					t+=(t&-t)-1;
					continue;
				}
				dp[next][t]=table[s][t]*dp[now][s];
			}
		}
	}
	
	
	double e=0;
	for(s=0;s<(1<<14);s++){
		e+=BitCnt(s)*dp[next][s];
	}
	printf("%.6f\n",2.0*e);
	
	for(s=0;s<(1<<14);s++){
		free(table[s]);
	}
	free(table);
	
	return 0;
}
0