結果
問題 |
No.1999 Lattice Teleportation
|
ユーザー |
![]() |
提出日時 | 2022-03-21 21:40:54 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 176 ms / 2,000 ms |
コード長 | 1,378 bytes |
コンパイル時間 | 1,019 ms |
コンパイル使用メモリ | 85,784 KB |
最終ジャッジ日時 | 2025-01-28 10:56:12 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 29 |
ソースコード
#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 gcd(ll a, ll b) { a = abs(a); b = abs(b); while (b != 0) { ll r = a % b; a = b; b = r; } return a; } ll get_mod(ll n) { if (n < 0) { return mod - (-n) % mod; } return n % mod; } ll area(V<ll>& v1, V<ll>& v2) { ll x1 = get_mod(v1[0]), y1 = get_mod(v1[1]); ll x2 = get_mod(v2[0]), y2 = get_mod(v2[1]); return get_mod(x1 * y2 % mod - x2 * y1 % mod); } 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 tmp, n = 0; cin >> tmp; VV<ll> v(tmp, V<ll>(2)); rep(i, 0, tmp) { ll x, y; cin >> x >> y; if (x == 0 and y == 0) { continue; } if (y < 0) { x = -x; y = -y; } v[n] = { x,y }; n++; } v.resize(n); sort(all(v), arg_sort); ll S = 0, bh = 0; ll lx = 0, ly = 0; for (auto& l : v) { ll x = l[0], y = l[1]; V<ll> v1 = { lx,ly }; V<ll> v2 = { lx + x,ly + y }; S += area(v1, v2); S %= mod; bh += gcd(x, y); bh %= mod; lx = v2[0]; ly = v2[1]; } ll ans = (S + bh + 1) % mod; cout << ans << endl; return 0; }