結果
問題 | No.1396 Giri |
ユーザー |
![]() |
提出日時 | 2021-02-14 21:45:14 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 19 ms / 2,000 ms |
コード長 | 1,121 bytes |
コンパイル時間 | 1,749 ms |
コンパイル使用メモリ | 168,648 KB |
実行使用メモリ | 15,064 KB |
最終ジャッジ日時 | 2024-07-22 09:14:39 |
合計ジャッジ時間 | 2,497 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,n) for(int i=0;i<n;++i) #define rep1(i,n) for(int i=1;i<=n;++i) #define debug(output) cout<<#output<<"= "<<output<<endl using lint = long long; typedef pair<int,int> P; const lint linf=1e18+7; const lint inf=1e9+7; const int MOD=998244353; long long myPow(long long x,long long n,long long rest){ long long res=1; long long cur=x; while(n>0){ if(n&1){ res*=cur; res%=rest; } n/=2;cur*=cur; cur%=rest; } return res; } signed main(){ int n;cin>>n; vector<int> memo(n+5,0); auto cnt=memo; vector<int> primes(n+5,0); for(int i=2;i<=n;++i){ if(memo[i]==0){ primes[i]=1; for(int j=1;j*i<=n;++j){ memo[i*j]=1; } } } for(int i=n;i>1;--i){ if(primes[i]){ primes[i]=0; break; } } lint res=1; rep1(i,n){ if(primes[i]>0){ int m=0; lint cur=1; rep1(j,n){ cur*=i; if(cur>n)break; ++m; } res*=myPow(i,m,MOD); res%=MOD; } } cout<<res<<"\n"; return 0; }