結果

問題 No.162 8020運動
ユーザー kuuso1kuuso1
提出日時 2015-03-21 01:46:49
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 1,269 ms / 5,000 ms
コード長 2,878 bytes
コンパイル時間 2,986 ms
コンパイル使用メモリ 113,252 KB
実行使用メモリ 95,488 KB
最終ジャッジ日時 2024-06-28 23:52:51
合計ジャッジ時間 33,862 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 663 ms
93,056 KB
testcase_01 AC 937 ms
94,080 KB
testcase_02 AC 1,095 ms
94,592 KB
testcase_03 AC 1,244 ms
95,360 KB
testcase_04 AC 1,240 ms
95,360 KB
testcase_05 AC 1,269 ms
95,232 KB
testcase_06 AC 1,177 ms
95,488 KB
testcase_07 AC 1,229 ms
95,232 KB
testcase_08 AC 1,178 ms
95,360 KB
testcase_09 AC 679 ms
92,672 KB
testcase_10 AC 1,061 ms
94,720 KB
testcase_11 AC 915 ms
94,080 KB
testcase_12 AC 779 ms
93,568 KB
testcase_13 AC 1,207 ms
95,360 KB
testcase_14 AC 702 ms
92,800 KB
testcase_15 AC 712 ms
93,184 KB
testcase_16 AC 833 ms
93,696 KB
testcase_17 AC 739 ms
92,928 KB
testcase_18 AC 1,159 ms
95,232 KB
testcase_19 AC 1,071 ms
94,592 KB
testcase_20 AC 1,102 ms
94,848 KB
testcase_21 AC 1,135 ms
94,976 KB
testcase_22 AC 1,052 ms
94,464 KB
testcase_23 AC 1,131 ms
94,848 KB
testcase_24 AC 1,112 ms
94,848 KB
testcase_25 AC 1,173 ms
94,976 KB
testcase_26 AC 681 ms
92,672 KB
testcase_27 AC 940 ms
94,080 KB
testcase_28 AC 1,160 ms
95,360 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;
using System.Collections;
using System.Collections.Generic;
using System.Threading.Tasks;
class TEST{
	static void Main(){
		Sol mySol =new Sol();
		mySol.Solve();
	}
}

class Sol{
	public void Solve(){
		
		//状態 s から推移できる状態数は 2^(sのビット数)
		//推移できる状態もstate[s]に列挙しておく
		double[][] table=new double[1<<14][];
		int[][] state=new int[1<<14][];
		for(int s=0;s<1<<14;s++){
			table[s]=new double[1<<BitCnt(s)];
			state[s]=new int[1<<BitCnt(s)];
		}
		
		//状態 s -> 状態 t になる確率を計算する。
		Parallel.For(0,1<<14,s=>{
			int ptr=0;
			for(int t=0;t<=s;t++){
				if((s|t)!=s){
					t+=(t&-t)-1;	//余計なビットが立った時にループを一部省略する。
					continue;
				}
				double p=1.0;
				for(int 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][ptr]=p;
				state[s][ptr]=t;
				ptr++;
			}
		});
		
		//あとは普通にdpする。
		double[][] dp=new double[2][];
		dp[0]=new double[1<<14];
		
		dp[0][(1<<14)-1]=1.0;
		int now=0,next=1;
		for(int i=0;i<80-A;i++){
			now=i%2;
			next=now^1;
			dp[next]=new double[1<<14];
			for(int s=0;s<(1<<14);s++){
				for(int t=0;t<table[s].Length;t++){
					dp[next][state[s][t]]+=table[s][t]*dp[now][s];
				}
			}
		}
		
		double e=0;
		for(int s=0;s<(1<<14);s++){
			e+=BitCnt(s)*dp[next][s];
		}
		Console.WriteLine(2.0*e);
		
	}
	
	int BitCnt(int s){
		int ret=0;
		for(int i=0;i<20;i++){
			if( ((s>>i)&1)>0 )ret++;
		}
		return ret;
	}
	
	int A;
	double P0,P1,P2;
	public Sol(){
		A=ri();
		var d=rda();
		P0=d[0]/100.0;P1=d[1]/100.0;P2=d[2]/100.0;
	}

	static String rs(){return Console.ReadLine();}
	static int ri(){return int.Parse(Console.ReadLine());}
	static long rl(){return long.Parse(Console.ReadLine());}
	static double rd(){return double.Parse(Console.ReadLine());}
	static String[] rsa(){return Console.ReadLine().Split(' ');}
	static int[] ria(){return Array.ConvertAll(Console.ReadLine().Split(' '),e=>int.Parse(e));}
	static long[] rla(){return Array.ConvertAll(Console.ReadLine().Split(' '),e=>long.Parse(e));}
	static double[] rda(){return Array.ConvertAll(Console.ReadLine().Split(' '),e=>double.Parse(e));}
}
0