結果

問題 No.75 回数の期待値の問題
コンテスト
ユーザー piyoko_212
提出日時 2015-12-28 16:49:43
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 5 ms / 5,000 ms
コード長 1,029 bytes
コンパイル時間 380 ms
コンパイル使用メモリ 47,684 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-19 07:36:13
合計ジャッジ時間 1,170 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:34:20: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   34 |         int a;scanf("%d",&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;
}
int main(){
	int a;scanf("%d",&a);
	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);
	printf("%.12f\n",x[0]);
}
0