結果
問題 | No.1659 Product of Divisors |
ユーザー | ytft |
提出日時 | 2021-08-30 12:52:30 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 25 ms / 2,000 ms |
コード長 | 1,027 bytes |
コンパイル時間 | 6,541 ms |
コンパイル使用メモリ | 413,488 KB |
実行使用メモリ | 4,500 KB |
最終ジャッジ日時 | 2023-08-15 10:12:49 |
合計ジャッジ時間 | 7,808 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,384 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 1 ms
4,376 KB |
testcase_03 | AC | 1 ms
4,500 KB |
testcase_04 | AC | 18 ms
4,380 KB |
testcase_05 | AC | 12 ms
4,376 KB |
testcase_06 | AC | 2 ms
4,380 KB |
testcase_07 | AC | 2 ms
4,380 KB |
testcase_08 | AC | 24 ms
4,380 KB |
testcase_09 | AC | 24 ms
4,376 KB |
testcase_10 | AC | 24 ms
4,376 KB |
testcase_11 | AC | 25 ms
4,380 KB |
testcase_12 | AC | 24 ms
4,380 KB |
testcase_13 | AC | 2 ms
4,380 KB |
testcase_14 | AC | 23 ms
4,380 KB |
testcase_15 | AC | 2 ms
4,380 KB |
testcase_16 | AC | 2 ms
4,376 KB |
testcase_17 | AC | 2 ms
4,376 KB |
testcase_18 | AC | 2 ms
4,380 KB |
testcase_19 | AC | 2 ms
4,380 KB |
testcase_20 | AC | 13 ms
4,380 KB |
testcase_21 | AC | 19 ms
4,380 KB |
testcase_22 | AC | 18 ms
4,376 KB |
testcase_23 | AC | 22 ms
4,376 KB |
testcase_24 | AC | 13 ms
4,384 KB |
ソースコード
#include <bits/stdc++.h> #include <boost/multiprecision/cpp_int.hpp> using namespace std; namespace mp = boost::multiprecision; int main(){ mp::cpp_int mod=1000000007; mp::cpp_int N,K; cin>>N>>K; vector<int> prime(0); int upper=(int)sqrt(N)+1; for(int i=2;i<=upper;++i){ if(N%i==0){ prime.push_back(0); while(N%i==0){ N/=i; ++prime[prime.size()-1]; } } } if(N>1){ prime.push_back(1); } mp::cpp_int ans=1; for(int i:prime){ mp::cpp_int temp1=1,temp2=1; for(int j=1;j<=i;++j){ temp1=(temp1*(j+K))%mod; temp2=(temp2*j)%mod; } mp::cpp_int remain=mod-2; mp::cpp_int d=1; while(remain){ if(remain%2){ d=(d*temp2)%mod; } remain/=2; temp2=(temp2*temp2)%mod; } temp1=(temp1*d)%mod; ans=(ans*temp1)%mod; } cout<<((int)ans)<<endl; }