結果
問題 | No.2435 Order All Company |
ユーザー |
|
提出日時 | 2023-08-26 01:32:37 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,170 bytes |
コンパイル時間 | 4,723 ms |
コンパイル使用メモリ | 260,056 KB |
最終ジャッジ日時 | 2025-02-16 14:32:53 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 28 WA * 8 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using ll = long long; #define rep(i,n) for(int i=0;i<(int)(n);i++) using mint = atcoder::modint998244353; int n; using mat=vector<vector<mint>>; void tri(mat &a){ rep(i,n-1){ if(a.at(i).at(i)==0){ for(int j=i+1;j<n;j++){ if(a.at(j).at(i)!=0){ rep(k,n) swap(a.at(i).at(k),a.at(j).at(k)); rep(k,n) a.at(i).at(k)*=-1; break; } } if(a.at(i).at(i)==0) continue; } rep(j,n){ if(i==j) continue; mint p=a.at(j).at(i)/a.at(i).at(i); rep(k,n) a.at(j).at(k)-=a.at(i).at(k)*p; } } } int main(){ int k; cin>>n>>k; vector<mat> g(k,mat(n,vector<mint>(n,0))); rep(i,k){ int t; cin>>t; rep(lp,t){ int a,b; cin>>a>>b; a--; b--; g.at(i).at(a).at(a)++; g.at(i).at(b).at(b)++; g.at(i).at(a).at(b)--; g.at(i).at(b).at(a)--; } } mint ans=0; rep(bit,(1<<k)){ mat ad(n,vector<mint>(n,0)); rep(i,k){ if(bit>>i&1){ rep(j,n) rep(k,n) ad.at(j).at(k)+=g.at(i).at(j).at(k); } } tri(ad); mint tmp=1; rep(i,n-1) tmp*=ad.at(i).at(i); if(__builtin_popcount(bit)%2) ans+=tmp; else ans-=tmp; } cout<<ans.val()<<endl; }