結果
問題 | No.1999 Lattice Teleportation |
ユーザー | 👑 potato167 |
提出日時 | 2022-06-14 23:18:18 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 145 ms / 2,000 ms |
コード長 | 2,902 bytes |
コンパイル時間 | 2,588 ms |
コンパイル使用メモリ | 214,760 KB |
実行使用メモリ | 16,136 KB |
最終ジャッジ日時 | 2024-11-08 11:09:03 |
合計ジャッジ時間 | 5,515 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 2 ms
5,248 KB |
testcase_03 | AC | 2 ms
5,248 KB |
testcase_04 | AC | 2 ms
5,248 KB |
testcase_05 | AC | 4 ms
5,248 KB |
testcase_06 | AC | 3 ms
5,248 KB |
testcase_07 | AC | 3 ms
5,248 KB |
testcase_08 | AC | 2 ms
5,248 KB |
testcase_09 | AC | 2 ms
5,248 KB |
testcase_10 | AC | 2 ms
5,248 KB |
testcase_11 | AC | 3 ms
5,248 KB |
testcase_12 | AC | 73 ms
13,956 KB |
testcase_13 | AC | 142 ms
16,004 KB |
testcase_14 | AC | 2 ms
5,248 KB |
testcase_15 | AC | 141 ms
15,876 KB |
testcase_16 | AC | 58 ms
10,552 KB |
testcase_17 | AC | 98 ms
15,080 KB |
testcase_18 | AC | 88 ms
14,164 KB |
testcase_19 | AC | 92 ms
14,620 KB |
testcase_20 | AC | 142 ms
15,880 KB |
testcase_21 | AC | 145 ms
16,004 KB |
testcase_22 | AC | 142 ms
16,136 KB |
testcase_23 | AC | 9 ms
5,248 KB |
testcase_24 | AC | 36 ms
8,044 KB |
testcase_25 | AC | 87 ms
14,152 KB |
testcase_26 | AC | 69 ms
11,744 KB |
testcase_27 | AC | 57 ms
10,424 KB |
testcase_28 | AC | 101 ms
15,264 KB |
testcase_29 | AC | 48 ms
9,372 KB |
testcase_30 | AC | 21 ms
5,780 KB |
testcase_31 | AC | 16 ms
5,248 KB |
testcase_32 | AC | 64 ms
11,156 KB |
ソースコード
#include <bits/stdc++.h> #pragma GCC optimize("Ofast") #define _GLIBCXX_DEBUG using namespace std; using ll=long long; using ld=long double; ll ILL=2167167167167167167; const int INF=2100000000; const ll mod=1e9+7; #define rep(i,a) for (ll i=0;i<a;i++) #define all(p) p.begin(),p.end() template<class T> using _pq = priority_queue<T, vector<T>, greater<T>>; template<class T> ll LB(vector<T> &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();} template<class T> ll UB(vector<T> &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();} template<class T> bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}else return 0;} template<class T> bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}else return 0;} template<class T> void So(vector<T> &v) {sort(v.begin(),v.end());} template<class T> void Sore(vector<T> &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});} void yneos(bool a){if(a) cout<<"Yes\n"; else cout<<"No\n";} template<class T> void vec_out(vector<T> &p){for(int i=0;i<(int)(p.size());i++){if(i) cout<<" ";cout<<p[i];}cout<<"\n";} bool _2d_comp(pair<ll,ll> &l,pair<ll,ll> &r){ if(l.first==0&&r.first==0) return l.second>r.second; if(l.first==0){ return l.second>0||r.first<0; } if(r.first==0){ return !(r.second>0||l.first<0); } if((l.first<0)^(r.first<0)) return l.first>0; return l.first*r.second<l.second*r.first; } bool _2d_same(pair<ll,ll> &l,pair<ll,ll> &r){ if((l.first<0)^(r.first<0)) return false; if((l.second<0)^(r.first<0)) return false; if((l.first==0)&&(r.second==0)) return false; if((l.second==0)&&(r.first==0)) return false; return l.first*r.second==l.second*r.first; } void solve(); // KITAKU BEATS int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int t=1; //cin>>t; rep(i,t) solve(); } void solve(){ int N; cin>>N; vector<pair<pair<ll,ll>,pair<ll,ll>>> p; vector<ll> X={0},Y={0}; vector<ll> A(N),B(N); rep(i,N){ ll a,b; cin>>a>>b; A[i]=a,B[i]=b; if(a==0&&b==0) continue; pair<ll,ll> S={-b,a},T={b,-a}; p.push_back({S,{1,i}}); p.push_back({T,{-1,i}}); if(!_2d_comp(S,T)){ X[0]+=a,Y[0]+=b; } } if(p.empty()){ cout<<"1\n"; return ; } sort(all(p),[&](auto l,auto r){ return _2d_comp(l.first,r.first); }); int L=p.size(),p_ind=0,a_ind=0; while(L!=p_ind){ X.push_back(X[a_ind]); Y.push_back(Y[a_ind]); a_ind++; do{ ll s=p[p_ind].second.first; X[a_ind]+=A[p[p_ind].second.second]*s; Y[a_ind]+=B[p[p_ind].second.second]*s; p_ind++; }while(p_ind!=L&&_2d_same(p[p_ind].first,p[p_ind-1].first)); } if(a_ind==2){ cout<<(1+__gcd(abs(X[0]-X[1]),abs(Y[0]-Y[1])))%mod<<"\n"; return ; } //vec_out(X); //vec_out(Y); ll ans=0; rep(i,a_ind){ ans+=__gcd(abs(X[i]-X[i+1]),abs(Y[i]-Y[i+1])); ans%=mod; ans+=((X[i+1]%mod)*(Y[i]%mod))%mod; ans-=((X[i]%mod)*(Y[i+1]%mod))%mod; ans%=mod; } ans+=2; ans=(ans%mod+mod)%mod; cout<<(ans*((mod+1)/2))%mod<<"\n"; }