結果

問題 No.886 Direct
コンテスト
ユーザー kotamanegi
提出日時 2019-09-13 23:15:19
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
AC  
実行時間 661 ms / 4,000 ms
コード長 2,667 bytes
コンパイル時間 1,355 ms
コンパイル使用メモリ 151,832 KB
実行使用メモリ 81,640 KB
最終ジャッジ日時 2024-11-30 14:57:02
合計ジャッジ時間 13,049 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define  _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <algorithm>
#include <utility>
#include <functional>
#include <cstring>
#include <queue>
#include <stack>
#include <math.h>
#include <iterator>
#include <vector>
#include <string>
#include <set>
#include <math.h>
#include <iostream>
#include <random>
#include<map>
#include <iomanip>
#include <time.h>
#include <stdlib.h>
#include <list>
#include <typeinfo>
#include <list>
#include <set>
#include <cassert>
#include<fstream>
#include <unordered_map>
#include <cstdlib>
#include <complex>
#include <cctype>
#include <bitset>
using namespace std;
typedef string::const_iterator State;
#define Ma_PI 3.141592653589793
#define eps 1e-5
#define LONG_INF 1e18
#define GOLD 1.61803398874989484820458
#define MAX_MOD 1000000007LL
#define GYAKU 500000004LL
#define MOD 998244353LL
#define seg_size 262144*4
#define REP(a,b) for(long long a = 0;a < b;++a)
int gcd(int a, int b) {
	if (b == 0) return a;
	return gcd(b, a % b);
}
long long cnter[4000000] = {};
long long wow[4000000] = {};
long long calc(long long start, long long hi) {
	long long kousuu = (start+hi-1) / hi;
	long long ende = start - (kousuu-1LL) * hi;
	long long hoge = ende + start;
	hoge %= MAX_MOD;
	hoge *= kousuu;
	hoge %= MAX_MOD;
	hoge *= GYAKU;
	hoge %= MAX_MOD;
	return hoge;
}
int main(){
	iostream::sync_with_stdio(false);
	cin.tie(0);
	long long h, w;
	cin >> h >> w;
	long long ans = (h - 1) * w + (w - 1) * h;
	ans %= MAX_MOD;
	vector<long long> sosuu;
	vector<pair<long long, long long>> gogo;
	for (long long i = 2; i < max(h,w)+1; ++i) {
		if (wow[i] == 0) {
			sosuu.push_back(i);
			for (int q = 1; i * q <= 3000000; ++q) {
				if (gcd(i, q) == 1) {
					wow[i * q]++;
				}
				else {
					wow[i * q] += 100000;
				}
			}
		}
	}
	for (int i = 2; i < max(h,w)+1; ++i) {
		if (wow[i] < 100000) {
			gogo.push_back(make_pair(i, wow[i]));
		}
	}
	for (long long i = 1; i < h; ++i) {
		cnter[i] = (w - 1LL) * w / 2LL;
		cnter[i] %= MAX_MOD;
	}
	for (int i = 0; i < gogo.size(); ++i) {
		if (gogo[i].first > h||gogo[i].first > w) break;
		long long downing = calc(w - gogo[i].first, gogo[i].first);
		if (gogo[i].second % 2 == 1) {
			downing = MAX_MOD - downing;
		}
		for (int t = 1; gogo[i].first * t < h; ++t) {
			cnter[gogo[i].first * t] += downing;
			cnter[gogo[i].first * t] %= MAX_MOD;
		}
	}
	long long tmp = 0;
	for (long long i = 1; i < h; ++i) {
		tmp += cnter[i] * 2LL * (h - i);
		tmp %= MAX_MOD;
	}
	ans += tmp;
	ans %= MAX_MOD;
	cout << ans<< endl;
	/*
	for (int i = 1; i < h; ++i) {
		for (int q = 1; q < w; ++q) {
			if (gcd(i, q) == 1) {
				ans += 2 * (h - i) * (w - q);
			}
		}
	}
	cout << ans << endl;
	*/
}
0