結果

問題 No.890 移調の限られた旋法
ユーザー kyort0nkyort0n
提出日時 2019-09-20 22:08:05
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 279 ms / 2,000 ms
コード長 2,925 bytes
コンパイル時間 2,071 ms
コンパイル使用メモリ 180,624 KB
実行使用メモリ 237,920 KB
最終ジャッジ日時 2024-09-14 17:49:30
合計ジャッジ時間 10,963 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 222 ms
237,788 KB
testcase_01 AC 222 ms
237,860 KB
testcase_02 AC 279 ms
237,824 KB
testcase_03 AC 220 ms
237,824 KB
testcase_04 AC 222 ms
237,892 KB
testcase_05 AC 221 ms
237,696 KB
testcase_06 AC 220 ms
237,824 KB
testcase_07 AC 220 ms
237,688 KB
testcase_08 AC 221 ms
237,824 KB
testcase_09 AC 221 ms
237,824 KB
testcase_10 AC 221 ms
237,824 KB
testcase_11 AC 222 ms
237,680 KB
testcase_12 AC 222 ms
237,824 KB
testcase_13 AC 221 ms
237,824 KB
testcase_14 AC 221 ms
237,696 KB
testcase_15 AC 222 ms
237,680 KB
testcase_16 AC 225 ms
237,816 KB
testcase_17 AC 222 ms
237,876 KB
testcase_18 AC 220 ms
237,824 KB
testcase_19 AC 221 ms
237,700 KB
testcase_20 AC 225 ms
237,696 KB
testcase_21 AC 222 ms
237,788 KB
testcase_22 AC 220 ms
237,776 KB
testcase_23 AC 222 ms
237,824 KB
testcase_24 AC 221 ms
237,696 KB
testcase_25 AC 222 ms
237,824 KB
testcase_26 AC 222 ms
237,696 KB
testcase_27 AC 222 ms
237,780 KB
testcase_28 AC 223 ms
237,696 KB
testcase_29 AC 223 ms
237,696 KB
testcase_30 AC 223 ms
237,696 KB
testcase_31 AC 223 ms
237,836 KB
testcase_32 AC 223 ms
237,920 KB
testcase_33 AC 222 ms
237,716 KB
testcase_34 AC 222 ms
237,824 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll, ll> l_l;
typedef pair<int, int> i_i;
template<class T>
inline bool chmax(T &a, T b) {
    if(a < b) {
        a = b;
        return true;
    }
    return false;
}

template<class T>
inline bool chmin(T &a, T b) {
    if(a > b) {
        a = b;
        return true;
    }
    return false;
}

#define EPS (1e-7)
#define INF (1e9)
#define PI (acos(-1))
const ll mod = 1000000007;
ll inv[10000100];
ll FactorialInv[10000100];
ll Factorial[10000100];
ll beki(ll a, ll b){
    if(b == 0){
        return 1;
    }
    ll ans = beki(a, b / 2);
    ans = ans * ans % mod;
    if(b % 2 == 1){
        ans = ans * a % mod;
    }
    return ans;
}
void init_combination(){
    const int MAX = 10000002;
    Factorial[0] = 1;
    inv[0] = 1;
    for(int i = 1; i <= MAX; i++){
        Factorial[i] = Factorial[i - 1] * i % mod;
    }
    FactorialInv[MAX] = beki(Factorial[MAX], mod - 2);
    for(ll i = MAX - 1; i >= 0; i--) {
        FactorialInv[i] = FactorialInv[i+1] * (i+1) % mod;
    }
    for(int i = 1; i <= MAX; i++) {
        inv[i] = FactorialInv[i] * Factorial[i-1] % mod;
    }
}
ll combination(ll a, ll b){
    if((a == b) || (b == 0)){
        return 1;
    }
    if(a < b) return 0;
    ll ans = Factorial[a] * FactorialInv[b] % mod;
    ans = ans * FactorialInv[a - b] % mod;
    return ans;
}
ll N, K;
map<ll, ll> mp;

vector<ll> f(ll x) {
    vector<ll> v;
    v.push_back(1);
    for(ll i = 2; i * i <= N; i++) {
        if(x % i != 0) continue;
        v.push_back(i);
        if(i * i != N) v.push_back(N / i);
    }
    sort(v.begin(), v.end());
    return v;
}

vector<ll> g(ll x) {
    vector<ll> ret;
    for(ll i = 2; i * i <= x; i++) {
        if(x % i == 0) {
            ret.push_back(i);
            while(x % i == 0) x /= i;
        }
    }
    if(x != 1) ret.push_back(x);
    return ret;
}

int counter(ll a) {
    if(a == 0) return 0;
    return (a & 1) + counter(a >> 1);
}

int main() {
    init_combination();
    cin >> N >> K;
    if(K == 1) {
        cout << 0 << endl;
        return 0;
    }
    vector<ll> divisors = f(N);
    ll ans = 0;
    for(auto range : divisors) {
        ll val = N / range;
        if(K % val != 0) continue;
        ll now = combination(range,  K / val);
        mp[range] = now;
        //cerr << val << " " << range << " " << now << endl;
        vector<ll> p = g(range);
        for(ll bit = 1; bit < (1LL << (int)p.size()); bit++) {
            ll nowrange = range;
            for(int i = 0; i < p.size(); i++) {
                if(bit & (1LL << i)) nowrange /= p[i];
            }
            if(counter(bit) & 1LL) now -= mp[nowrange];
            else now += mp[nowrange];
            now += mod;
            now %= mod;
        }
        now %= mod;
        now += mod;
        now %= mod;
        ans = (ans + now) % mod;
    }
    cout << ans << endl;
    return 0;
}
0