結果

問題 No.1688 Veterinarian
コンテスト
ユーザー kotatsugame
提出日時 2021-09-24 21:33:00
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 40 ms / 3,000 ms
コード長 984 bytes
コンパイル時間 518 ms
コンパイル使用メモリ 72,448 KB
実行使用メモリ 5,632 KB
最終ジャッジ日時 2024-07-05 10:08:12
合計ジャッジ時間 1,225 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:6:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
    6 | main()
      | ^~~~

ソースコード

diff # 

#include<iostream>
#include<iomanip>
using namespace std;
int A,B,C,N;
double dp[2][51][51][51];
main()
{
	cin>>A>>B>>C>>N;
	int now=0;
	dp[now][A][B][C]=1;
	for(int _=0;_<N;_++)
	{
		int nxt=1-now;
		for(int i=1;i<=A;i++)for(int j=1;j<=B;j++)for(int k=1;k<=C;k++)dp[nxt][i][j][k]=0;
		for(int i=1;i<=A;i++)for(int j=1;j<=B;j++)for(int k=1;k<=C;k++)
		{
			int n=i+j+k;
			int all=n*(n-1)/2;
			int usd=all;
			{
				int a=i*(i-1)/2;
				dp[nxt][i-1][j][k]+=dp[now][i][j][k]*a/all;
				usd-=a;
			}
			{
				int b=j*(j-1)/2;
				dp[nxt][i][j-1][k]+=dp[now][i][j][k]*b/all;
				usd-=b;
			}
			{
				int c=k*(k-1)/2;
				dp[nxt][i][j][k-1]+=dp[now][i][j][k]*c/all;
				usd-=c;
			}
			dp[nxt][i][j][k]+=dp[now][i][j][k]*usd/all;
		}
		now=nxt;
	}
	double a=0,b=0,c=0;
	for(int i=1;i<=A;i++)for(int j=1;j<=B;j++)for(int k=1;k<=C;k++)
	{
		a+=dp[now][i][j][k]*(A-i);
		b+=dp[now][i][j][k]*(B-j);
		c+=dp[now][i][j][k]*(C-k);
	}
	cout<<fixed<<setprecision(16)<<a<<" "<<b<<" "<<c<<endl;
}
0