結果

問題 No.1116 Cycles of Dense Graph
ユーザー mfbgjsczmfbgjscz
提出日時 2020-07-18 19:38:51
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 503 ms / 2,000 ms
コード長 3,616 bytes
コンパイル時間 2,288 ms
コンパイル使用メモリ 196,552 KB
実行使用メモリ 238,196 KB
最終ジャッジ日時 2023-08-20 20:59:19
合計ジャッジ時間 21,069 ms
ジャッジサーバーID
(参考情報)
judge11 / judge14
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 389 ms
237,824 KB
testcase_01 AC 432 ms
237,892 KB
testcase_02 AC 395 ms
237,848 KB
testcase_03 AC 396 ms
237,892 KB
testcase_04 AC 394 ms
237,840 KB
testcase_05 AC 408 ms
237,864 KB
testcase_06 AC 411 ms
237,976 KB
testcase_07 AC 415 ms
237,852 KB
testcase_08 AC 395 ms
237,840 KB
testcase_09 AC 400 ms
238,052 KB
testcase_10 AC 395 ms
237,920 KB
testcase_11 AC 407 ms
237,968 KB
testcase_12 AC 397 ms
237,836 KB
testcase_13 AC 392 ms
237,892 KB
testcase_14 AC 399 ms
237,832 KB
testcase_15 AC 402 ms
237,912 KB
testcase_16 AC 398 ms
237,888 KB
testcase_17 AC 378 ms
237,860 KB
testcase_18 AC 385 ms
237,828 KB
testcase_19 AC 371 ms
237,828 KB
testcase_20 AC 396 ms
237,864 KB
testcase_21 AC 414 ms
237,932 KB
testcase_22 AC 418 ms
238,100 KB
testcase_23 AC 394 ms
237,756 KB
testcase_24 AC 475 ms
237,928 KB
testcase_25 AC 390 ms
237,860 KB
testcase_26 AC 488 ms
237,936 KB
testcase_27 AC 400 ms
237,924 KB
testcase_28 AC 419 ms
238,076 KB
testcase_29 AC 402 ms
237,840 KB
testcase_30 AC 380 ms
237,792 KB
testcase_31 AC 376 ms
237,704 KB
testcase_32 AC 385 ms
237,852 KB
testcase_33 AC 396 ms
237,780 KB
testcase_34 AC 397 ms
237,776 KB
testcase_35 AC 396 ms
238,000 KB
testcase_36 AC 503 ms
238,196 KB
testcase_37 AC 463 ms
237,916 KB
testcase_38 AC 471 ms
237,904 KB
testcase_39 AC 429 ms
237,844 KB
testcase_40 AC 446 ms
237,980 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
/*#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/multiprecision/cpp_int.hpp>
*/#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")    
using namespace std;
/*namespace mp = boost::multiprecision;
using Bint = mp::cpp_int;
using Real = mp::number<mp::cpp_dec_float<1024>>;
*/#define lli long long int
#define uli unsigned long long int
#define INF 9999999999999999
#define rep(i,m,n) for(lli i = m;i < n;i++)
#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(),n.end())
#define Out(S) cout << S << endl
#define NeOut(S) cout << S
#define HpOut(S) cout << setprecision(50) << S << endl
#define Vec(K,L,N,S) vector<L> K(N,S)
#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))
#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))
#define pint pair<lli,lli>
#define paf(L,R) pair<L,R>
#define mod 998244353
#define MAX 10000000
#define ALL(a)  a.begin(),a.end()
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
long long fac[MAX], finv[MAX], inv[MAX];
void COMinit() {
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (int i = 2; i < MAX; i++){
        fac[i] = fac[i - 1] * i % mod;
        inv[i] = mod - inv[mod%i] * (mod / i) % mod;
        finv[i] = finv[i - 1] * inv[i] % mod;
    }
}
lli nCr(lli n, lli r){
    if (n < r) return 0;
    if (n < 0 || r < 0) return 0;
    return fac[n] * (finv[r] * finv[n - r] % mod) % mod;
}
lli modpow(lli n,lli k,lli m){
  if(k==0)return 1;
  else if(k%2==1)return modpow(n,k-1,m)*n%m;
  else{
    lli temp=modpow(n,k/2,m);
    return temp*temp%m;
  }
}
int main(){
  ios::sync_with_stdio(false);
  cin.tie(0);
  lli A,B,X,Y,Z,N,M,sum=0,flag=0;
  cin >> N >> M;
  COMinit();
  Vec(G,pint,0,pint(0,0));
  map<int,int>m1;
  while(M--){
    cin >> A >> B;
    G.pb(pint(A-1,B-1));
    m1[A-1]++;
    m1[B-1]++;
  }
  if(3>N){
    Out(0);
    return 0;
  }
  Vec(pk,pint,0,pint(0,0));
  for(auto v:m1)pk.pb(pint(pk.size(),v.first));
  map<int,int>m2;
  for(auto v:pk)m2[v.second]=v.first;
  M=G.size();
  rep(i,0,M){
    G[i].first=m2[G[i].first];
    G[i].second=m2[G[i].second];
  }
  DV(dp,lli,31,15,0);
  Vec(dp2,lli,N+1,1);
  rep(i,2,N+1)dp2[i]=i*(dp2[i-1])%mod;
  rep(i,3,N+1)sum=(sum+dp2[i-1]*nCr(N,i)%mod*inv[2])%mod;
  rep(i,1,N-1)dp[2][0]=(dp[2][0]+dp2[i]*nCr(N-2,i))%mod;
  rep(i,3,31)rep(j,0,N-i+1)rep(k,0,15)dp[i][k]=(dp[i][k]+dp2[j]*nCr(N-i,j)%mod*nCr(j+k,k)%mod)%mod;
  for(lli i=1;i<(1<<M);i++){
    DV(Gp,int,30,0,0);
    X=0,Y=0,Z=0;
    flag=1;
    Vec(dist,lli,30,-1);
    rep(j,0,M)if(i&(1<<j)){
      Gp[G[j].first].pb(G[j].second);
      Gp[G[j].second].pb(G[j].first);
      dist[G[j].first]=0;
      dist[G[j].second]=0;
      X++;  
    }
    rep(j,0,30)Y+=dist[j]+1;
    rep(j,0,30)if(Gp[j].size()>2)flag=0;
    if(!flag)continue;
    rep(j,0,30){
      if(!Gp[j].size())continue;
      if(dist[j])continue;
      dist[j]=1;
      queue<int>que;
      que.push(j);
      while(!que.empty()){
        auto v=que.front();que.pop();
        for(auto nv:Gp[v]){
          if(dist[nv])continue;
          dist[nv]=dist[v];
          que.push(nv);
        }
      }
      Z++;
    }
    if(Y-Z>=X){
      if(X==1)sum+=pow(-1,X%2)*dp[2][0];
      else sum+=pow(-1,X%2)*(dp2[Z-1]*modpow(2,Z-1,mod)%mod*dp[Y][Z-1]%mod);
    }
    else if(X==Y&&Z==1)sum+=pow(-1,X%2);
    sum=(sum+mod)%mod;
  }
  Out(sum%mod);
}
0