結果

問題 No.1036 Make One With GCD 2
ユーザー 沙耶花
提出日時 2020-04-25 19:57:00
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 636 ms / 2,000 ms
コード長 1,188 bytes
コンパイル時間 2,838 ms
コンパイル使用メモリ 197,340 KB
最終ジャッジ日時 2025-01-10 01:27:32
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:52:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   52 |                 scanf("%lld",&A[i]);
      |                 ~~~~~^~~~~~~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define modulo 1000000007
#define mod(mod_x) ((((long long)mod_x+modulo))%modulo)
#define Inf 1000000000000000001

template <typename T>
struct sparse_table{
	const T init_value = 0;
	
	vector<T> v;
	int n;
	int sz;
	sparse_table(vector<T> &x){
		sz = x.size();
		for(int i=0;true;i++){
			if((1<<i)>sz){
				n = i;
				break;
			}
		}
		v.resize(sz*n,init_value);
		
		for(int i=0;i<sz;i++){
			v[i] = x[i];
		}
		
		for(int i=1;i<n;i++){
			for(int j=0;j<sz+1-(1<<i);j++){
				int temp = j+(1<<(i-1));
				v[i*sz+j] = func(v[(i-1)*sz+j],v[(i-1)*sz+temp]);
			}
		}
	}
	
	T query(int l,int r){
		int x = 31-__builtin_clz(r-l);
		return func(v[sz*x+l],v[sz*x+r-(1<<x)]);
	}
	
	T func(T a,T b){
		return gcd(a,b);//min(a,b);
	}
};

int main(){
	
	long long N;
	cin>>N;
	vector<long long> A(N);
	for(int i=0;i<N;i++){
		scanf("%lld",&A[i]);
	}
	
	sparse_table<long long> seg(A);
	
	long long ans = 0;
	
	for(int i=0;i<N;i++){
		if(seg.query(i,N)>1)break;
		int ok = i;
		int ng = N+1;
		
		while(ng-ok>1){
			int mid = (ok+ng)/2;
			if(seg.query(i,mid)==1)ng = mid;
			else ok=mid;
		}
		
		ans += N+1-ng;
	}
	
	cout<<ans<<endl;
		
	return 0;
}
0