結果

問題 No.248 ミラー君の宿題
コンテスト
ユーザー 👑 kmjp
提出日時 2015-07-21 23:13:47
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 1,710 ms / 5,000 ms
コード長 2,106 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,307 ms
コンパイル使用メモリ 162,188 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-23 14:15:20
合計ジャッジ時間 17,597 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;

#undef _P
#define _P(...) (void)printf(__VA_ARGS__)
#define FOR(x,to) for(x=0;x<to;x++)
#define FORR(x,arr) for(auto& x:arr)
#define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
//-------------------------------------------------------

int T;
int N;
ll P[20];
int bi[20];
int minbi[1<<15];
double bip[20][65];
double bipsum[20][65];
double memo[1<<15];

double dodo(int pm) {
	
	if(__builtin_popcount(pm)<=1) return __builtin_popcount(pm);
	if(memo[pm]>=0) return memo[pm];
	
	vector<int> id;
	double ret=0;
	double left=1, p2b=1;
	int mask,i,x;
	
	FOR(i,N) if(pm&(1<<i)) id.push_back(i), p2b *= (P[i]-1.0)/P[i];
	int n=id.size();
	
	
	FOR(mask,1<<n) {
		double p1=1;
		int nm[2]={};
		FOR(i,n) {
			if(mask&(1<<i)) p1 *= 1.0/P[id[i]], nm[0] |= 1<<id[i];
			else p1 *= (P[id[i]]-1.0)/P[id[i]], nm[1] |= 1<<id[i];
		}
		if(nm[0] && nm[1]) {
			int mb=minbi[nm[0]];
			
			double p2=0;
			for(x=1;x<=63;x++) {
				double hoge=1;
				FOR(i,n) {
					if(mask&(1<<i)) hoge *= bip[id[i]][x];
					else hoge *= bipsum[id[i]][x-1];
				}
				p2+=hoge;
			}
			
			double pat = p1 + p2*p2b;
			left -= pat;
			ret += pat * (dodo(nm[0])+dodo(nm[1]));
		}
	}
	
	return memo[pm]=(1+ret)/(1-left);
}

void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>N;
	FOR(i,1<<N) memo[i]=-1;
	
	ZERO(bip);
	ZERO(bipsum);
	FOR(i,N) {
		cin>>P[i];
		ll p=P[i]-1;
		bi[i]=0;
		while(p%2==0) bi[i]++, p/=2;
		
		bipsum[i][0]=bip[i][0]=pow(0.5,bi[i]);
		for(x=1;x<=64;x++) {
			if(x<=bi[i]) bip[i][x]=pow(0.5,bi[i]+1-x);
			bipsum[i][x]=bipsum[i][x-1]+bip[i][x];
		}
		
	}
	
	FOR(i,1<<N) {
		minbi[i]=100;
		FOR(j,N) if(i & (1<<j)) minbi[i]=min(minbi[i],bi[j]);
	}
	
	_P("%.12lf\n",dodo((1<<N)-1));
	
	
}


int main(int argc,char** argv){
	string s;int i;
	if(argc==1) ios::sync_with_stdio(false);
	FOR(i,argc-1) s+=argv[i+1],s+='\n';
	FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);
	cin>>T;
	FOR(i,T) solve();
	return 0;
}
0