結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2020-04-24 21:52:52 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 426 ms / 2,000 ms |
コード長 | 1,913 bytes |
コンパイル時間 | 2,697 ms |
コンパイル使用メモリ | 197,776 KB |
最終ジャッジ日時 | 2025-01-09 23:35:11 |
ジャッジサーバーID (参考情報) |
judge4 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
コンパイルメッセージ
main.cpp: In function ‘int main()’: main.cpp:98:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 98 | scanf("%lld",&A[i]); | ~~~~~^~~~~~~~~~~~~~
ソースコード
#include <bits/stdc++.h>using namespace std;#define modulo 1000000007#define mod(mod_x) ((((long long)mod_x+modulo))%modulo)#define Inf 1000000000000000001template <typename T,typename F>struct segtree{//元データx[i]はv[n+i]//v[i]の親はv[i/2],子はv[i*2]とv[i*2+1]F func;vector<T> v;int n;T init_value;segtree(int sz,F f,T iv):func(f){init_value = iv;n=1;while(true){if(n>=sz)break;n*=2;}v.resize(2*n,init_value);for(int i=n-1;i>=0;i--){v[i]=func(v[i<<1],v[(i<<1)+1]);}}segtree(vector<T> &x,F f,T iv):func(f){init_value = iv;n=1;while(true){if(n>=x.size())break;n*=2;}v.resize(2*n,init_value);for(int i=0;i<x.size();i++){v[n+i]=x[i];}for(int i=n-1;i>=0;i--){v[i]=func(v[i<<1],v[(i<<1)+1]);}}void update(int x,T val){x+=n;v[x]=val;while(x>0){x>>=1;v[x]=func(v[x<<1],v[(x<<1)+1]);}}//区間[l,r)におけるクエリ処理T query(int l,int r){if(l>=r)return init_value;l+=n;r+=n;T res1 = init_value;T res2 = init_value;while(true){if(l&1){res1=func(res1,v[l++]);}if(r&1){res2=func(v[--r],res2);}if(l>=r)break;l>>=1;r>>=1;}return func(res1,res2);}void show(){int n = 1;for(int i=1;i<v.size();i++){for(int j=0;j<n;j++){if(j!=0)cout<<' ';cout<<v[i+j];}cout<<endl;i+=n-1;n*=2;}}};int main(){long long N;cin>>N;vector<long long> A(N);for(int i=0;i<N;i++){scanf("%lld",&A[i]);}auto f = [](long long a,long long b){return gcd(a,b);};segtree<long long,decltype(f)> seg(A,f,0LL);long long last = 1;long long ans = 0;for(int i=0;i<N;i++){long long X = seg.query(i,last);while(X!=1&&last!=N){X = gcd(X,A[last]);last++;}//cout<<last<<endl;if(X==1)ans += N+1-last;}cout<<ans<<endl;return 0;}