結果

問題 No.1036 Make One With GCD 2
ユーザー 👑 obakyanobakyan
提出日時 2021-08-02 23:11:56
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 410 ms / 2,000 ms
コード長 1,676 bytes
コンパイル時間 605 ms
コンパイル使用メモリ 71,668 KB
実行使用メモリ 11,232 KB
最終ジャッジ日時 2023-10-14 20:59:40
合計ジャッジ時間 11,022 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 102 ms
9,428 KB
testcase_01 AC 85 ms
9,116 KB
testcase_02 AC 110 ms
11,152 KB
testcase_03 AC 13 ms
7,416 KB
testcase_04 AC 21 ms
7,744 KB
testcase_05 AC 2 ms
5,464 KB
testcase_06 AC 2 ms
5,460 KB
testcase_07 AC 32 ms
7,612 KB
testcase_08 AC 26 ms
7,536 KB
testcase_09 AC 232 ms
9,252 KB
testcase_10 AC 216 ms
8,960 KB
testcase_11 AC 237 ms
9,220 KB
testcase_12 AC 218 ms
8,792 KB
testcase_13 AC 405 ms
8,952 KB
testcase_14 AC 410 ms
8,948 KB
testcase_15 AC 385 ms
8,972 KB
testcase_16 AC 387 ms
8,768 KB
testcase_17 AC 399 ms
8,944 KB
testcase_18 AC 2 ms
5,364 KB
testcase_19 AC 3 ms
5,444 KB
testcase_20 AC 2 ms
5,408 KB
testcase_21 AC 2 ms
5,448 KB
testcase_22 AC 379 ms
8,792 KB
testcase_23 AC 279 ms
7,748 KB
testcase_24 AC 393 ms
8,884 KB
testcase_25 AC 359 ms
8,504 KB
testcase_26 AC 373 ms
8,660 KB
testcase_27 AC 2 ms
5,704 KB
testcase_28 AC 2 ms
5,388 KB
testcase_29 AC 2 ms
5,380 KB
testcase_30 AC 2 ms
5,516 KB
testcase_31 AC 2 ms
5,520 KB
testcase_32 AC 2 ms
5,560 KB
testcase_33 AC 2 ms
5,520 KB
testcase_34 AC 2 ms
5,532 KB
testcase_35 AC 2 ms
5,668 KB
testcase_36 AC 2 ms
5,452 KB
testcase_37 AC 2 ms
5,412 KB
testcase_38 AC 106 ms
11,232 KB
testcase_39 AC 108 ms
9,400 KB
testcase_40 AC 278 ms
7,884 KB
testcase_41 AC 270 ms
9,140 KB
testcase_42 AC 271 ms
9,184 KB
testcase_43 AC 236 ms
10,144 KB
testcase_44 AC 257 ms
9,656 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <vector>
#include <cassert>
#include <iostream>
#include <cstdint>
#include <chrono>
using namespace std;
#define I64 int64_t
#define curtime chrono::system_clock::time_point
#define getcurtime chrono::system_clock::now

I64 getgcd(I64 x, I64 y)
{
  I64 t;
  while(0 < x){
    t = x;
    x = y % x;
    y = t;
  }
  return y;
}

I64 a[500010];
int suf_spos;
I64 pre_sum[500010];
I64 suf_sum;
int pre_idx_len;
int suf_len;

void push(int idx)
{
  suf_len++;
  if(suf_len == 1){
    suf_spos = idx;
    suf_sum = a[idx];
  } else {
    suf_sum = getgcd(suf_sum, a[idx]);
  }
}
void pop()
{
  if(pre_idx_len == 0){
    I64 v = a[suf_spos + suf_len - 1];
    pre_sum[0] = v;
    for(int i = 2; i < suf_len; i++){
      v = getgcd(v, a[suf_spos + suf_len - i]);
      pre_sum[i - 1] = v;
    }
    pre_idx_len = suf_len - 1;
    suf_len = 0;
  } else {
    pre_idx_len--;
  }
}
I64 query()
{
  if(pre_idx_len == 0){
    return suf_sum;
  } else if(suf_len == 0){
    return pre_sum[pre_idx_len - 1];
  } else {
    return getgcd(suf_sum, pre_sum[pre_idx_len - 1]);
  }
}

int main(void)
{
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  int n;
  cin >> n;
  for(int i = 0; i < n; i++){
    I64 v;
    cin >> v;
    a[i] = v;
  }
  a[n] = 1;
  suf_spos = 0;
  suf_sum = 0;
  pre_idx_len = 0;
  suf_len = 0;

  I64 ret = 0;
  I64 cur = a[0];
  int right = 0;
  push(0);
  for(int i = 0; i < n; i++){
    if(right < i){
      push(i);
      right = i;
      cur = a[i];
    }
    cur = query();
    while(1 < cur){
      right++;
      push(right);
      cur = getgcd(cur, a[right]);
    }
    ret += n - right;
    pop();
  }
  cout << ret << endl;
}
0