結果
問題 | No.890 移調の限られた旋法 |
ユーザー |
![]() |
提出日時 | 2019-09-20 22:45:23 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 106 ms / 2,000 ms |
コード長 | 905 bytes |
コンパイル時間 | 933 ms |
コンパイル使用メモリ | 75,932 KB |
最終ジャッジ日時 | 2025-01-07 18:47:50 |
ジャッジサーバーID (参考情報) |
judge4 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> using namespace std; typedef long long ll; const int MAX = 1000010; ll inv[MAX],fact[MAX],fi[MAX],mod = 1000000007; void inverse(){ int i; inv[1] = 1; for(i=2;i<MAX;i++){ inv[i] = mod - (mod/i)*inv[mod%i]%mod; } fact[0] = fi[0] = 1; for(i=1;i<MAX;i++){ fact[i] = fact[i-1]*i%mod; fi[i] = (fi[i-1]*inv[i])%mod; } } ll comb(int n,int k){ if(n<0 || k<0 || n<k){ return 0; }else{ return fact[n]%mod*fi[k]%mod*fi[n-k]%mod; } } ll ans = 0,cnt[1000010] = {}; vector<ll> v; int main(){ int i,j,n,k; cin >> n >> k; for(i=1;i<n;i++){ if(n%i==0){ v.push_back(i); } } inverse(); for(i=0;i<v.size();i++){ ll x = v[i]; if(k%(n/x)!=0) continue; cnt[x] = comb(x,k/(n/x)); for(j=1;j<x;j++){ if(x%j==0){ cnt[x] -= cnt[j]; if(cnt[x]<0) cnt[x] += mod; } } (ans += cnt[x]) %= mod; } cout << ans << endl; }