結果

問題 No.3674 Zero Sum Game
コンテスト
ユーザー butsurizuki
提出日時 2026-09-04 23:14:02
言語 C++23
(gcc 15.3.0 + boost 1.92.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
TLE  
実行時間 -
コード長 1,635 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,972 ms
コンパイル使用メモリ 345,080 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2026-09-04 23:14:10
合計ジャッジ時間 6,884 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge5_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other TLE * 1 -- * 38
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>

using namespace std;
using ld=long double;

void fix(vector<ld> &v){
  ld s=0.0;
  for(auto &nx : v){
    s+=nx;
  }
  s=sqrt(s);
  for(auto &nx : v){
    nx/=s;
  }
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  int t;
  cin >> t;
  while(t--){
    int n,m;
    cin >> n >> m;
    vector<vector<ld>> a(n,vector<ld>(m,0));
    for(auto &nx : a){
      for(auto &ny : nx){cin >> ny;}
    }
    ld alpha=1.0;
    ld beta=0.99;
    vector<ld> alice(n,((ld)1.0)/((ld)n));
    vector<ld> bob(m,((ld)1.0)/((ld)m));
    for(int tr=0;tr<10000;tr++){
      {
        // alice
        vector<ld> as(n,0.0);
        for(int i=0;i<n;i++){
          for(int j=0;j<m;j++){
            as[i]+=bob[j]*a[i][j];
          }
        }
        int tg=0;
        for(int i=0;i<n;i++){
          if(as[tg]<as[i]){tg=i;}
        }
        alice[tg]+=(1.0-alice[tg])*alpha;
        for(int i=0;i<n;i++){
          if(i!=tg){alice[i]*=(1.0-alpha);}
        }
      }
      fix(alice);
      {
        // bob
        vector<ld> bs(m,0.0);
        for(int i=0;i<n;i++){
          for(int j=0;j<m;j++){
            bs[j]+=alice[i]*a[i][j];
          }
        }
        int tg=0;
        for(int i=0;i<m;i++){
          if(bs[tg]>bs[i]){tg=i;}
        }
        bob[tg]+=(1.0-bob[tg])*alpha;
        for(int i=0;i<m;i++){
          if(i!=tg){bob[i]*=(1.0-alpha);}
        }
      }
      alpha*=beta;
      fix(bob);
    }
    ld res=0.0;
    for(int i=0;i<n;i++){
      for(int j=0;j<m;j++){
        res+=alice[i]*bob[j]*a[i][j];
      }
    }
    cout << fixed << setprecision(20) << res << "\n";
  }
  return 0;
}
0