結果

問題 No.309 シャイな人たち (1)
ユーザー piyoko_212
提出日時 2015-12-25 23:34:56
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 321 ms / 4,000 ms
コード長 1,548 bytes
コンパイル時間 498 ms
コンパイル使用メモリ 44,748 KB
実行使用メモリ 18,372 KB
最終ジャッジ日時 2024-09-19 00:05:02
合計ジャッジ時間 4,956 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 13
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:15:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   15 |         int a,b;scanf("%d%d",&a,&b);
      |                 ~~~~~^~~~~~~~~~~~~~
main.cpp:16:54: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   16 |         for(int i=0;i<a;i++)for(int j=0;j<b;j++)scanf("%d",&p[i][j]);
      |                                                 ~~~~~^~~~~~~~~~~~~~~
main.cpp:17:54: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   17 |         for(int i=0;i<a;i++)for(int j=0;j<b;j++)scanf("%d",&q[i][j]);
      |                                                 ~~~~~^~~~~~~~~~~~~~~

ソースコード

diff #
プレゼンテーションモードにする

#include<stdio.h>
#include<algorithm>
#include<queue>
using namespace std;
int p[20][20];
int q[20][20];
double dp[13][1<<11];
int st[20];
double pr[1<<11];
int v[20];
int si[1<<22];
int mask[1<<11];
double EPS=1e-14;
int main(){
    int a,b;scanf("%d%d",&a,&b);
    for(int i=0;i<a;i++)for(int j=0;j<b;j++)scanf("%d",&p[i][j]);
    for(int i=0;i<a;i++)for(int j=0;j<b;j++)scanf("%d",&q[i][j]);
    for(int i=0;i<(1<<b);i++){
        for(int j=0;j<b;j++){
            if(!(i&(1<<j))){
                mask[i]+=(1<<(j*2))+(1<<(j*2+1));
            }
        }
        //printf("%d: %d\n",i,mask[i]);
    }
    for(int i=0;i<(1<<(b*2));i++){
        for(int j=0;j<b;j++){
            st[j]=3-(i>>(j*2))%4;
        }
        for(int j=0;j<b-1;j++){
            if(st[j]>=3)st[j+1]++;
        }
        for(int j=b-1;j>0;j--){
            if(st[j]>=3)st[j-1]++;
        }
        int to=0;
        for(int j=0;j<b;j++)if(st[j]>=3)to+=(1<<j);
        si[i]=to;
        
    }
    dp[0][0]=1;
    for(int i=0;i<a;i++){
        for(int j=0;j<(1<<b);j++){
            double th=1;
            for(int k=0;k<b;k++){
                if(j&(1<<k)){
                    th=th*p[i][k]/100;
                }else{
                    th=th*(100-p[i][k])/100;
                }
            }
            pr[j]=th;
        }
        for(int k=0;k<(1<<b);k++){
            if(dp[i][k]<EPS)continue;
            int req=0;
            for(int l=b-1;l>=0;l--){
                req*=4;
                int tt=q[i][l];
                if(q[i][l]&&(k&(1<<l)))tt--;
                if(tt==4)tt--;
                req+=tt;
            }
            for(int l=0;l<(1<<b);l++){
                if(pr[l]<EPS)continue;
                dp[i+1][si[req|mask[l]]]+=pr[l]*dp[i][k];
            }
        }
    }
    double ret=0;
    for(int i=1;i<=a;i++){
        for(int j=0;j<(1<<b);j++){
    //      if(dp[i][j]>EPS)printf("%d %d: %f\n",i,j,dp[i][j]);
            ret+=dp[i][j]*(__builtin_popcount(j));
        }
    }
    printf("%.12f\n",ret);
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0