結果
問題 | No.890 移調の限られた旋法 |
ユーザー | tekihei2317 |
提出日時 | 2019-09-21 09:11:53 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 104 ms / 2,000 ms |
コード長 | 919 bytes |
コンパイル時間 | 1,641 ms |
コンパイル使用メモリ | 170,660 KB |
実行使用メモリ | 34,484 KB |
最終ジャッジ日時 | 2024-09-17 06:38:07 |
合計ジャッジ時間 | 5,086 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 40 ms
26,752 KB |
testcase_01 | AC | 41 ms
26,804 KB |
testcase_02 | AC | 37 ms
26,752 KB |
testcase_03 | AC | 42 ms
26,968 KB |
testcase_04 | AC | 44 ms
26,752 KB |
testcase_05 | AC | 42 ms
26,792 KB |
testcase_06 | AC | 46 ms
26,772 KB |
testcase_07 | AC | 44 ms
26,756 KB |
testcase_08 | AC | 42 ms
26,880 KB |
testcase_09 | AC | 42 ms
26,752 KB |
testcase_10 | AC | 49 ms
26,752 KB |
testcase_11 | AC | 44 ms
26,880 KB |
testcase_12 | AC | 44 ms
26,880 KB |
testcase_13 | AC | 97 ms
34,304 KB |
testcase_14 | AC | 92 ms
34,432 KB |
testcase_15 | AC | 88 ms
34,484 KB |
testcase_16 | AC | 92 ms
34,304 KB |
testcase_17 | AC | 96 ms
33,536 KB |
testcase_18 | AC | 101 ms
33,684 KB |
testcase_19 | AC | 63 ms
30,720 KB |
testcase_20 | AC | 57 ms
29,292 KB |
testcase_21 | AC | 50 ms
27,136 KB |
testcase_22 | AC | 83 ms
32,492 KB |
testcase_23 | AC | 87 ms
33,972 KB |
testcase_24 | AC | 69 ms
30,720 KB |
testcase_25 | AC | 49 ms
27,628 KB |
testcase_26 | AC | 88 ms
33,996 KB |
testcase_27 | AC | 77 ms
34,236 KB |
testcase_28 | AC | 67 ms
31,612 KB |
testcase_29 | AC | 59 ms
29,912 KB |
testcase_30 | AC | 103 ms
33,664 KB |
testcase_31 | AC | 73 ms
30,652 KB |
testcase_32 | AC | 104 ms
33,152 KB |
testcase_33 | AC | 95 ms
33,408 KB |
testcase_34 | AC | 92 ms
33,536 KB |
ソースコード
#include<bits/stdc++.h> #define int long long using namespace std; const int mod=1e9+7; const int sz=1000010; int fact[sz],inv[sz],ifact[sz]; void make(){ fact[0]=fact[1]=inv[1]=ifact[0]=ifact[1]=1; for(int i=2;i<sz;i++){ fact[i]=fact[i-1]*i%mod; inv[i]=inv[mod%i]*(mod-mod/i)%mod; ifact[i]=ifact[i-1]*inv[i]%mod; } } int comb(int n,int k){ if(n<k) return 0; return fact[n]*ifact[k]%mod*ifact[n-k]%mod; } signed main() { cin.tie(0); ios::sync_with_stdio(false); make(); int N,K; cin>>N>>K; vector<int> dp(N); for(int i=1;i<N;i++) if(N%i==0){ int j=N/i; if(K%j==0){ dp[i]=comb(i,K/j); } } int ans=0; for(int i=1;i<N;i++) if(N%i==0){ ans=(ans+dp[i])%mod; for(int j=i*2;j<N;j+=i) if(N%j==0){ dp[j]=(dp[j]-dp[i]+mod)%mod; } } cout<<ans<<endl; return 0; }