結果

問題 No.644 G L C C D M
コンテスト
ユーザー e869120
提出日時 2018-02-02 21:55:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 950 bytes
コンパイル時間 889 ms
コンパイル使用メモリ 82,904 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-12-31 07:31:09
合計ジャッジ時間 2,098 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23 WA * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <cmath>
#include <queue>
#include <functional>
#include <map>
using namespace std;

long long n, m, mod = 1000000007; bool primes[1000009];

long long prime(long long p) {
	if (p == 0) return 0;
	if (p <= 2500) {
		long long R = 0;
		for (int i = 1; i <= p; i++) {
			for (int j = 1; j <= p; j++) {
				if (i == j) continue;
				if (__gcd(i, j) == 1) R++;
			}
		}
		return R;
	}
	long long R = p*p;
	for (int i = 2; i <= p; i++) {
		if (primes[i] == true)R -= 1LL*(p / i)*(p / i);
	}
	return R - 1;
}

int main() {
	for (int i = 2; i <= 1000008; i++) primes[i] = true;
	for (int i = 2; i*i <= 1000008; i++) {
		for (int j = i*i; j <= 1000008; j += i) primes[j] = false;
	}
	cin >> n >> m;
	long long G = prime(n / m);
	long long ret = 1;
	for (int i = 1; i <= n - 2; i++) { ret *= i; ret %= mod; }
	ret *= (G % mod); ret %= mod;
	cout << ret << endl;
	return 0;
}
0