結果

問題 No.1036 Make One With GCD 2
ユーザー fura
提出日時 2021-05-28 17:29:31
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 304 ms / 2,000 ms
コード長 1,255 bytes
コンパイル時間 2,216 ms
コンパイル使用メモリ 202,336 KB
最終ジャッジ日時 2025-01-21 18:57:35
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:50:21: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   50 |         int n; scanf("%d",&n);
      |                ~~~~~^~~~~~~~~
main.cpp:52:23: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   52 |         rep(i,n) scanf("%lld",&a[i]);
      |                  ~~~~~^~~~~~~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;
using lint=long long;

template<class B>
class sparse_table{
	vector<vector<B>> st;
	vector<int> h;
public:
	sparse_table()=default;
	template<class T>
	sparse_table(const vector<T>& a){ build(a); }

	template<class T>
	void build(const vector<T>& a){
		int n=a.size();
		h.assign(n+1,0);
		for(int i=2;i<=n;i++) h[i]=h[i>>1]+1;
		st.resize(h[n]+1);
		st[0].resize(n);
		rep(i,n) st[0][i]=a[i];
		for(int i=1;i<h[n]+1;i++){
			st[i].resize(n-(1<<i)+1);
			rep(j,n-(1<<i)+1) st[i][j]=st[i-1][j]*st[i-1][j+(1<<(i-1))];
		}
	}

	B product(int l,int r)const{
		int i=h[r-l];
		return st[i][l]*st[i][r-(1<<i)];
	}
};

template<class T>
class gcd_monoid{
	T a;
public:
	gcd_monoid(const T& val=T()):a(val){}
	gcd_monoid operator*(const gcd_monoid& x)const{
		return gcd(a,x.a);
	}
	T& get(){ return a; }
	const T& get()const{ return a; }
};

int main(){
	int n; scanf("%d",&n);
	vector<lint> a(n);
	rep(i,n) scanf("%lld",&a[i]);

	sparse_table<gcd_monoid<lint>> ST(a);

	lint ans=0;
	for(int l=0,r=0;l<n;l++){
		r=max(r,l+1);
		while(r<n && ST.product(l,r).get()>1) r++;
		if(r==n && ST.product(l,r).get()>1) continue;
		ans+=n-r+1;
	}
	printf("%lld\n",ans);

	return 0;
}
0