結果

問題 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,442 ms
コンパイル使用メモリ 147,388 KB
実行使用メモリ 19,116 KB
最終ジャッジ日時 2023-10-12 20:35:43
合計ジャッジ時間 3,147 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 11 ms
12,860 KB
testcase_01 AC 11 ms
12,848 KB
testcase_02 AC 11 ms
12,920 KB
testcase_03 AC 11 ms
12,816 KB
testcase_04 AC 11 ms
12,824 KB
testcase_05 AC 11 ms
12,856 KB
testcase_06 AC 11 ms
12,844 KB
testcase_07 AC 10 ms
12,948 KB
testcase_08 AC 11 ms
12,876 KB
testcase_09 AC 11 ms
12,856 KB
testcase_10 AC 11 ms
12,872 KB
testcase_11 AC 11 ms
12,844 KB
testcase_12 AC 10 ms
12,856 KB
testcase_13 AC 27 ms
18,988 KB
testcase_14 AC 27 ms
18,944 KB
testcase_15 AC 27 ms
18,944 KB
testcase_16 AC 27 ms
19,076 KB
testcase_17 AC 26 ms
18,972 KB
testcase_18 AC 26 ms
18,976 KB
testcase_19 AC 20 ms
18,948 KB
testcase_20 AC 16 ms
16,888 KB
testcase_21 AC 12 ms
13,376 KB
testcase_22 AC 23 ms
18,992 KB
testcase_23 AC 27 ms
19,080 KB
testcase_24 AC 20 ms
19,080 KB
testcase_25 AC 13 ms
13,836 KB
testcase_26 AC 26 ms
19,080 KB
testcase_27 AC 27 ms
19,060 KB
testcase_28 AC 21 ms
19,116 KB
testcase_29 AC 18 ms
16,928 KB
testcase_30 AC 26 ms
18,976 KB
testcase_31 AC 19 ms
17,000 KB
testcase_32 AC 24 ms
19,080 KB
testcase_33 AC 25 ms
19,076 KB
testcase_34 AC 26 ms
18,992 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