結果

問題 No.1036 Make One With GCD 2
ユーザー 👑 jupirojupiro
提出日時 2020-04-24 21:44:23
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 805 ms / 2,000 ms
コード長 3,456 bytes
コンパイル時間 1,603 ms
コンパイル使用メモリ 134,640 KB
実行使用メモリ 15,556 KB
最終ジャッジ日時 2024-04-24 17:26:16
合計ジャッジ時間 13,824 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 805 ms
15,556 KB
testcase_01 AC 203 ms
15,500 KB
testcase_02 AC 103 ms
15,476 KB
testcase_03 AC 22 ms
8,840 KB
testcase_04 AC 38 ms
14,092 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 64 ms
9,176 KB
testcase_08 AC 47 ms
8,796 KB
testcase_09 AC 249 ms
15,352 KB
testcase_10 AC 240 ms
15,220 KB
testcase_11 AC 251 ms
15,408 KB
testcase_12 AC 248 ms
15,104 KB
testcase_13 AC 327 ms
15,232 KB
testcase_14 AC 347 ms
15,360 KB
testcase_15 AC 304 ms
15,092 KB
testcase_16 AC 344 ms
15,232 KB
testcase_17 AC 310 ms
15,232 KB
testcase_18 AC 2 ms
5,376 KB
testcase_19 AC 2 ms
5,376 KB
testcase_20 AC 3 ms
5,376 KB
testcase_21 AC 2 ms
5,376 KB
testcase_22 AC 304 ms
15,104 KB
testcase_23 AC 238 ms
14,148 KB
testcase_24 AC 330 ms
15,192 KB
testcase_25 AC 287 ms
14,848 KB
testcase_26 AC 303 ms
14,956 KB
testcase_27 AC 2 ms
5,376 KB
testcase_28 AC 1 ms
5,376 KB
testcase_29 AC 2 ms
5,376 KB
testcase_30 AC 2 ms
5,376 KB
testcase_31 AC 2 ms
5,376 KB
testcase_32 AC 2 ms
5,376 KB
testcase_33 AC 2 ms
5,376 KB
testcase_34 AC 1 ms
5,376 KB
testcase_35 AC 2 ms
5,376 KB
testcase_36 AC 1 ms
5,376 KB
testcase_37 AC 2 ms
5,376 KB
testcase_38 AC 94 ms
15,544 KB
testcase_39 AC 658 ms
15,488 KB
testcase_40 AC 244 ms
14,080 KB
testcase_41 AC 586 ms
15,344 KB
testcase_42 AC 506 ms
15,428 KB
testcase_43 AC 536 ms
15,552 KB
testcase_44 AC 580 ms
15,440 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <iostream>
#include <string>
#include <sstream>
#include <stack>
#include <algorithm>
#include <cmath>
#include <queue>
#include <map>
#include <set>
#include <cstdlib>
#include <bitset>
#include <tuple>
#include <assert.h>
#include <deque>
#include <bitset>
#include <iomanip>
#include <limits>
#include <chrono>
#include <random>
#include <array>
#include <unordered_map>
#include <functional>
#include <complex>
#include <numeric>

template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }

constexpr long long MAX = 5100000;
constexpr long long INF = 1LL << 60;
constexpr int inf = 1000000007;
constexpr long long mod = 1000000007LL;
//constexpr long long mod = 998244353LL;

using namespace std;
typedef unsigned long long ull;
typedef long long ll;

template< typename Monoid >
struct SegmentTree {
    using F = function< Monoid(Monoid, Monoid) >;

    int sz;
    vector< Monoid > seg;

    const F f;
    const Monoid M1;

    SegmentTree(int n, const F f, const Monoid& M1) : f(f), M1(M1) {
        sz = 1;
        while (sz < n) sz <<= 1;
        seg.assign(2 * sz, M1);
    }

    void set(int k, const Monoid& x) {
        seg[k + sz] = x;
    }

    void build() {
        for (int k = sz - 1; k > 0; k--) {
            seg[k] = f(seg[2 * k + 0], seg[2 * k + 1]);
        }
    }

    void update(int k, const Monoid& x) {
        k += sz;
        seg[k] = x;
        while (k >>= 1) {
            seg[k] = f(seg[2 * k + 0], seg[2 * k + 1]);
        }
    }

    Monoid query(int a, int b) {
        Monoid L = M1, R = M1;
        if (a >= b) return M1;
        for (a += sz, b += sz; a < b; a >>= 1, b >>= 1) {
            if (a & 1) L = f(L, seg[a++]);
            if (b & 1) R = f(seg[--b], R);
        }
        return f(L, R);
    }

    template< typename C >
    int find_subtree(int a, const C& check, Monoid& M, bool type) {
        while (a < sz) {
            Monoid nxt = type ? f(seg[2 * a + type], M) : f(M, seg[2 * a + type]);
            if (check(nxt)) a = 2 * a + type;
            else M = nxt, a = 2 * a + 1 - type;
        }
        return a - sz;
    }


    template< typename C >
    int find_first(int a, const C& check) {
        Monoid L = M1;
        if (a <= 0) {
            if (check(f(L, seg[1]))) return find_subtree(1, check, L, false);
            return -1;
        }
        int b = sz;
        for (a += sz, b += sz; a < b; a >>= 1, b >>= 1) {
            if (a & 1) {
                Monoid nxt = f(L, seg[a]);
                if (check(nxt)) return find_subtree(a, check, L, false);
                L = nxt;
                ++a;
            }
        }
        return -1;
    }
};

ll f(ll a, ll b) {
    return gcd(a, b);
}

bool C(ll a) {
    return a == 1;
}

int main()
{
	/*
	cin.tie(nullptr);
	ios::sync_with_stdio(false);
	*/

	ll N; scanf("%lld", &N);
	vector<ll> a(N); for (int i = 0; i < N; i++) scanf("%lld", &a[i]);
    SegmentTree<ll> seg(N, f, 0);
    for (int i = 0; i < N; i++) seg.set(i, a[i]);
    seg.build();
    ll res = 0;
    for (int i = 0; i < N; i++) {
        if (a[i] == 1) {
            res += N - i;
            continue;
        }
        if (seg.query(i, N) != 1) break;
        ll right = seg.find_first(i, C);
        res += N - right;
        //cout << res << endl;
    }
    cout << res << endl;
	return 0;
}
0