結果

問題 No.1036 Make One With GCD 2
ユーザー Mi
提出日時 2020-04-24 22:39:03
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 951 ms / 2,000 ms
コード長 2,855 bytes
コンパイル時間 1,843 ms
コンパイル使用メモリ 174,476 KB
実行使用メモリ 15,372 KB
最終ジャッジ日時 2024-09-16 13:24:13
合計ジャッジ時間 17,615 ms
ジャッジサーバーID
(参考情報)
judge2 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll,ll> P;
typedef pair<int,int> Pi;
#define rep(i,n) for(ll i=0;i<n;i++)
#define FOR(i,a,b) for(ll i=a;i<b;i++)
#define fi first
#define se second
#define endl "\n"
template<typename T> inline bool chmax(T &a, T b){if(a<b){a=b;return true;}return false;}
template<typename T> inline bool chmin(T &a, T b){if(a>b){a=b;return true;}return false;}
template<typename T> ostream& operator<<(ostream& s,const complex<T>& d) {return s<<"("<<d.real()<<", "<<d.imag()<< ")";}
template<typename T1, typename T2> ostream& operator<<(ostream& s,const pair<T1,T2>& d) {return s<<"("<<d.first<<", "<<d.second<<")";}
template<typename T> ostream& operator<<(ostream& s, const vector<T>& d){int len=d.size();rep(i,len){s<<d[i];if(i<len-1) s<<" ";}return s;}
template<typename T> ostream& operator<<(ostream& s,const vector<vector<T>>& d){int len=d.size();rep(i,len){s<<d[i]<<endl;}return s;}
template<typename T> ostream& operator<<(ostream& s,const set<T>& v){s<<"{ ";for(auto itr=v.begin();itr!=v.end();++itr) {if (itr!=v.begin()) {s<< ", "
    ;}s<<(*itr);}s<<" }";return s;}
template<typename T1, typename T2> ostream& operator<<(ostream& s,const map<T1,T2>& m){s<<"{"<<endl;for(auto itr=m.begin();itr!=m.end();++itr){s<<" "
    <<(*itr).first<<" : "<<(*itr).second<<endl;}s<<"}"<<endl;return s;}
const ll mod=1e9+7;
const ll inf=1e17;
const int INF=1e9;
const double PI=acos(-1);
const double EPS=1e-10;
//a,b
template <class T>
T gcd(T a,T b) {
if (b==0) return a;
else return gcd(b,a%b);
}
template<typename Monoid,typename F>
struct SegmentTree{
int sz;
vector<Monoid> node;
const F op;//
const Monoid e;//
SegmentTree(int n,const F op,const Monoid &e):op(op),e(e){
sz=1;
while(sz<n) sz*=2;
node.assign(2*sz,e);
}
void build(const vector<Monoid> &v){
int n=v.size();
for(int i=0;i<n;i++) node[i+sz-1]=v[i];
for(int i=sz-2;i>=0;i--) node[i]=op(node[2*i+1],node[2*i+2]);
}
void update(int k,const Monoid &x){
k+=sz-1;
node[k]=x;
while(k>0){
k=(k-1)/2;
node[k]=op(node[2*k+1],node[2*k+2]);
}
}
Monoid query(int a,int b,int k=0,int l=0,int r=-1){
if(r<0) r=sz;
if(r<=a || b<=l) return e;
if(a<=l && r<=b) return node[k];
else{
Monoid vl=query(a,b,2*k+1,l,(l+r)/2);
Monoid vr=query(a,b,2*k+2,(l+r)/2,r);
return op(vl,vr);
}
}
Monoid operator[](const int &k)const{
return node[k+sz-1];
}
};
int main(){
cin.tie(0);ios::sync_with_stdio(false);
ll n;
cin>>n;
vector<ll> a(n);
rep(i,n){
cin>>a[i];
}
auto op=[](ll a,ll b){return gcd(a,b);};
SegmentTree<ll,decltype(op)> seg(n,op,0);
seg.build(a);
ll ans=0,r=0,now=0;
rep(l,n){
while(r<n && now!=1){
now=gcd(now,a[r]);
r++;
}
if(now==1){
ans+=n-r+1;
}
now=seg.query(l+1,r);
}
cout<<ans<<endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0