結果
問題 | No.890 移調の限られた旋法 |
ユーザー | mtsd |
提出日時 | 2019-09-20 21:56:56 |
言語 | C++14 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 67 ms / 2,000 ms |
コード長 | 1,892 bytes |
コンパイル時間 | 1,240 ms |
コンパイル使用メモリ | 118,008 KB |
実行使用メモリ | 32,444 KB |
最終ジャッジ日時 | 2023-10-12 19:00:52 |
合計ジャッジ時間 | 4,591 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge15 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 46 ms
28,144 KB |
testcase_01 | AC | 42 ms
28,084 KB |
testcase_02 | AC | 45 ms
28,084 KB |
testcase_03 | AC | 46 ms
28,224 KB |
testcase_04 | AC | 45 ms
28,168 KB |
testcase_05 | AC | 44 ms
28,156 KB |
testcase_06 | AC | 45 ms
28,012 KB |
testcase_07 | AC | 45 ms
28,028 KB |
testcase_08 | AC | 45 ms
28,148 KB |
testcase_09 | AC | 48 ms
28,096 KB |
testcase_10 | AC | 48 ms
28,152 KB |
testcase_11 | AC | 43 ms
28,080 KB |
testcase_12 | AC | 47 ms
28,092 KB |
testcase_13 | AC | 56 ms
32,344 KB |
testcase_14 | AC | 51 ms
28,076 KB |
testcase_15 | AC | 57 ms
28,056 KB |
testcase_16 | AC | 56 ms
30,148 KB |
testcase_17 | AC | 55 ms
28,228 KB |
testcase_18 | AC | 63 ms
30,180 KB |
testcase_19 | AC | 55 ms
30,100 KB |
testcase_20 | AC | 48 ms
28,068 KB |
testcase_21 | AC | 50 ms
28,172 KB |
testcase_22 | AC | 51 ms
28,084 KB |
testcase_23 | AC | 46 ms
28,044 KB |
testcase_24 | AC | 49 ms
28,032 KB |
testcase_25 | AC | 50 ms
28,248 KB |
testcase_26 | AC | 53 ms
28,172 KB |
testcase_27 | AC | 50 ms
28,160 KB |
testcase_28 | AC | 48 ms
28,068 KB |
testcase_29 | AC | 49 ms
28,024 KB |
testcase_30 | AC | 57 ms
28,076 KB |
testcase_31 | AC | 52 ms
28,352 KB |
testcase_32 | AC | 67 ms
30,280 KB |
testcase_33 | AC | 62 ms
32,444 KB |
testcase_34 | AC | 56 ms
30,216 KB |
ソースコード
#include <algorithm> #include <bitset> #include <cassert> #include <chrono> #include <climits> #include <cmath> #include <complex> #include <cstring> #include <deque> #include <functional> #include <iostream> #include <list> #include <map> #include <numeric> #include <queue> #include <random> #include <set> #include <stack> #include <unordered_map> #include <unordered_set> #include <vector> #include <cstdint> using namespace std; typedef long long ll; #define MP make_pair #define PB push_back #define inf 1000000007 #define mod 1000000007 #define rep(i,n) for(int i = 0; i < (int)(n); ++i) template<typename T> vector<T> divisor(T n) { vector<T> res; for(T i=1;(long long)i*i<=n;i++){ if(n%i==0){ res.push_back(i); if(i != n/i){ res.push_back(n/i); } } } sort(res.rbegin(),res.rend()); return res; } #define MAX_N 2000010 #define MOD 1000000007 int inv[MAX_N],fac[MAX_N],finv[MAX_N]; void make() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for(int i=2;i<MAX_N;i++){ inv[i] = MOD - (long long)inv[MOD%i] * (MOD/i) % MOD; fac[i] = (long long)fac[i-1] * i % MOD; finv[i] = (long long)finv[i-1] * inv[i] % MOD; } } int comb(int a,int b) { if(a < b){ return 0; } return fac[a] * ((long long)finv[b] * finv[a-b] % MOD) % MOD; } ll dp[2000000]={}; int main(){ make(); int n,k; cin >> n >> k; vector<int> res = divisor(n); ll sm = 0; for(int x:res){ if(x==1)continue; if(k%x==0){ int c = k/x; int p = n/x; dp[x] = comb(p,c); for(int i=x+x;i<=n;i+=x){ dp[x] += mod-dp[i]; dp[x] %= mod; } sm += dp[x]; sm %= mod; } } cout << sm << endl; return 0; }