結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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 4000000000000000001
mint 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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0