結果

問題 No.301 サイコロで確率問題 (1)
コンテスト
ユーザー piyoko_212
提出日時 2015-12-28 16:54:10
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 18 ms / 1,000 ms
コード長 1,230 bytes
コンパイル時間 389 ms
コンパイル使用メモリ 47,524 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-19 07:36:17
合計ジャッジ時間 969 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:53:20: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   53 |         int T;scanf("%d",&T);
      |               ~~~~~^~~~~~~~~
main.cpp:55:34: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   55 |                 long long a;scanf("%lld",&a);
      |                             ~~~~~^~~~~~~~~~~

ソースコード

diff #
raw source code 

#include<stdio.h>
#include<algorithm>
#include<vector>
using namespace std;
const double EPS=1e-8;
typedef vector<double>vec;
typedef vector<vec>mat;
inline double ABS(double a){return max(a,-a);}
vec gauss_jordan(const mat &A,const vec &b){
	int n=A.size();
	mat B(n,vec(n+1));
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)B[i][j]=A[i][j];
	for(int i=0;i<n;i++)B[i][n]=b[i];
	for(int i=0;i<n;i++){
		int pivot=i;
		for(int j=i;j<n;j++){
			if(ABS(B[j][i])>ABS(B[pivot][i]))pivot=j;
		}
		swap(B[i],B[pivot]);
		if(ABS(B[i][i])<EPS)return vec();
		for(int j=i+1;j<=n;j++)B[i][j]/=B[i][i];
		for(int j=0;j<n;j++){
			if(i!=j){
				for(int k=i+1;k<=n;k++)B[j][k]-=B[j][i]*B[i][k];
			}
		}
	}
	vec x(n);
	for(int i=0;i<n;i++)x[i]=B[i][n];
	return x;
}
double e[110];
int main(){
	for(int I=1;I<=100;I++){
		int a=I;
		mat A(a+1,vec(a+1));
		vec b(a+1);
		for(int i=0;i<=a;i++){
			A[i][i]=1;
			if(i<a){
				for(int j=1;j<=6;j++){
					int to=j+i;
					if(to>a)to=0;
					A[i][to]-=1.0/6;
				}
			}
			if(i<a)b[i]=1;
		}
		vec x=gauss_jordan(A,b);
		e[I]=x[0];
	}
	int T;scanf("%d",&T);
	while(T--){
		long long a;scanf("%lld",&a);
		if(a<=100){
			printf("%.18f\n",e[a]);
		}else{
			printf("%lld.6666666666666666\n",a+1);
		}
	}
}
0