結果

問題 No.162 8020運動
ユーザー kuuso1kuuso1
提出日時 2015-03-21 01:46:49
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 1,430 ms / 5,000 ms
コード長 2,878 bytes
コンパイル時間 3,446 ms
コンパイル使用メモリ 108,784 KB
実行使用メモリ 103,712 KB
最終ジャッジ日時 2023-09-11 09:12:01
合計ジャッジ時間 40,167 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 804 ms
101,200 KB
testcase_01 AC 1,093 ms
98,300 KB
testcase_02 AC 1,263 ms
97,392 KB
testcase_03 AC 1,422 ms
102,096 KB
testcase_04 AC 1,429 ms
103,712 KB
testcase_05 AC 1,427 ms
103,588 KB
testcase_06 AC 1,404 ms
99,816 KB
testcase_07 AC 1,430 ms
100,092 KB
testcase_08 AC 1,404 ms
101,984 KB
testcase_09 AC 804 ms
96,980 KB
testcase_10 AC 1,247 ms
101,288 KB
testcase_11 AC 1,084 ms
102,688 KB
testcase_12 AC 926 ms
99,732 KB
testcase_13 AC 1,403 ms
102,152 KB
testcase_14 AC 829 ms
97,592 KB
testcase_15 AC 824 ms
97,528 KB
testcase_16 AC 971 ms
100,056 KB
testcase_17 AC 870 ms
99,404 KB
testcase_18 AC 1,369 ms
101,808 KB
testcase_19 AC 1,240 ms
102,972 KB
testcase_20 AC 1,282 ms
101,328 KB
testcase_21 AC 1,304 ms
101,520 KB
testcase_22 AC 1,246 ms
97,384 KB
testcase_23 AC 1,299 ms
103,412 KB
testcase_24 AC 1,263 ms
101,232 KB
testcase_25 AC 1,339 ms
101,732 KB
testcase_26 AC 805 ms
99,104 KB
testcase_27 AC 1,118 ms
98,680 KB
testcase_28 AC 1,376 ms
99,516 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