結果

問題 No.493 とても長い数列と文字列(Long Long Sequence and a String）
コンテスト
ユーザー rickytheta
提出日時 2017-03-10 23:57:14
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 2 ms / 800 ms
コード長 3,066 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,414 ms
コンパイル使用メモリ 163,332 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-24 09:09:06
合計ジャッジ時間 4,208 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 115
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

// mod
const ll MOD = 1000000007ll;
#define FIX(a) ((a)%MOD+MOD)%MOD

ll k,l,r;

ll lensum[125];
ll wasum[125];
ll sekisum[125];
string val[125];

string tos(ll v){
  string ret = "";
  while(v){
    ret = (char)(v%10 + '0') + ret;
    v/=10;
  }
  return ret;
}

inline pll merge(const pll &a, const pll &b){
  return pll(a.first+b.first,a.second*b.second%MOD);
}

pll dfs(int level, ll lft){
  if(l<=lft && lft+lensum[level]<=r){
    pll ret = pll(wasum[level],sekisum[level]);
    return ret;
  }else if(r<=lft || lft+lensum[level]<=l){
    pll ret = pll(0,1);
    return ret;
  }
  pll ret(0,1);
  ret = merge(ret,dfs(level-1,lft));
  lft += lensum[level-1];
  REP(i,val[level].size()){
    if(l<=lft+i && lft+i<r){
      int yo = val[level][i]-'0';
      if(yo==0)yo=10;
      ret = merge(ret,pll(yo,yo));
    }
  }
  lft += val[level].size();
  ret = merge(ret,dfs(level-1,lft));
  return ret;
}

pll naive(){
  string s;
  FOR(i,1,k+1){
    string a = tos(i*i);
    s = s + a + s;
  }
  ll wa = 0, seki = 1;
  FOR(i,l,r){
    int v = s[i]-'0';
    if(v==0)v=10;
    wa += v;
    seki = seki*v%MOD;
  }
  return pll(wa,seki);
}

uint32_t xorshift(void) {
  static uint32_t x = 123456789;
  static uint32_t y = 362436069;
  static uint32_t z = 521288629;
  uint32_t t;

  t = (x ^ (x << 3)) ^ (y ^ (y >> 19)) ^ (z ^ (z << 6));
  x = y; y = z;
  return z = t;
}

int main(){
  cin>>k>>l>>r;
  lensum[0] = 0;
  wasum[0] = 0;
  sekisum[0] = 1;
  val[0] = "";
  FOR(i,1,k+1){
    ll v = i*i;
    val[i] = tos(v);
    lensum[i] = 2*lensum[i-1] + val[i].size();
    wasum[i] = 2*wasum[i-1];
    sekisum[i] = sekisum[i-1]*sekisum[i-1]%MOD;
    REP(j,val[i].size()){
      int yo = val[i][j]-'0';
      if(yo==0)yo=10;
      wasum[i] = wasum[i]+yo;
      sekisum[i] = sekisum[i]*yo%MOD;
    }
    if(lensum[i]>=r){
      k = i;
      break;
    }
  }
  if(lensum[k]<r){
    puts("-1");
    return 0;
  }
  --l;
  // DEBUG(lensum[k]);
  // while(true){
  //   int len = lensum[k];
  //   while(true){
  //     l = xorshift()%len;
  //     r = xorshift()%len;r++;
  //     if(l<r)break;
  //   }
  //   pll ans = dfs(k,0);
  //   pll yo = naive();
  //   if(ans!=yo){
  //     DEBUG(ans.first);
  //     DEBUG(ans.second);
  //     DEBUG(yo.first);
  //     DEBUG(yo.second);
  //     DEBUG(k);
  //     DEBUG(l);
  //     DEBUG(r);
  //     break;
  //   }
  // }
  pll ans = dfs(k,0);
  printf("%lld %lld\n",ans.first,ans.second);
  return 0;
}
0