結果

問題 No.1999 Lattice Teleportation
ユーザー kwm_tkwm_t
提出日時 2022-07-02 11:00:07
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 109 ms / 2,000 ms
コード長 1,961 bytes
コンパイル時間 5,279 ms
コンパイル使用メモリ 268,412 KB
実行使用メモリ 7,920 KB
最終ジャッジ日時 2023-08-18 02:04:42
合計ジャッジ時間 7,909 ms
ジャッジサーバーID
(参考情報)
judge11 / judge14
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 3 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 1 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 2 ms
4,380 KB
testcase_12 AC 53 ms
7,768 KB
testcase_13 AC 106 ms
7,920 KB
testcase_14 AC 2 ms
4,376 KB
testcase_15 AC 107 ms
7,892 KB
testcase_16 AC 39 ms
5,700 KB
testcase_17 AC 65 ms
7,588 KB
testcase_18 AC 57 ms
7,064 KB
testcase_19 AC 60 ms
7,180 KB
testcase_20 AC 109 ms
7,768 KB
testcase_21 AC 109 ms
7,728 KB
testcase_22 AC 108 ms
7,764 KB
testcase_23 AC 7 ms
4,376 KB
testcase_24 AC 24 ms
4,708 KB
testcase_25 AC 56 ms
6,996 KB
testcase_26 AC 47 ms
6,320 KB
testcase_27 AC 39 ms
5,628 KB
testcase_28 AC 65 ms
7,548 KB
testcase_29 AC 32 ms
5,224 KB
testcase_30 AC 14 ms
4,376 KB
testcase_31 AC 12 ms
4,376 KB
testcase_32 AC 43 ms
6,136 KB
権限があれば一括ダウンロードができます

ソースコード

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;
#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<long long,long long>
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; }
bool comp(const P& lh, const P& rh) {
	return lh.second * rh.first < rh.second * lh.first;
}
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);
	vector<P>v;
	int n; cin >> n;
	int c = 0;
	rep(i, n) {
		long long a, b; cin >> a >> b;
		if ((0 == a) && (0 == b))continue;
		if (a < 0) {
			a = -a;
			b = -b;
		}
		v.emplace_back(a, b);
	}
	n = v.size();
	if (0 == n) {
		cout << 1 << endl;
		return 0;
	}
	sort(all(v), comp);
	v.resize(n * 2);
	rep(i, n) v[i + n] = { -v[i].first, -v[i].second };
	vector<pair<mint,mint>>point(v.size());
	point[0] = { 0,0 };
	rep(i, v.size() - 1) {
		point[i + 1] = { point[i].first + v[i].first, point[i].second + v[i].second };
	}
	mint s2 = 0;
	rep(i, point.size()) {
		int j = (i + 1) % point.size();
		s2 += point[j].second * point[i].first - point[j].first * point[i].second;
	}
	// s2 = 辺上 + 内部 * 2 - 2
	mint on = 0, in = 0;
	rep(i, v.size()) on += abs(gcd(v[i].first, v[i].second));
	in = (s2 + 2 - on) / 2;
	mint ans = in + on;
	cout << ans.val() << endl;
	return 0;
}
0