結果

問題 No.890 移調の限られた旋法
コンテスト
ユーザー polylogK
提出日時 2019-09-20 22:14:50
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,249 bytes
コンパイル時間 3,062 ms
コンパイル使用メモリ 170,184 KB
実行使用メモリ 19,108 KB
最終ジャッジ日時 2024-09-14 17:59:42
合計ジャッジ時間 3,171 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15 WA * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std::literals::string_literals;
using i64 = long long;
using std::cout;
using std::endl;
using std::cin;

template<typename T>
std::vector<T> make_v(size_t a){return std::vector<T>(a);}

template<typename T,typename... Ts>
auto make_v(size_t a,Ts... ts){
  return std::vector<decltype(make_v<T>(ts...))>(a,make_v<T>(ts...));
}

int main() {
	int n, k; scanf("%d%d", &n, &k);

	const int MOD = 1e9 + 7;
	auto pow = [&](i64 x, i64 n) {
		i64 ret = 1;
		while(n) {
			if(n & 1) (ret *= x) %= MOD;
			(x *= x) %= MOD;
			n >>= 1;
		}
		return ret;
	};
	auto inv = [&](i64 n) {
		return pow(n, MOD - 2);
	};
	std::vector<i64> fact(n + 1, 1);
	for(i64 i = 0; i < n; i++) fact[i + 1] = fact[i] * (i + 1);
	auto comb = [&](int n, int r) {
		if(n < r) return 0LL;
	
		return fact[n] * inv(fact[r]) % MOD * inv(fact[n - r]) % MOD;
	};
	
	i64 ans = 0;
	std::vector<i64> latte(n, 0);
	for(int i = 1; i < n; i++) {
		if(n % i) continue;
	
		int cnt = n / i;
		if(k % cnt) continue;
		int loop = k / cnt;

		(latte[i] = comb(i, loop)) %= MOD;
		
		for(int j = 1; j < i; j++) {
			if(i % j) continue;

			latte[i] -= latte[j];
		}
	}

	for(auto v: latte) (ans += v % MOD + MOD) %= MOD;
	printf("%lld\n", ans);
	return 0;
}
0