結果

問題 No.1999 Lattice Teleportation
コンテスト
ユーザー MasKoaTS
提出日時 2022-02-22 23:15:30
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 1,779 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,028 ms
コンパイル使用メモリ 84,372 KB
最終ジャッジ日時 2025-01-28 01:28:44
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 27 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
#include <algorithm>
#define rep(i, l, n) for (int i = (l); i < (n); i++)
#define all(x) x.begin(), x.end()
using namespace std;
using ll = long long;
template <class T>
using V = vector<T>;
template <class T>
using VV = V<V<T> >;
const ll mod = 1000000007;


ll floor_sum(ll n, ll m, ll a, ll b) {
	if (m == 0) {
		return 0;
	}
	ll res = 0;
	while (true) {
		if (a >= m or a < 0) {
			res += n * (n - 1) * (a / m) / 2;
			a %= m;
		}
		if (b >= m or b < 0) {
			res += n * (b / m);
			b %= m;
		}
		ll y_max = a * n + b;
		if (y_max < m) {
			break;
		}
		n = y_max / m;	b = y_max % m;
		m = a;	a = m;
	}
	return res;
}

ll gcd(ll a, ll b) {
	a = abs(a);	b = abs(b);
	while (b != 0) {
		ll r = a % b;
		a = b;
		b = r;
	}
	return a;
}

bool arg_sort(const V<ll>& v1, const V<ll>& v2) {
	ll x1 = v1[0], y1 = v1[1];
	ll x2 = v2[0], y2 = v2[1];
	return (x1 * y2 > x2 * y1);
}


int main(void) {
	int n;	cin >> n;
	VV<ll> v(n, V<ll>(2));
	rep(i, 0, n) {
		ll x, y;	cin >> x >> y;
		if (y < 0) {
			x = -x;	y = -y;
		}
		v[i] = { x,y };
	}
	sort(all(v), arg_sort);

	ll ans = 1;
	for(auto& p : v) {
		ans += gcd(p[0], p[1]);
		ans %= mod;
	}

	ll last = 0;
	for (auto itr = rbegin(v); itr != rend(v); ++itr) {
		ll x = (*itr)[0], y = (*itr)[1];
		ll a = 1;
		ll s = 0;
		if (x < 0) {
			a = -1;
			x = -x;
		}
		s += floor_sum(x, x, y, 0) % mod + x * last % mod;
		s %= mod;
		ans += mod + a * s;
		ans %= mod;
		last += mod + y;
		last %= mod;
	}

	last = 0;
	for (auto& p : v) {
		ll x = p[0], y = p[1];
		ll a = -1;
		ll s = 0;
		if (x < 0) {
			a = 1;
			x = -x;
		}
		s += floor_sum(x, x, y, 0) % mod + x * last % mod;
		s %= mod;
		ans += mod + a * s;
		ans %= mod;
		last += mod + y;
		last %= mod;
	}

	cout << ans << endl;
	return 0;
}
0