結果

問題 No.1339 循環小数
ユーザー HaarHaar
提出日時 2021-01-15 23:34:31
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 577 ms / 2,000 ms
コード長 3,878 bytes
コンパイル時間 2,297 ms
コンパイル使用メモリ 210,616 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-05-05 02:03:07
合計ジャッジ時間 9,228 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 1 ms
5,376 KB
testcase_11 AC 4 ms
5,376 KB
testcase_12 AC 4 ms
5,376 KB
testcase_13 AC 4 ms
5,376 KB
testcase_14 AC 4 ms
5,376 KB
testcase_15 AC 4 ms
5,376 KB
testcase_16 AC 4 ms
5,376 KB
testcase_17 AC 4 ms
5,376 KB
testcase_18 AC 4 ms
5,376 KB
testcase_19 AC 4 ms
5,376 KB
testcase_20 AC 4 ms
5,376 KB
testcase_21 AC 305 ms
5,376 KB
testcase_22 AC 348 ms
5,376 KB
testcase_23 AC 332 ms
5,376 KB
testcase_24 AC 314 ms
5,376 KB
testcase_25 AC 339 ms
5,376 KB
testcase_26 AC 366 ms
5,376 KB
testcase_27 AC 370 ms
5,376 KB
testcase_28 AC 339 ms
5,376 KB
testcase_29 AC 281 ms
5,376 KB
testcase_30 AC 309 ms
5,376 KB
testcase_31 AC 506 ms
5,376 KB
testcase_32 AC 512 ms
5,376 KB
testcase_33 AC 332 ms
5,376 KB
testcase_34 AC 187 ms
5,376 KB
testcase_35 AC 577 ms
5,376 KB
testcase_36 AC 330 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#ifdef DEBUG
#include <Mylib/Debug/debug.cpp>
#else
#define dump(...) ((void)0)
#endif

template <typename T, typename U>
bool chmin(T &a, const U &b){
  return (a > b ? a = b, true : false);
}

template <typename T, typename U>
bool chmax(T &a, const U &b){
  return (a < b ? a = b, true : false);
}

template <typename T, size_t N, typename U>
void fill_array(T (&a)[N], const U &v){
  std::fill((U*)a, (U*)(a + N), v);
}

template <typename T, size_t N, size_t I = N>
auto make_vector(const std::array<int, N> &a, T value = T()){
  static_assert(I >= 1);
  static_assert(N >= 1);
  if constexpr (I == 1){
    return std::vector<T>(a[N - I], value);
  }else{
    return std::vector(a[N - I], make_vector<T, N, I - 1>(a, value));
  }
}

template <typename T>
std::ostream& operator<<(std::ostream &s, const std::vector<T> &a){
  for(auto it = a.begin(); it != a.end(); ++it){
    if(it != a.begin()) s << " ";
    s << *it;
  }
  return s;
}

template <typename T>
std::istream& operator>>(std::istream &s, std::vector<T> &a){
  for(auto &x : a) s >> x;
  return s;
}

std::string YesNo(bool value){return value ? "Yes" : "No";}
std::string YESNO(bool value){return value ? "YES" : "NO";}
std::string yesno(bool value){return value ? "yes" : "no";}

template <typename T>
void putl(const T &value){
  std::cout << value << "\n";
}

template <typename Head, typename ... Tail>
void putl(const Head head, const Tail &... tail){
  std::cout << head << " ";
  putl(tail ...);
}

namespace haar_lib {
  constexpr int64_t mod_pow(int64_t n, int64_t p, int64_t m){
    int64_t ret = 1;
    while(p > 0){
      if(p & 1) (ret *= n) %= m;
      (n *= n) %= m;
      p >>= 1;
    }
    return ret;
  }
}

namespace haar_lib {
  int64_t mod_inv(int64_t a, int64_t m){
    int64_t b = m, u = 1, v = 0;

    while(b){
      int64_t t = a / b;
      a -= t * b; std::swap(a, b);
      u -= t * v; std::swap(u, v);
    }

    u %= m;
    if(u < 0) u += m;

    return u;
  }
}

namespace haar_lib {
  std::optional<int64_t> mod_log(int64_t a, int64_t b, int64_t m){
    if(b == 1) return 0;

    int64_t d = 0;

    while(1){
      if(auto g = std::gcd(a, m); g != 1){
        if(b % g != 0) return {};

        d += 1;
        m /= g;
        b /= g;
        (b *= mod_inv(a / g, m)) %= m;

        if(b == 1) return d;
      }else{
        break;
      }
    }

    const int64_t sq = std::sqrt(m) + 1;

    std::unordered_map<int64_t, int64_t> mp;
    {
      int64_t t = 1 % m;

      for(int i = 0; i < sq; ++i){
        if(mp.find(t) == mp.end()) mp[t] = i;
        (t *= a) %= m;
      }
    }

    {
      int64_t A = mod_pow(mod_inv(a, m), sq, m);
      int64_t t = b % m;

      for(int i = 0; i < sq; ++i){
        if(mp.find(t) != mp.end()){
          int64_t ret = i * sq + mp[t] + d;

          return ret;
        }

        (t *= A) %= m;
      }
    }

    return {};
  }
}


namespace haar_lib {}

namespace solver {
  using namespace haar_lib;

  constexpr int m1000000007 = 1000000007;
  constexpr int m998244353 = 998244353;

  void init(){
    std::cin.tie(0);
    std::ios::sync_with_stdio(false);
    std::cout << std::fixed << std::setprecision(12);
    std::cerr << std::fixed << std::setprecision(12);
    std::cin.exceptions(std::ios_base::failbit);
  }

  int repeating_decimal(int N){
    while(N % 2 == 0) N /= 2;
    while(N % 5 == 0) N /= 5;

    if(N == 1){
      return 1;
    }

    return *mod_log(10, mod_inv(10, N), N) + 1;
  }

  void solve(){
    int T; std::cin >> T;
    while(T--){
      int N; std::cin >> N;
      std::cout << repeating_decimal(N) << "\n";
    }
  }
}

int main(){
  solver::init();
  while(true){
    try{
      solver::solve();
      std::cout << std::flush;
      std::cerr << std::flush;
    }catch(const std::istream::failure &e){
      break;
    }catch(...){
      break;
    }
  }
  return 0;
}
0