結果

問題 No.1036 Make One With GCD 2
ユーザー leaf_1415
提出日時 2021-04-16 05:14:23
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 2,000 bytes
コンパイル時間 1,162 ms
コンパイル使用メモリ 86,708 KB
実行使用メモリ 96,384 KB
最終ジャッジ日時 2024-07-02 11:51:04
合計ジャッジ時間 23,578 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 37 TLE * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>
#include <cassert>
#include <vector>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <bitset>
#include <string>
#include <algorithm>
#include <utility>
#define llint long long
#define inf 1e18
#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)
#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))
#define mod 1000000007

using namespace std;
typedef pair<llint, llint> P;

llint gcd(llint a, llint b)
{
	if(b == 0) return a;
	return gcd(b, a%b);
}

struct SparseTable{
	typedef llint T;
	
	int size;
	vector<vector<T> > vec;
	vector<int> msb;
	
	T Ident(){return 0;}  //単位元
	T ope(T a, T b){return gcd(a, b);} //作用素
	
	SparseTable(){}
	SparseTable(int size){
		this->size = size;
		vec.resize(size+1);
		for(int i = 0; i <= size; i++) vec[i].resize(1<<size);
		for(int i = 0; i < (1<<size); i++) vec[0][i] = Ident();
		msb.resize((1<<size));
		for(int i = 2; i < (1<<size); i++) msb[i] = msb[i>>1]+1;
	}
	void set(int i, T val){
		vec[0][i] = val;
	}
	void calc(){
		for(int i = 1; i <= size; i++){
			for(int j = 0; j < (1<<size); j++){
				if(j+(1<<(i-1)) < (1<<size)){
					vec[i][j] = ope(vec[i-1][j], vec[i-1][j+(1<<(i-1))]);
				}
			}
		}
	}
	T query(int l, int r)
	{
		if(l > r) return Ident();
		int k = msb[r-l+1];
		return ope(vec[k][l], vec[k][r-(1<<k)+1]);
	}
};


llint n;
llint a[500005];
SparseTable sp(19);

int main(void)
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	
	cin >> n;
	for(int i = 1; i <= n; i++) cin >> a[i];
	for(int i = 1; i <= n; i++) sp.set(i, a[i]);
	sp.calc();
	
	llint ans = 0;
	for(int i = 1; i <= n; i++){
		llint ub = n+1, lb = i-1, mid;
		while(ub-lb>1){
			mid = (ub+lb)/2;
			if(mid >= i && sp.query(i, mid) == 1) ub = mid;
			else lb = mid;
		}
		ans += n-ub+1;
	}
	cout << ans << endl;
	
	return 0;
}
0