結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
|
提出日時 | 2020-04-25 13:43:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,191 ms / 2,000 ms |
コード長 | 2,070 bytes |
コンパイル時間 | 1,505 ms |
コンパイル使用メモリ | 170,868 KB |
実行使用メモリ | 15,360 KB |
最終ジャッジ日時 | 2024-09-16 13:42:21 |
合計ジャッジ時間 | 17,674 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef long double ld;template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }const ll INF=1LL<<60;const int inf=1<<30;const int mod=1e9+7;const int MOD=998244353;ll gcd(ll a,ll b){return (b?gcd(b,a%b):a);}template <typename T, typename E, typename F, typename G>struct SegmentTree{int n;F f;G g;T ti;vector<T> dat;SegmentTree(){};SegmentTree(F f,G g,T ti):f(f),g(g),ti(ti){}void init(int n_){n=1;while(n<n_) n<<=1;dat.assign(n<<1,ti);}void build(const vector<T> &v){int n_=v.size();init(n_);for(int i=0;i<n_;i++) dat[n+i]=v[i];for(int i=n-1;i;i--)dat[i]=f(dat[(i<<1)|0],dat[(i<<1)|1]);}void update(int k,const E &x){k += n;dat[k] = g(dat[k], x);while(k>>=1)dat[k]=f(dat[(k<<1)|0],dat[(k<<1)|1]);}T operator [](int k) const { return dat[k+n]; }// [a,b)T query(int a,int b) const {T vl=ti,vr=ti;for(int l=a+n,r=b+n;l<r;l>>=1,r>>=1) {if(l&1) vl=f(vl,dat[l++]);if(r&1) vr=f(dat[--r],vr);}return f(vl,vr);}};int main(){cin.tie(0);ios::sync_with_stdio(false);ll n;cin >> n;using T = ll; // type Tusing E = ll; // type Eauto f = [](T a, T b){ // return type T valuereturn gcd(a,b);};auto g = [](T a, E b){ // return type T valuereturn b; // return b;};T ti = 0; // identify elementSegmentTree<T,E,decltype(f),decltype(g)> seg(f,g,ti);vector<ll> a(n);seg.build(a);for(int i=0;i<n;i++){ll t;cin >> t;seg.update(i,t);}int r=0;ll ans=0;for(int i=0;i<n;i++){while(r<n&&seg.query(i,r+1)!=1){r++;}ans+=n-r;}cout << ans << endl;}