結果

問題 No.1688 Veterinarian
コンテスト
ユーザー tails
提出日時 2021-09-25 16:24:34
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 22 ms / 3,000 ms
コード長 635 bytes
コンパイル時間 330 ms
コンパイル使用メモリ 31,104 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-05 12:06:51
合計ジャッジ時間 1,048 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます
コンパイルメッセージ
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 'scanf' [-Wimplicit-function-declaration]
   11 |         scanf("%ld%ld%ld%ld",&a,&b,&c,&n);
      |         ^~~~~
main.c:1:1: note: include '<stdio.h>' or provide a declaration of 'scanf'
  +++ |+#include <stdio.h>
    1 | #pragma GCC optimize("Ofast")
main.c:11:9: warning: incompatible implicit declaration of built-in function 'scanf' [-Wbuiltin-declaration-mismatch]
   11 |         scanf("%ld%ld%ld%ld",&a,&b,&c,&n);
      |         ^~~~~
main.c:11:9: note: include '<stdio.h>' or provide a declaration of 'scanf'
main.c:31:9: warning: implicit declaration of function 'printf' [-Wimplicit-function-declaration]
   31 |         printf("%.8f %.8f %.8f",d[a][b][c],d[b][c][a],d[c][a][b]);
      |         ^~~~~~
main.c:31:9: note: include '<stdio.h>' or provide a declaration of 'printf'
main.c:31:9: warning: incompatible implicit declaration of built-in function 'printf' [-Wbuiltin-declaration-mismatch]
main.c:31:9: note: include '<stdio.h>' or provide a declaration of 'printf'

ソースコード

diff # 

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

#define rep(v,e) for(long v=0;v<e;++v)
#define rrep3(v,s,e) for(long v=e;--v>=s;)

double d[51][51][51];

main(){
	long a,b,c,n;
	scanf("%ld%ld%ld%ld",&a,&b,&c,&n);
	rep(i,n){
		rrep3(x,2,51){
			double xx=x*(x-1);
			rrep3(y,1,51){
				double yy=y*(y-1);
				rrep3(z,1,51){
					double zz=z*(z-1);
					long s=x+y+z;
					double ss=s*(s-1);
					double a=d[x][y][z];
					d[x][y][z]=a+(double)(
						xx*(d[x-1][y][z]-a+1)+
						yy*(d[x][y-1][z]-a)+
						zz*(d[x][y][z-1]-a)
					)/(double)ss;
				}
			}
		}
	}
	printf("%.8f %.8f %.8f",d[a][b][c],d[b][c][a],d[c][a][b]);
}
0