結果

問題 No.162 8020運動
コンテスト
ユーザー kuuso1
提出日時 2015-03-08 18:42:19
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 1,495 ms / 5,000 ms
コード長 2,848 bytes
コンパイル時間 916 ms
コンパイル使用メモリ 114,788 KB
実行使用メモリ 102,412 KB
最終ジャッジ日時 2024-06-24 16:16:28
合計ジャッジ時間 38,996 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
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;
 
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 になる確率を計算する。
		for(int s=0;s<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