結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2020-04-24 21:51:44 |
言語 | C++11 (gcc 13.3.0) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,816 bytes |
コンパイル時間 | 739 ms |
コンパイル使用メモリ | 84,276 KB |
実行使用メモリ | 94,252 KB |
最終ジャッジ日時 | 2024-11-07 20:33:51 |
合計ジャッジ時間 | 29,558 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 37 TLE * 4 |
ソースコード
#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 1000000007using 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{int size;vector<vector<llint> > vec;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] = 0; //}void set(int i, llint 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] = gcd(vec[i-1][j], vec[i-1][j+(1<<(i-1))]); //}}}}llint query(int l, int r){int k = 0;for(int i = 0; r-l+1 >= (1<<i); i++) k = i;return gcd(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;}