結果

問題 No.644 G L C C D M
ユーザー ianCK
提出日時 2020-08-12 16:16:01
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 6 ms / 2,000 ms
コード長 828 bytes
コンパイル時間 279 ms
コンパイル使用メモリ 37,888 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-09 18:27:15
合計ジャッジ時間 1,295 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 27
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:35:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   35 |         scanf("%d%d", &n, &m);
      |         ~~~~~^~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <stdio.h>
#include <vector>
using namespace std;
#define PB push_back

typedef long long int ll;
constexpr int kMod = int(1E9 + 7), kN = int(2E5 + 10);

ll mu[kN];
vector<int> prime;

void pre() {
	mu[1] = 1;
	for (int i = 2; i < kN; i++) mu[i] = -2;
	for (int i = 2; i < kN; i++) {
		if (mu[i] == -2) {
			mu[i] = kMod - 1;
			prime.PB(i);
		}
		for (int j : prime) {
			if (i * j >= kN) break;
			mu[i * j] = kMod - mu[i];
			if (i % j == 0) {
				mu[i * j] = 0;
				break;
			}
		}
	}
	return ;
}

int main() {
	int n, m;
	ll ans = 0;
	scanf("%d%d", &n, &m);
	pre();
	if (n < m) printf("0\n");
	else {
		for (int k = 1; k <= n / m; k++) 
			ans += mu[k] * (n / k / m) * (n / k / m) % kMod;
		ans = (ans + kMod - 1) % kMod;
		for (int i = 1; i <= n - 2; i++) ans = ans * i % kMod;
		printf("%lld\n", ans% kMod);
	}
}
0