結果

問題 No.1036 Make One With GCD 2
ユーザー IKyopro
提出日時 2020-04-24 21:58:47
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,527 ms / 2,000 ms
コード長 2,947 bytes
コンパイル時間 2,973 ms
コンパイル使用メモリ 198,048 KB
最終ジャッジ日時 2025-01-09 23:43:27
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

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

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
template <class T, class U> using Pa = pair<T, U>;
template <class T> using vec = vector<T>;
template <class T> using vvec = vector<vec<T>>;
template<typename Monoid,typename F>
class SegmentTree{
private:
int sz;
vector<Monoid> seg;
const F op;
const Monoid e;
public:
SegmentTree(int n,const F op,const Monoid &e):op(op),e(e){
sz = 1;
while(sz<=n) sz <<= 1;
seg.assign(2*sz,e);
}
void set(int k, const Monoid &x){
seg[k+sz] = x;
}
void build(){
for(int i=sz-1;i>0;i--){
seg[i] = op(seg[2*i],seg[2*i+1]);
}
}
void update(int k,const Monoid &x){
k += sz;
seg[k] = x;
while(k>>=1){
seg[k] = op(seg[2*k],seg[2*k+1]);
}
}
Monoid query(int l,int r){
Monoid L = e,R = e;
for(l+=sz,r+=sz;l<r;l>>=1,r>>=1){
if(l&1) L = op(L,seg[l++]);
if(r&1) R = op(seg[--r],R);
}
return op(L,R);
}
Monoid operator[](const int &k)const{
return seg[k+sz];
}
template<typename C>
int find_subtree(int a,const C &check,Monoid &M,bool type) {
while(a<sz) {
Monoid nxt = type ? op(seg[2*a+type],M) : op(M,seg[2*a+type]);
if(check(nxt)) a = 2*a+type;
else M = nxt, a = 2*a+1-type;
}
return a - sz;
}
template<typename C>
int find_first(int a,const C &check) {
Monoid L = e;
if(a <= 0) {
if(check(op(L,seg[1]))) return find_subtree(1,check,L,false);
return -1;
}
int b = sz;
for(a+=sz,b+=sz;a<b;a>>=1,b>>=1) {
if(a&1) {
Monoid nxt = op(L,seg[a]);
if(check(nxt)) return find_subtree(a,check,L,false);
L = nxt;
++a;
}
}
return -1;
}
template<typename C>
int find_last(int b,const C &check) {
Monoid R = e;
if(b >= sz) {
if(check(op(seg[1], R))) return find_subtree(1, check, R, true);
return -1;
}
int a = sz;
for(b+=sz; a<b; a>>=1,b>>=1) {
if(b&1) {
Monoid nxt = op(seg[--b],R);
if(check(nxt)) return find_subtree(b,check,R,true);
R = nxt;
}
}
return -1;
}
};
int main(){
cin.tie(0);
ios::sync_with_stdio(false);
int N;
cin >> N;
vec<ll> A(N);
for(int i=0;i<N;i++) cin >> A[i];
auto op = [&](ll a,ll b){return gcd(a,b);};
SegmentTree<ll,decltype(op)> seg(N,op,0);
for(int i=0;i<N;i++) seg.set(i,A[i]);
seg.build();
ll ans = 0;
for(int i=0;i<N;i++){
if(seg.query(i,N)>1) continue;
int r = seg.find_first(i,[&](ll b){return gcd(A[i],b)==1;});
ans += N-r;
}
cout << ans << "\n";
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0