結果
問題 |
No.1036 Make One With GCD 2
|
ユーザー |
![]() |
提出日時 | 2020-04-25 20:17:37 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 326 ms / 2,000 ms |
コード長 | 2,105 bytes |
コンパイル時間 | 2,220 ms |
コンパイル使用メモリ | 202,488 KB |
最終ジャッジ日時 | 2025-01-10 01:28:52 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define ALL(x) x.begin(),x.end() #define rep(i,n) for(int i=0;i<(n);i++) #define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl; #define INF 1000000000 #define mod 1000000007 using ll=long long; const ll LINF=1001002003004005006ll; int dx[]={1,0,-1,0}; int dy[]={0,1,0,-1}; // ll gcd(ll a,ll b){return b?gcd(b,a%b):a;} template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;} template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;} // Monoid Sliding-Window-AGregation template<typename Monoid> struct SWAG{ using F=function<Monoid(Monoid,Monoid)>; const F f; const Monoid gen; stack<pair<Monoid,Monoid>> fr,ba; SWAG(const F f,const Monoid &gen):f(f),gen(gen){} Monoid get(){ ll a=fr.empty()?gen:fr.top().second; ll b=ba.empty()?gen:ba.top().second; return f(a,b); } void push(Monoid x){ if(ba.empty()) ba.push(make_pair(x,x)); else ba.push(make_pair(x,f(ba.top().second,x))); } void pop(){ if(empty()) return ; if(!fr.empty()) fr.pop(); else{ while(!ba.empty()){ pair<Monoid,Monoid> p=ba.top();ba.pop(); if(fr.empty()) fr.push(make_pair(p.first,p.first)); else fr.push(make_pair(p.first,f(p.first,fr.top().second))); } fr.pop(); } } size_t size(){ return fr.size()+ba.size(); } bool empty(){ return size()==0; } }; signed main(){ cin.tie(0); ios::sync_with_stdio(0); int n;cin>>n; vector<ll> a(n); rep(i,n){ cin>>a[i]; } function<ll(ll,ll)> g=[](ll a,ll b){return gcd(a,b);}; SWAG<ll> que(g,0); ll r=0,res=0; for(ll l=0;l<n;l++){ while(r<n and g(que.get(),a[r])!=1){ que.push(a[r]); r++; } res+=(ll)que.size(); if(r==l){ r++; }else{ que.pop(); } } cout<<(ll)n*(n+1)/2-res<<endl; return 0; }