結果

問題 No.162 8020運動
コンテスト
ユーザー kuuso1
提出日時 2015-03-08 17:23:33
言語 C#(csc)
(csc 3.9.0)
結果
MLE  
実行時間 -
コード長 2,483 bytes
コンパイル時間 3,589 ms
コンパイル使用メモリ 108,288 KB
実行使用メモリ 824,832 KB
最終ジャッジ日時 2024-06-24 16:04:28
合計ジャッジ時間 4,602 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample MLE * 1 -- * 2
other -- * 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(){
		
		double[][] table=new double[1<<14][];
		for(int i=0;i<1<<14;i++){
			table[i]=new double[1<<14];
		}
		
		for(int s=0;s<1<<14;s++){
			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][t]=p;
			}
		}
		
		
		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<=s;t++){
					if((s|t)!=s){
						t+=(t&-t)-1;
						continue;
					}
					dp[next][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