結果
問題 | 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;}