結果
問題 |
No.1310 量子アニーリング
|
ユーザー |
|
提出日時 | 2020-12-07 16:25:56 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 87 ms / 2,000 ms |
コード長 | 1,238 bytes |
コンパイル時間 | 2,326 ms |
コンパイル使用メモリ | 198,040 KB |
最終ジャッジ日時 | 2025-01-16 18:58:13 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
// #pragma GCC optimize("Ofast") // #pragma GCC optimize("unroll-loops") // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native") #include <bits/stdc++.h> using namespace std; typedef long long int ll; constexpr ll mod=998244353; ll mod_pow(ll a,ll b){ a%=mod; if(b==0)return 1; if(b==1)return a; ll res=mod_pow(a,b/2)%mod; res*=res; res%=mod; if(b%2)res*=a; return res%mod; } struct perm{ private: int sz; vector<ll> p,invp; public: perm(int n){ sz=n+1; p.resize(sz),invp.resize(sz); p[0]=1; for(int i=1;i<=sz-1;i++){ p[i]=p[i-1]*i%mod; } invp[sz-1]=mod_pow(p[sz-1],mod-2); for(int i=sz-2;i>=0;i--){ invp[i]=invp[i+1]*(i+1)%mod; } } ll comb(ll x,ll y){ if(x<y||y<0)return 0; return (p[x]*invp[x-y]%mod)*invp[y]%mod; } }; perm p(1<<20); int main(){ cin.tie(nullptr); ios::sync_with_stdio(false); ll n; cin >> n; ll res=mod_pow(2,n)*2%mod; for(int i=2;i<=n;i+=2){ // 連結成分は何個か? ll k=abs(n-2*i); ll a=p.comb(n,i)*2%mod; res+=a*mod_pow(2,k)%mod; res%=mod; } cout << res << endl; }