結果

問題 No.1627 三角形の成立
コンテスト
ユーザー kwm_t
提出日時 2021-08-29 01:57:11
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,480 bytes
コンパイル時間 4,214 ms
コンパイル使用メモリ 232,060 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2024-11-22 01:32:07
合計ジャッジ時間 5,063 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 1
other AC * 3 WA * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
#include "atcoder/all"
using namespace std;
using namespace atcoder;
using mint = modint1000000007;
const int mod = 1000000007;
//using mint = modint998244353;
//const int mod = 998244353;
//const int INF = 1e9;
//const long long LINF = 1e18;
//const bool debug = false;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
// f(k) = sum_{GCD(i, j) = k} g(i)h(j)
// F(k) = sum_{k | i} f(i)
template<class T> struct FastGCDConvolution {
	int n;
	FastGCDConvolution(int n) : n(n) {}

	void zeta(vector<mint> &v) {
		int n = v.size();
		for (int k = 1; k < n; ++k) {
			for (int i = k * 2; i < n; i += k) {
				v[k] += v[i];
			}
		}
	}

	void mebius(vector<T> &v) {
		int n = v.size();
		for (int k = n - 1; k >= 1; --k) {
			for (int i = k * 2; i < n; i += k) {
				v[k] -= v[i];
			}
		}
	}

	vector<T> mult(vector<T> &a, vector<T> &b) {
		vector<T> c(n, 0);
		vector<T> subc(n, 0);
		rep2(i, 1, a.size())subc[i] += a[i] * b[0];
		rep2(i, 1, b.size())subc[i] += a[0] * b[i];
		zeta(a); zeta(b);
		for (int i = 0; i < min(a.size(), b.size()); ++i) c[i] = a[i] * b[i];
		//rep(i, c.size()) cout << c[i].val() << " ";
		//cout << endl;
		mebius(c);
		rep(i, n)c[i] += subc[i];
		//rep(i, c.size()) cout << c[i].val() << " ";
		//cout << endl;
		return c;
	}
};
long long gcd(long long a, long long b) {
	if (b == 0) return a;
	else return gcd(b, a%b);
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int n, m; cin >> n >> m;
	mint nm = ((mint)n) * m;
	mint ans = nm * (nm - 1)*(nm - 2) / 6;
	vector<mint>v0(n), v1(m);

	rep(i, n) v0[i] = n - i;
	rep(i, m) v1[i] = m - i;
	/*vector<mint>vv(max(n, m), 0);
	rep(i, n)rep(j, m) {
		vv[gcd(i, j)] += v0[i] * v1[j];
	}
	rep(i, vv.size()) cout << vv[i].val() << " ";
	cout << endl;
	rep2(i, 1, vv.size()) {
		ans -= vv[i] * (i - 1);
	}*/
	FastGCDConvolution<mint> fz(max(n, m));
	auto v = fz.mult(v0, v1);
	rep2(i, 1, max(n, m)) ans -= v[i] * (i - 1);
	cout << ans.val() << endl;
	return 0;
}
0