結果

問題 No.1151 チャレンジゲーム
ユーザー tails
提出日時 2020-11-11 11:37:04
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 10 ms / 2,000 ms
コード長 1,157 bytes
コンパイル時間 1,438 ms
コンパイル使用メモリ 34,048 KB
実行使用メモリ 17,996 KB
最終ジャッジ日時 2024-07-22 18:23:12
合計ジャッジ時間 2,348 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 50
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:9:1: warning: return type defaults to 'int' [-Wimplicit-int]
    9 | main(){
      | ^~~~
main.c: In function 'main':
main.c:11:9: warning: implicit declaration of function 'read' [-Wimplicit-function-declaration]
   11 |         read(0,rbuf,sizeof rbuf);
      |         ^~~~
main.c:61:9: warning: implicit declaration of function 'printf' [-Wimplicit-function-declaration]
   61 |         printf("%f",d[0][0][0]);
      |         ^~~~~~
main.c:1:1: note: include '<stdio.h>' or provide a declaration of 'printf'
  +++ |+#include <stdio.h>
    1 | #pragma GCC optimize("Ofast")
main.c:61:9: warning: incompatible implicit declaration of built-in function 'printf' [-Wbuiltin-declaration-mismatch]
   61 |         printf("%f",d[0][0][0]);
      |         ^~~~~~
main.c:61:9: note: include '<stdio.h>' or provide a declaration of 'printf'

ソースコード

diff #

#pragma GCC optimize("Ofast")
#pragma GCC target("avx2")

#define RD(v) int v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}

int a[1<<10];
double d[1<<10][1<<10][2];

main(){
	char rbuf[512];
	read(0,rbuf,sizeof rbuf);
	char*rp=rbuf;
	RD(n);
	for(int i=0;i<n;++i){
		RD(t);
		a[1<<i]=t;
	}

	int mm=(1<<n)-1;
	{
		for(int mi=mm;mi>=0;--mi){
			int s=0;
			for(int i=0;i<n;++i){
				int ai=a[1<<i];
				s+=mi&1<<i?ai:-ai;
			}
			int mj=mm^mi;
			if(s>0){
				d[mi][mj][0]=1;
				d[mi][mj][1]=1;
			}else{
				d[mi][mj][0]=0;
				d[mi][mj][1]=0;
			}
		}
	}
	for(int m=mm-1;m>=0;--m){
		int mi=0;
		do{
			int mj=m^mi;
			double b0ij=0;
			double b1ij=1;
			for(int xi=mm&~m;xi;xi&=xi-1){
				int i=xi&-xi;
				double b0j=1;
				double b1i=0;
				for(int xj=mm&~m;xj;xj&=xj-1){
					int j=xj&-xj;
					double b0=(a[j]*d[mi|i][mj][1]+(a[i]-1)*d[mi][mj|j][0])/(a[i]+a[j]-1);
					double b1=((a[i]-1)*d[mi|j][mj][1]+a[j]*d[mi][mj|i][0])/(a[i]+a[j]-1);
					if(b0j>b0) b0j=b0;
					if(b1i<b1) b1i=b1;
				}
				if(b0ij<b0j) b0ij=b0j;
				if(b1ij>b1i) b1ij=b1i;
			}
			d[mi][mj][0]=b0ij;
			d[mi][mj][1]=b1ij;
		}while(mi=mi-1&m);
	}
	printf("%f",d[0][0][0]);
}
0