結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2020-04-24 21:39:42 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 2,877 bytes |
コンパイル時間 | 3,926 ms |
コンパイル使用メモリ | 203,692 KB |
最終ジャッジ日時 | 2025-01-09 23:21:23 |
ジャッジサーバーID (参考情報) |
judge3 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 37 TLE * 4 |
ソースコード
#include <bits/stdc++.h>using namespace std;using int64 = long long;const int mod = 1e9 + 7;// const int mod = 998244353;const int64 infll = (1LL << 62) - 1;const int inf = (1 << 30) - 1;struct IoSetup {IoSetup() {cin.tie(nullptr);ios::sync_with_stdio(false);cout << fixed << setprecision(10);cerr << fixed << setprecision(10);}} iosetup;template< typename T1, typename T2 >ostream &operator<<(ostream &os, const pair< T1, T2 > &p) {os << p.first << " " << p.second;return os;}template< typename T1, typename T2 >istream &operator>>(istream &is, pair< T1, T2 > &p) {is >> p.first >> p.second;return is;}template< typename T >ostream &operator<<(ostream &os, const vector< T > &v) {for(int i = 0; i < (int) v.size(); i++) {os << v[i] << (i + 1 != v.size() ? " " : "");}return os;}template< typename T >istream &operator>>(istream &is, vector< T > &v) {for(T &in : v) is >> in;return is;}template< typename T1, typename T2 >inline bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); }template< typename T1, typename T2 >inline bool chmin(T1 &a, T2 b) { return a > b && (a = b, true); }template< typename T = int64 >vector< T > make_v(size_t a) {return vector< T >(a);}template< typename T, typename... Ts >auto make_v(size_t a, Ts... ts) {return vector< decltype(make_v< T >(ts...)) >(a, make_v< T >(ts...));}template< typename T, typename V >typename enable_if< is_class< T >::value == 0 >::type fill_v(T &t, const V &v) {t = v;}template< typename T, typename V >typename enable_if< is_class< T >::value != 0 >::type fill_v(T &t, const V &v) {for(auto &e : t) fill_v(e, v);}template< typename F >struct FixPoint : F {FixPoint(F &&f) : F(forward< F >(f)) {}template< typename... Args >decltype(auto) operator()(Args &&... args) const {return F::operator()(*this, forward< Args >(args)...);}};template< typename F >inline decltype(auto) MFP(F &&f) {return FixPoint< F >{forward< F >(f)};}int main() {int N;cin >> N;vector< int64 > A(N);cin >> A;A.emplace_back(1);int64 ret = 0;auto dp = make_v< int64 >(21, N + 1);auto nxt = make_v< int >(21, N + 1);fill_v(nxt, N + 1);for(int i = 0; i <= N; i++) {dp[0][i] = A[i];nxt[0][i] = i + 1;}for(int i = 0; i < 20; i++) {for(int j = 0; j <= N; j++) {nxt[i + 1][j] = nxt[i][j];dp[i + 1][j] = dp[i][j];if(nxt[i][j] != N + 1) {nxt[i + 1][j] = nxt[i][nxt[i][j]];dp[i + 1][j] = __gcd(dp[i][j], dp[i][nxt[i][j]]);}}}for(int i = 0; i < N; i++) {int64 po = i, cur = A[i];for(int j = 19; j >= 0; j--) {if(__gcd(cur, dp[j][po]) > 1) {cur = __gcd(cur, dp[j][po]);po = nxt[j][po];}}ret += N - po;}cout << ret << endl;}