結果

問題 No.1339 循環小数
ユーザー mugen_1337mugen_1337
提出日時 2021-01-15 23:06:57
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 744 ms / 2,000 ms
コード長 2,728 bytes
コンパイル時間 2,279 ms
コンパイル使用メモリ 206,468 KB
最終ジャッジ日時 2025-01-17 20:25:14
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define ALL(x) begin(x),end(x)
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}

struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
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==(int)v.size()?"":" ");
    return os;
}
template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}


bool isPrime(ll num){
    if(num<2)           return false;
    else if(num==2)     return true;
    else if(num%2==0)   return false;
    for(ll i=3;i*i<=num;i+=2) if(num%i==0) return false;
    return true;
}

ll pow_mod(ll a,ll n,ll m){
    ll ret=1;
    while(n){
        if(n&1) ret=ret*a%m;
        a=a*a%m;
        n>>=1;
    }
    return ret;
}

template<typename T>
T mod_log(T a,T b,T MOD){
  using ll = long long;
  ll g=1;
  {
    ll m=MOD;
    while(m){
      g=(ll)g*a%MOD;
      m>>=1;
    }
  }
  g=__gcd(g,(ll)MOD);
  ll c=0,t=1;
  while(t%g){
    if(t==b) return c;
    t=t*a%MOD;
    c++;
  }
  // not found
  if(b%g) return MOD;
  t/=g;b/=g;
  const ll n=MOD/g;
  ll h=0,gs=1;
  while(h*h<n){
    gs=gs*a%n;
    ++h;
  }
  unordered_map<ll, ll> bs;
  {
    ll s=0,e=b;
    while(s<h){
      e=e*a%n;
      bs[e]=++s;
    }
  }
  ll ret=LINF;
  {
    ll s=0,e=t;
    while(s<n){
      e=e*gs%n;
      s+=h;
      if(bs.count(e)){
        //   cout<<c+s-bs[e]<<endl;
        // return c+s-bs[e];
        if(c+s-bs[e]>0) chmin(ret,c+s-bs[e]);
      }
    }
  }
  return ret;
}

vector<ll> divisor(ll n) {
    vector<ll> ret;
    for(ll i=1; i*i<=n; i++) {
        if(n % i == 0) {
            ret.push_back(i);
            if(n / i != i) ret.push_back(n/i);
        }
    }
    sort(ALL(ret));
    return ret;
}


void solve(){
    ll n;cin>>n;
    while(n%2==0) n/=2;
    while(n%5==0) n/=5;
    if(n==1){
        cout<<1<<endl;
        return ;
    }

    ll r=mod_log(10ll,1ll,n);
    if(r==LINF){
        cout<<1<<endl;
        return ;
    }


    for(auto x:divisor(r)){
        if(pow_mod(10,x,n)==1){
            cout<<x<<endl;
            return ;
        }
    }
}

signed main(){
    int q;cin>>q;
    while(q--) solve();
    return 0;
}
0