結果

問題 No.890 移調の限られた旋法
ユーザー Y17Y17
提出日時 2019-09-20 22:49:20
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 27 ms / 2,000 ms
コード長 1,312 bytes
コンパイル時間 1,346 ms
コンパイル使用メモリ 161,300 KB
実行使用メモリ 18,944 KB
最終ジャッジ日時 2024-09-14 18:56:28
合計ジャッジ時間 2,995 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 10 ms
11,136 KB
testcase_01 AC 11 ms
11,136 KB
testcase_02 AC 10 ms
11,392 KB
testcase_03 AC 11 ms
11,264 KB
testcase_04 AC 10 ms
11,264 KB
testcase_05 AC 10 ms
11,264 KB
testcase_06 AC 11 ms
11,264 KB
testcase_07 AC 11 ms
11,392 KB
testcase_08 AC 10 ms
11,264 KB
testcase_09 AC 11 ms
11,264 KB
testcase_10 AC 10 ms
11,264 KB
testcase_11 AC 11 ms
11,264 KB
testcase_12 AC 10 ms
11,264 KB
testcase_13 AC 26 ms
18,892 KB
testcase_14 AC 27 ms
18,944 KB
testcase_15 AC 26 ms
18,944 KB
testcase_16 AC 25 ms
18,944 KB
testcase_17 AC 26 ms
18,304 KB
testcase_18 AC 26 ms
18,176 KB
testcase_19 AC 19 ms
15,360 KB
testcase_20 AC 16 ms
13,952 KB
testcase_21 AC 11 ms
11,776 KB
testcase_22 AC 21 ms
17,152 KB
testcase_23 AC 25 ms
18,688 KB
testcase_24 AC 18 ms
15,488 KB
testcase_25 AC 13 ms
12,288 KB
testcase_26 AC 26 ms
18,816 KB
testcase_27 AC 26 ms
18,944 KB
testcase_28 AC 20 ms
16,256 KB
testcase_29 AC 17 ms
14,336 KB
testcase_30 AC 25 ms
18,304 KB
testcase_31 AC 18 ms
15,360 KB
testcase_32 AC 23 ms
17,664 KB
testcase_33 AC 24 ms
18,176 KB
testcase_34 AC 25 ms
18,048 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

#define int long long
#define MOD 1000000007

int pow_mod(int n, int m){
    int ans = 1;
    while(m > 0){
        if(m & 1) ans = (ans * n) % MOD;
        n = (n * n) % MOD;
        m >>= 1;
    }
    return ans;
}

int kai[1000010];
void init(){
    kai[0] = 1;
    for(int i = 1;i < 1000010;i++){
        kai[i] = (kai[i-1] * i) % MOD;
    }
    return;
}

int combi(int n, int r){
    return ((kai[n] * pow_mod(kai[r], MOD-2) % MOD) * pow_mod(kai[n-r], MOD-2)) % MOD;
}

int so[1000010] = {};

int n, k;

vector<int> vec;

int dfs(int i, int num, int cnt){
    if(i == vec.size()){
        if(cnt == 0) return 0;
        if(cnt % 2 == 1) return combi(n/num, k/num);
        else return -combi(n/num, k/num);
    }

    int ans1 = dfs(i+1, num, cnt);
    int ans2 = dfs(i+1, num*vec[i], cnt+1);
    int ans = (MOD*2ll + ans1 + ans2) % MOD;
    return ans;
}

signed main(){
    cin >> n >> k;

    init();

    for(int i = 2;i <= n;i++){
        if(so[i] == 0){
            for(int j = i*2;j < n;j+=i){
                so[j] = 1;
            }
        }
    }

    int ans = 0;
    for(int i = 2;i <= n;i++){
        if(so[i] == 0 && n % i == 0 && k % i == 0){
            vec.push_back(i);
        }
    }

    cout << dfs(0, 1, 0) << endl;
    return 0;
}
0