結果

問題 No.886 Direct
ユーザー kotamanegikotamanegi
提出日時 2019-09-13 23:15:19
言語 C++17(clang)
(17.0.6 + boost 1.83.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 57 ms
26,880 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 99 ms
28,096 KB
testcase_03 AC 360 ms
28,408 KB
testcase_04 AC 55 ms
28,052 KB
testcase_05 AC 93 ms
28,720 KB
testcase_06 AC 92 ms
27,336 KB
testcase_07 AC 54 ms
27,040 KB
testcase_08 AC 91 ms
28,088 KB
testcase_09 AC 94 ms
27,348 KB
testcase_10 AC 100 ms
27,688 KB
testcase_11 AC 97 ms
27,124 KB
testcase_12 AC 94 ms
26,988 KB
testcase_13 AC 92 ms
27,392 KB
testcase_14 AC 94 ms
27,084 KB
testcase_15 AC 93 ms
27,348 KB
testcase_16 AC 95 ms
27,844 KB
testcase_17 AC 306 ms
28,336 KB
testcase_18 AC 337 ms
27,488 KB
testcase_19 AC 353 ms
27,840 KB
testcase_20 AC 315 ms
27,184 KB
testcase_21 AC 336 ms
28,308 KB
testcase_22 AC 342 ms
28,080 KB
testcase_23 AC 459 ms
53,340 KB
testcase_24 AC 472 ms
51,692 KB
testcase_25 AC 403 ms
45,988 KB
testcase_26 AC 425 ms
51,300 KB
testcase_27 AC 523 ms
69,136 KB
testcase_28 AC 571 ms
75,332 KB
testcase_29 AC 504 ms
81,336 KB
testcase_30 AC 431 ms
63,096 KB
testcase_31 AC 661 ms
81,640 KB
testcase_32 AC 637 ms
81,168 KB
testcase_33 AC 656 ms
81,248 KB
testcase_34 AC 648 ms
81,036 KB
testcase_35 AC 640 ms
80,708 KB
権限があれば一括ダウンロードができます

ソースコード

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