結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
|
提出日時 | 2020-04-25 11:03:23 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,037 ms / 2,000 ms |
コード長 | 2,042 bytes |
コンパイル時間 | 3,339 ms |
コンパイル使用メモリ | 214,508 KB |
最終ジャッジ日時 | 2025-01-10 00:59:11 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef long long ll; typedef long double ld; typedef pair<int,int> P; typedef pair<ll,ll> l_l; typedef vector<int> vi; typedef vector<vi> vvi; typedef vector<ll> vl; typedef vector<P> vp; #define pb push_back #define sz(x) (int)(x).size() #define fi first #define se second #define rep(i,n) for(int i = 0; i < (n); ++i) #define rrep(i,n) for(int i = 1; i <= (n); ++i) #define drep(i,n) for(int i = (n)-1; i >= 0; --i) const int INF=1001001000; const int mINF=-1001001000; const ll LINF=1010010010010010000; template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } //グリッド:(典型)dp,dfs,bfs,最短経路,その他 ll gcd(ll a, ll b) { return b?gcd(b,a%b):a;} struct SegmentTree{ int n; vector<ll> node; SegmentTree(vector<ll> v){ int sz=v.size(); n=1; while(n<sz) n*=2; node.resize(2*n-1,-1); for(int i=0;i<sz;i++) node[i+n-1]=v[i]; for(int i=n-2;i>=0;i--){ if(node[2*i+1]==-1&&node[2*i+2]==-1)node[i]=-1; else if(node[2*i+1]==-1&&node[2*i+2]!=-1) node[i]=node[2*i+2]; else if(node[2*i+1]!=-1&&node[2*i+2]==-1)node[i]=node[2*i+1]; else node[i]=gcd(node[2*i+1],node[2*i+2]); } } void update(ll x,ll val){ x+=(n-1); node[x]=val; while(x>0){ x=(x-1)/2; node[x]=min(node[2*x+1],node[2*x+2]); } } ll getgcd(ll a,ll b,ll k=0,ll l=0,ll r=-1){ if(r<0) r=n; if(r<=a||b<=l) return -1; if(a<=l&&r<=b) return node[k]; ll v1=getgcd(a,b,2*k+1,l,(r+l)/2); ll v2=getgcd(a,b,2*k+2,(r+l)/2,r); if(v1==-1&&v2!=-1)return v2; if(v1!=-1&&v2==-1)return v1; return gcd(v1,v2); } }; int main(){ ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0); ll n;cin>>n; vl a(n);rep(i,n)cin>>a[i]; SegmentTree seg(a); ll j=1; ll ans=0; rep(i,n){ while(seg.getgcd(i,j)!=1&&j<=n){ j++; } if(j==n&&seg.getgcd(i,n)!=1)break; ans+=n-j+1; } cout<<ans<<endl; return 0; }