結果

問題 No.890 移調の限られた旋法
ユーザー polylogKpolylogK
提出日時 2019-09-20 22:09:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,360 bytes
コンパイル時間 1,933 ms
コンパイル使用メモリ 168,620 KB
実行使用メモリ 18,944 KB
最終ジャッジ日時 2024-09-14 17:51:35
合計ジャッジ時間 2,885 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14 WA * 18
権限があれば一括ダウンロードができます

ソースコード

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);
	if((n ^ k) & 1) {
		printf("0\n");
		return 0;
	}

	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;

		(ans += comb(i, loop)) %= MOD;
		(latte[i] = comb(i, loop)) %= MOD;
	
		for(int j = 2; j * j <= i; j++) {
			if(i % j) continue;
			
			ans = (ans - latte[j] + MOD) % MOD;
			if(j != i / j) ans = (ans - latte[i / j] + MOD) % MOD;
		}
	}
	printf("%lld\n", ans);
	return 0;
}
0