#691351 (C++23(draft)) No.1653 Squarefree

• 問題: No.1653 Squarefree
• ユーザー: KoD
• 提出日時: 2021-08-20 22:48:56
• 言語: C++23(draft) (gcc 13.2.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 1,140 ms / 2,000 ms
• コード長: 3,680 bytes
• コンパイル時間: 3,851 ms
• コンパイル使用メモリ: 277,248 KB
• 実行使用メモリ: 6,944 KB
• 最終ジャッジ日時: 2024-04-22 09:13:44
• 合計ジャッジ時間: 34,258 ms
• ジャッジサーバーID （参考情報）: judge4 / judge1

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1,131 ms; 5,248 KB
testcase_01	AC	23 ms; 5,248 KB
testcase_02	AC	23 ms; 5,376 KB
testcase_03	AC	23 ms; 5,376 KB
testcase_04	AC	23 ms; 5,376 KB
testcase_05	AC	23 ms; 5,376 KB
testcase_06	AC	23 ms; 5,376 KB
testcase_07	AC	23 ms; 5,376 KB
testcase_08	AC	23 ms; 5,376 KB
testcase_09	AC	23 ms; 5,376 KB
testcase_10	AC	24 ms; 5,376 KB
testcase_11	AC	293 ms; 5,376 KB
testcase_12	AC	292 ms; 5,376 KB
testcase_13	AC	1,125 ms; 5,376 KB
testcase_14	AC	1,122 ms; 5,376 KB
testcase_15	AC	1,134 ms; 5,376 KB
testcase_16	AC	1,129 ms; 5,376 KB
testcase_17	AC	1,131 ms; 5,376 KB
testcase_18	AC	1,136 ms; 5,376 KB
testcase_19	AC	1,130 ms; 5,376 KB
testcase_20	AC	1,134 ms; 5,376 KB
testcase_21	AC	1,138 ms; 5,376 KB
testcase_22	AC	1,134 ms; 5,376 KB
testcase_23	AC	1,122 ms; 5,376 KB
testcase_24	AC	1,128 ms; 5,376 KB
testcase_25	AC	1,131 ms; 5,376 KB
testcase_26	AC	1,123 ms; 5,376 KB
testcase_27	AC	1,132 ms; 5,376 KB
testcase_28	AC	1,120 ms; 5,376 KB
testcase_29	AC	1,131 ms; 6,940 KB
testcase_30	AC	1,121 ms; 6,940 KB
testcase_31	AC	1,129 ms; 6,940 KB
testcase_32	AC	1,140 ms; 6,944 KB
testcase_33	AC	1,126 ms; 6,940 KB
testcase_34	AC	1,131 ms; 6,940 KB
testcase_35	AC	1,130 ms; 6,944 KB
testcase_36	AC	23 ms; 6,944 KB
testcase_37	AC	23 ms; 6,944 KB
testcase_38	AC	23 ms; 6,940 KB
testcase_39	AC	23 ms; 6,940 KB
testcase_40	AC	23 ms; 6,940 KB

ソースコード

#include <bits/stdc++.h>

using i32 = std::int32_t;
using u32 = std::uint32_t;
using i64 = std::int64_t;
using u64 = std::uint64_t;
using i128 = __int128_t;
using u128 = __uint128_t;
using isize = std::ptrdiff_t;
using usize = std::size_t;

class rep {
    struct Iter {
        usize itr;
        constexpr Iter(const usize pos) noexcept : itr(pos) {}
        constexpr void operator++() noexcept { ++itr; }
        constexpr bool operator!=(const Iter& other) const noexcept { return itr != other.itr; }
        constexpr usize operator*() const noexcept { return itr; }
    };
    const Iter first, last;

  public:
    explicit constexpr rep(const usize first, const usize last) noexcept : first(first), last(std::max(first, last)) {}
    constexpr Iter begin() const noexcept { return first; }
    constexpr Iter end() const noexcept { return last; }
};

constexpr u64 ceil_log2(const u64 x) {
    u64 e = 0;
    while (((u64)1 << e) < x) ++e;
    return e;
}

template <class F> class AutoReallocation {
    using R = typename decltype(std::declval<F>()((usize)0))::value_type;

    F func;
    mutable std::vector<R> data;

  public:
    explicit AutoReallocation(F&& f) : func(std::forward<F>(f)), data() {}

    void reserve(const usize size) const {
        if (data.size() < size) data = func(((usize)1 << ceil_log2(size)));
    }
    R operator[](const usize i) const {
        reserve(i + 1);
        return data[i];
    }
};

template <class F> decltype(auto) auto_realloc(F&& f) {
    using G = std::decay_t<F>;
    return AutoReallocation<G>(std::forward<G>(f));
}

struct PrimeSieve {
    static inline const auto min_prime = auto_realloc([](const usize n) {
        std::vector<usize> ret(n);
        std::iota(ret.begin(), ret.end(), (usize)0);
        std::vector<usize> list;
        for (const usize i : rep(2, n)) {
            if (ret[i] == i) list.push_back(i);
            for (const usize p : list) {
                if (p * i >= n || p > ret[i]) break;
                ret[p * i] = p;
            }
        }
        return ret;
    });
    static bool is_prime(const usize n) {
        if (n <= 1) return false;
        return min_prime[n] == n;
    }
    template <class T> static std::vector<std::pair<T, usize>> factorize(T x) {
        assert(x > 0);
        std::vector<std::pair<T, usize>> ret;
        while (x != 1) {
            const usize p = min_prime[x];
            ret.emplace_back((T)p, 0);
            while (min_prime[x] == p) {
                ret.back().second++;
                x /= p;
            }
        }
        return ret;
    }
};

template <class T> using Vec = std::vector<T>;

i64 isqrt(const i64 X) {
    if (X <= 1) {
        return X;
    }
    i64 ok = 0, ng = X;
    while (ng - ok > 1) {
        const i128 md = (ok + ng) / 2;
        (md * md > (i128)X ? ng : ok) = md;
    }
    return ok;
}

bool is_sq(const i64 x) {
    const i64 y = isqrt(x);
    return y != 1 and y * y == x;
}

void main_() {
    i64 L, R;
    std::cin >> L >> R;
    Vec<char> ok(R - L + 1, true);
    for (i64 k = 2; k <= 1000000; ++k) {
        const i64 x = k * k;
        for (i64 i = (L + x - 1) / x * x; i <= R; i += x) {
            ok[i - L] = false;
        }
    }
    for (i64 k = 1; k <= 1000000; ++k) {
        for (i64 i = (L + k - 1) / k * k; i <= R; i += k) {
            if (ok[i - L] and is_sq(i / k)) {
                ok[i - L] = false;
            }
        }
    }
    usize cnt = 0;
    for (const auto c : ok) {
        if (c) {
            cnt += 1;
        }
    }
    std::cout << cnt << '\n';
}

int main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(nullptr);
    main_();
    return 0;
}