結果

問題 No.1036 Make One With GCD 2
ユーザー mugen_1337mugen_1337
提出日時 2020-04-25 09:52:16
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 265 ms / 2,000 ms
コード長 2,121 bytes
コンパイル時間 2,247 ms
コンパイル使用メモリ 202,920 KB
最終ジャッジ日時 2025-01-10 00:56:42
ジャッジサーバーID
（参考情報） 		judge1 / judge6
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;}

// Semigroup Sliding-Window-AGregation
template<typename Semigroup>
struct SWAG{
    using F=function<Semigroup(Semigroup,Semigroup)>;
    const F f;
    stack<pair<Semigroup,Semigroup>> fr,ba;
    SWAG(const F f):f(f){}
    Semigroup get(){
        if(fr.empty()) return ba.top().second;
        if(ba.empty()) return fr.top().second;
        return f(fr.top().second,ba.top().second);
    }
    void push(Semigroup 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<Semigroup,Semigroup> 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];

    auto f=[](ll a,ll b){return gcd(a,b);};
    SWAG<ll> que(f);
    ll res=0;
    for(int l=0,r=0;l<n;l++){
        while(r<n){
            if(que.empty() and a[r]==1) break;
            if(!que.empty() and f(que.get(),a[r])==1) break;
            que.push(a[r]);r++;
        }
        res+=(ll)que.size();
        if(l==r) r++;
        else que.pop();
    }
    cout<<(ll)n*(n+1)/2-res<<endl;
    return 0;
}
0