#157266 (C++14) No.493 とても長い数列と文字列(Long Long Sequence and a String）

提出ソース

結果

問題	No.493 とても長い数列と文字列(Long Long Sequence and a String）
ユーザー	tossy
提出日時	2017-03-11 11:32:10
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	2 ms / 800 ms
コード長	2,177 bytes
コンパイル時間	1,529 ms
コンパイル使用メモリ	170,436 KB
実行使用メモリ	4,384 KB
最終ジャッジ日時	2023-09-06 17:34:02
合計ジャッジ時間	5,083 ms
ジャッジサーバーID （参考情報）	judge11 / judge13

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 4,380 KB
testcase_01	AC	1 ms 4,380 KB
testcase_02	AC	2 ms 4,376 KB
testcase_03	AC	1 ms 4,376 KB
testcase_04	AC	2 ms 4,376 KB
testcase_05	AC	2 ms 4,376 KB
testcase_06	AC	2 ms 4,376 KB
testcase_07	AC	1 ms 4,380 KB
testcase_08	AC	1 ms 4,376 KB
testcase_09	AC	2 ms 4,376 KB
testcase_10	AC	2 ms 4,380 KB
testcase_11	AC	2 ms 4,380 KB
testcase_12	AC	2 ms 4,380 KB
testcase_13	AC	1 ms 4,380 KB
testcase_14	AC	1 ms 4,376 KB
testcase_15	AC	2 ms 4,376 KB
testcase_16	AC	2 ms 4,376 KB
testcase_17	AC	1 ms 4,376 KB
testcase_18	AC	2 ms 4,376 KB
testcase_19	AC	1 ms 4,376 KB
testcase_20	AC	2 ms 4,380 KB
testcase_21	AC	1 ms 4,376 KB
testcase_22	AC	2 ms 4,376 KB
testcase_23	AC	1 ms 4,380 KB
testcase_24	AC	2 ms 4,376 KB
testcase_25	AC	1 ms 4,376 KB
testcase_26	AC	1 ms 4,380 KB
testcase_27	AC	2 ms 4,376 KB
testcase_28	AC	1 ms 4,380 KB
testcase_29	AC	2 ms 4,376 KB
testcase_30	AC	2 ms 4,384 KB
testcase_31	AC	1 ms 4,376 KB
testcase_32	AC	1 ms 4,380 KB
testcase_33	AC	2 ms 4,380 KB
testcase_34	AC	1 ms 4,380 KB
testcase_35	AC	2 ms 4,376 KB
testcase_36	AC	2 ms 4,376 KB
testcase_37	AC	1 ms 4,384 KB
testcase_38	AC	2 ms 4,376 KB
testcase_39	AC	2 ms 4,376 KB
testcase_40	AC	1 ms 4,380 KB
testcase_41	AC	1 ms 4,380 KB
testcase_42	AC	1 ms 4,380 KB
testcase_43	AC	1 ms 4,376 KB
testcase_44	AC	2 ms 4,380 KB
testcase_45	AC	1 ms 4,376 KB
testcase_46	AC	1 ms 4,376 KB
testcase_47	AC	2 ms 4,380 KB
testcase_48	AC	2 ms 4,376 KB
testcase_49	AC	1 ms 4,380 KB
testcase_50	AC	2 ms 4,380 KB
testcase_51	AC	1 ms 4,380 KB
testcase_52	AC	1 ms 4,380 KB
testcase_53	AC	1 ms 4,380 KB
testcase_54	AC	1 ms 4,380 KB
testcase_55	AC	2 ms 4,376 KB
testcase_56	AC	1 ms 4,376 KB
testcase_57	AC	1 ms 4,376 KB
testcase_58	AC	2 ms 4,380 KB
testcase_59	AC	2 ms 4,380 KB
testcase_60	AC	1 ms 4,380 KB
testcase_61	AC	1 ms 4,376 KB
testcase_62	AC	1 ms 4,380 KB
testcase_63	AC	1 ms 4,380 KB
testcase_64	AC	1 ms 4,384 KB
testcase_65	AC	2 ms 4,376 KB
testcase_66	AC	1 ms 4,376 KB
testcase_67	AC	1 ms 4,380 KB
testcase_68	AC	2 ms 4,380 KB
testcase_69	AC	2 ms 4,380 KB
testcase_70	AC	2 ms 4,380 KB
testcase_71	AC	1 ms 4,380 KB
testcase_72	AC	2 ms 4,376 KB
testcase_73	AC	1 ms 4,376 KB
testcase_74	AC	2 ms 4,376 KB
testcase_75	AC	2 ms 4,380 KB
testcase_76	AC	2 ms 4,376 KB
testcase_77	AC	1 ms 4,380 KB
testcase_78	AC	1 ms 4,376 KB
testcase_79	AC	1 ms 4,376 KB
testcase_80	AC	2 ms 4,376 KB
testcase_81	AC	1 ms 4,376 KB
testcase_82	AC	1 ms 4,376 KB
testcase_83	AC	1 ms 4,380 KB
testcase_84	AC	1 ms 4,380 KB
testcase_85	AC	1 ms 4,376 KB
testcase_86	AC	1 ms 4,380 KB
testcase_87	AC	1 ms 4,380 KB
testcase_88	AC	1 ms 4,380 KB
testcase_89	AC	1 ms 4,384 KB
testcase_90	AC	1 ms 4,376 KB
testcase_91	AC	2 ms 4,380 KB
testcase_92	AC	2 ms 4,376 KB
testcase_93	AC	2 ms 4,380 KB
testcase_94	AC	2 ms 4,376 KB
testcase_95	AC	1 ms 4,380 KB
testcase_96	AC	2 ms 4,380 KB
testcase_97	AC	2 ms 4,380 KB
testcase_98	AC	1 ms 4,376 KB
testcase_99	AC	2 ms 4,376 KB
testcase_100	AC	2 ms 4,376 KB
testcase_101	AC	2 ms 4,376 KB
testcase_102	AC	2 ms 4,376 KB
testcase_103	AC	1 ms 4,380 KB
testcase_104	AC	1 ms 4,376 KB
testcase_105	AC	1 ms 4,380 KB
testcase_106	AC	1 ms 4,380 KB
testcase_107	AC	1 ms 4,380 KB
testcase_108	AC	2 ms 4,376 KB
testcase_109	AC	2 ms 4,376 KB
testcase_110	AC	2 ms 4,380 KB
testcase_111	AC	1 ms 4,380 KB
testcase_112	AC	2 ms 4,376 KB
testcase_113	AC	1 ms 4,376 KB
testcase_114	AC	2 ms 4,380 KB
testcase_115	AC	2 ms 4,380 KB
testcase_116	AC	2 ms 4,380 KB
testcase_117	AC	1 ms 4,380 KB
testcase_118	AC	2 ms 4,376 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std;
#define repl(i,a,b) for(int i=(int)(a);i<(int)(b);i++)
#define rep(i,n) repl(i,0,n)
#define mp(a,b) make_pair((a),(b))
#define pb(a) push_back((a))
#define all(x) (x).begin(),(x).end()
#define uniq(x) sort(all(x)),(x).erase(unique(all(x)),end(x))
#define fi first
#define se second
#define dbg(x) cout<<#x" = "<<((x))<<endl
template<class T,class U> ostream& operator<<(ostream& o, const pair<T,U> &p){o<<"("<<p.fi<<","<<p.se<<")";return o;}
template<class T> ostream& operator<<(ostream& o, const vector<T> &v){o<<"[";for(T t:v){o<<t<<",";}o<<"]";return o;}

#define INF 2147483600
#define MOD 1000000007

// x^n mod p
long mod_pow(long x, long n, long p=MOD){
  long res=1;
  x %= p;
  while(n>0){
    if(n&1) res=res*x%p;
    x=x*x%p;
    n>>=1;
  }
  return res;
}
long mod_inv(long x, long p=MOD){ return mod_pow(x%p, p-2, p); }


#define LIM 61
long len[LIM],add[LIM],mul[LIM];
vector<long> Nstr[LIM]; // i*i を桁ごと

// n*nの桁数返し，nの桁のaddとmulも計算
long num(int n, long *a, long *m, long maxL=INF){
  long l = 0;
  *a = 0;
  *m = 1;
  for(auto d : Nstr[n]){
    if(l>=maxL) return l;
    *a += d;
    *m *= d;
    l++;
  }
  return l;
}

// 1からx桁目までの和と積
pair<long,long> calc(long x){
  long a=0, m=1;
  while(x>0){
    int idx = upper_bound(len, len+LIM, x) - len;
    idx--;
    a = (a + add[idx]);
    m = (m * mul[idx])%MOD;
    x -= len[idx];
    long p,q,r;
    p = num(idx+1, &q, &r, x);
    a = a+q;
    m = (m*r)%MOD;
    x -= p;
  }
  return mp(a,m);
}

int main(){
  long k,l,r;
  cin>>k>>l>>r;
  l--;

  repl(i,1,LIM){
    long x = i*i;
    while(x>0){
      long d = x%10;
      if(d==0) d=10;
      Nstr[i].pb(d);
      x /= 10;
    }
    reverse(all(Nstr[i]));
  }

  len[0]=0;
  add[0]=0;
  mul[0]=1;

  repl(i,1,LIM){
    long ll,a,m;
    ll = num(i, &a, &m);
    len[i] = len[i-1]*2+ll;
    add[i] = add[i-1]*2+a;
    mul[i] = (mul[i-1]*mul[i-1]%MOD * m)%MOD;
  }

  if(k<LIM && len[k] < r){
    cout<<"-1"<<endl; return 0;
  }

  auto lb = calc(l);
  auto ub = calc(r);

  cout << ub.fi-lb.fi << " ";
  cout << ub.se*mod_inv(lb.se)%MOD << endl;

  return 0;
}

yukicoder

結果

テストケース

ソースコード