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