結果

問題 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 #
プレゼンテーションモードにする

#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);
        }
    }
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0