結果
問題 | No.890 移調の限られた旋法 |
ユーザー | otera |
提出日時 | 2019-09-21 20:44:18 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 95 ms / 2,000 ms |
コード長 | 2,356 bytes |
コンパイル時間 | 1,174 ms |
コンパイル使用メモリ | 110,796 KB |
実行使用メモリ | 36,480 KB |
最終ジャッジ日時 | 2024-09-19 03:06:48 |
合計ジャッジ時間 | 4,394 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 77 ms
36,224 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 78 ms
36,216 KB |
testcase_03 | AC | 79 ms
36,224 KB |
testcase_04 | AC | 73 ms
36,344 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 76 ms
36,272 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 77 ms
36,300 KB |
testcase_10 | AC | 2 ms
5,376 KB |
testcase_11 | AC | 2 ms
5,376 KB |
testcase_12 | AC | 77 ms
36,352 KB |
testcase_13 | AC | 88 ms
36,224 KB |
testcase_14 | AC | 13 ms
5,376 KB |
testcase_15 | AC | 90 ms
36,480 KB |
testcase_16 | AC | 91 ms
36,224 KB |
testcase_17 | AC | 90 ms
36,352 KB |
testcase_18 | AC | 89 ms
36,480 KB |
testcase_19 | AC | 81 ms
36,480 KB |
testcase_20 | AC | 6 ms
5,376 KB |
testcase_21 | AC | 3 ms
5,376 KB |
testcase_22 | AC | 84 ms
36,224 KB |
testcase_23 | AC | 13 ms
5,376 KB |
testcase_24 | AC | 87 ms
36,352 KB |
testcase_25 | AC | 4 ms
5,376 KB |
testcase_26 | AC | 95 ms
36,352 KB |
testcase_27 | AC | 91 ms
36,224 KB |
testcase_28 | AC | 9 ms
5,376 KB |
testcase_29 | AC | 76 ms
36,352 KB |
testcase_30 | AC | 85 ms
36,480 KB |
testcase_31 | AC | 83 ms
36,260 KB |
testcase_32 | AC | 93 ms
36,480 KB |
testcase_33 | AC | 93 ms
36,352 KB |
testcase_34 | AC | 91 ms
36,480 KB |
ソースコード
#include<iostream> #include<string> #include<cstdio> #include<cstring> #include<vector> #include<cmath> #include<algorithm> #include<functional> #include<iomanip> #include<queue> #include<deque> #include<ciso646> #include<random> #include<map> #include<set> #include<complex> #include<bitset> #include<stack> #include<unordered_map> #include<utility> #include<cassert> using namespace std; typedef long long ll; typedef unsigned long long ul; typedef unsigned int ui; typedef long double ld; const int inf=1e9+7; const ll longinf=1LL<<60 ; const ll mod=1e9+7 ; #define rep(i,n) for(int i=0;i<n;i++) #define per(i,n) for(int i=n-1;i>=0;i--) #define Rep(i,sta,n) for(int i=sta;i<n;i++) #define rep1(i,n) for(int i=1;i<=n;i++) #define per1(i,n) for(int i=n;i>=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef complex<ld> Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair<int, int> P; typedef pair<ld, ld> LDP; typedef pair<ll, ll> LP; #define fr first #define sc second #define all(c) c.begin(),c.end() #define pb push_back #define debug(x) cout << #x << " = " << (x) << endl; template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } #define int long long ll comb(ll N_, ll C_) { const int NUM_=1400001; static ll fact[NUM_+1],factr[NUM_+1],inv[NUM_+1]; if (fact[0]==0) { inv[1]=fact[0]=factr[0]=1; for (int i=2;i<=NUM_;++i) inv[i] = inv[mod % i] * (mod - mod / i) % mod; for (int i=1;i<=NUM_;++i) fact[i]=fact[i-1]*i%mod, factr[i]=factr[i-1]*inv[i]%mod; } if(C_<0 || C_>N_) return 0; return factr[C_]*fact[N_]%mod*factr[N_-C_]%mod; } void solve() { ll n, k; cin >> n >> k; ll f[1001000]; rep1(i, n - 1) { if(n % i == 0) { ll p = n / i; if(k % p == 0) f[i] = comb(i, k / p); } } ll ans = 0; rep1(i, n - 1) { if(f[i]) { (ans += f[i]) %= mod; for(int j = i * 2; j <= n; j+= i) { if(f[j]) { (f[j] += mod - f[i]) %= mod; } } } } cout << ans << endl; } signed main() { ios::sync_with_stdio(false); cin.tie(0); //cout << fixed << setprecision(10); //init(); solve(); //cout << "finish" << endl; return 0; }