結果
問題 | No.2365 Present of good number |
ユーザー |
![]() |
提出日時 | 2023-06-30 22:00:33 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 1,183 bytes |
コンパイル時間 | 4,133 ms |
コンパイル使用メモリ | 250,372 KB |
最終ジャッジ日時 | 2025-02-15 03:52:49 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 39 |
ソースコード
#include <stdio.h>#include <atcoder/all>#include <bits/stdc++.h>using namespace std;using namespace atcoder;using mint = modint1000000007;#define rep(i,n) for (int i = 0; i < (n); ++i)#define Inf32 1000000001#define Inf64 4000000000000000001mint get(long long p,long long b,long long K){if(p==1)return 1;//cout<<p<<','<<b<<','<<K<<endl;//if(mp.count({p,b,K}))return mp[{p,b,K}];if(K==0)return mint(p).pow(b);if(p==2||p==3){if(K%2==1){if(p==2){return get(3,b,K-1);}else{return get(2,b*2%1000000006,K-1);}}else{//cout<<p<<','<<K<<','<<b<<endl;mint ret = p;ret = ret.pow(pow_mod(2,K/2,1000000006)*b);return ret;}}p++;{mint ret = 1;for(long long i=2;i*i<=p;i++){if(p%i==0){long long c = 0;while(p%i==0){p /= i;c++;}ret *= get(i,c*b%1000000006,K-1);}}if(p!=1){ret *= get(p,b,K-1);}return ret;}}int main(){int N;long long K;cin>>N>>K;mint ans = 1;for(long long i=2;i<=N;i++){if(N%i==0){int c = 0;while(N%i==0){N /= i;c++;}ans *= get(i,c,K);}}cout<<ans.val()<<endl;return 0;}