結果

問題 No.644 G L C C D M
コンテスト
ユーザー kurenai3110
提出日時 2018-02-02 23:36:39
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 847 bytes
コンパイル時間 606 ms
コンパイル使用メモリ 57,444 KB
最終ジャッジ日時 2024-11-14 20:19:50
合計ジャッジ時間 1,197 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
main.cpp: In function ‘int solve(int, int)’:
main.cpp:29:63: error: ‘pow’ was not declared in this scope
   29 |                 for (int i = 0; i < prime.size(); i++)tmp *= (pow(prime[i].first, prime[i].second) - pow(prime[i].first, prime[i].second - 1));
      |                                                               ^~~

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <string>
using namespace std;

#define MOD (int)(1e9+7)

int solve(int N,int M) {
	long long cnt = N / M;

	long long ans = 0;
	for (int i = 2; i <= cnt; i++) {
		int x = i;
		vector<pair<int, int>>prime;
		for (int j = 2; j*j <= x; j++) {
			if (x%j == 0) {
				int tmp = 0;
				while (x%j == 0) {
					x /= j;
					tmp++;
				}
				prime.push_back(make_pair(j, tmp));
			}
		}
		if (x != 1)prime.push_back(make_pair(x, 1));

		long long tmp = 1;
		for (int i = 0; i < prime.size(); i++)tmp *= (pow(prime[i].first, prime[i].second) - pow(prime[i].first, prime[i].second - 1));
		ans += 2*tmp;
	}
	ans %= MOD;
	for (int i = N - 2; i >= 1; i--) {
		ans *= i;
		ans %= MOD;
	}
	return ans;
}

int main()
{
	int N, M; cin >> N >> M;

	cout << solve(N, M) << endl;

    return 0;
}
0