結果

問題 No.890 移調の限られた旋法
コンテスト
ユーザー Y17
提出日時 2019-09-20 22:49:20
言語 C++11（廃止可能性あり）
(gcc 13.3.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

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