結果
問題 | 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;}