結果

問題 No.315 世界のなんとか3.5
ユーザー rickythetarickytheta
提出日時 2015-12-08 17:28:13
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 336 ms / 2,000 ms
コード長 4,171 bytes
コンパイル時間 1,827 ms
コンパイル使用メモリ 180,632 KB
実行使用メモリ 37,660 KB
最終ジャッジ日時 2024-09-14 20:21:45
合計ジャッジ時間 12,197 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 185 ms
36,864 KB
testcase_01 AC 187 ms
36,992 KB
testcase_02 AC 191 ms
36,812 KB
testcase_03 AC 189 ms
36,864 KB
testcase_04 AC 189 ms
36,864 KB
testcase_05 AC 191 ms
36,736 KB
testcase_06 AC 194 ms
36,864 KB
testcase_07 AC 205 ms
36,736 KB
testcase_08 AC 197 ms
36,736 KB
testcase_09 AC 190 ms
36,864 KB
testcase_10 AC 189 ms
36,864 KB
testcase_11 AC 188 ms
36,864 KB
testcase_12 AC 223 ms
37,120 KB
testcase_13 AC 228 ms
37,248 KB
testcase_14 AC 256 ms
37,264 KB
testcase_15 AC 254 ms
37,248 KB
testcase_16 AC 256 ms
37,268 KB
testcase_17 AC 255 ms
37,268 KB
testcase_18 AC 223 ms
36,992 KB
testcase_19 AC 225 ms
36,992 KB
testcase_20 AC 222 ms
37,120 KB
testcase_21 AC 227 ms
37,120 KB
testcase_22 AC 255 ms
37,136 KB
testcase_23 AC 255 ms
37,264 KB
testcase_24 AC 255 ms
37,260 KB
testcase_25 AC 253 ms
37,140 KB
testcase_26 AC 256 ms
37,260 KB
testcase_27 AC 260 ms
37,136 KB
testcase_28 AC 257 ms
37,260 KB
testcase_29 AC 316 ms
37,656 KB
testcase_30 AC 317 ms
37,532 KB
testcase_31 AC 319 ms
37,524 KB
testcase_32 AC 321 ms
37,528 KB
testcase_33 AC 255 ms
37,504 KB
testcase_34 AC 320 ms
37,660 KB
testcase_35 AC 336 ms
37,656 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<vvl> vvvl;
typedef vector<vvvl> vvvvl;
typedef vector<vvvvl> vvvvvl;
typedef vector<vvvvvl> vvvvvvl;
typedef complex<double> P;
typedef pair<int,int> pii;
#define REP(i,n) for(ll i=0;i<n;++i)
#define REPR(i,n) for(ll i=1;i<n;++i)
#define FOR(i,a,b) for(ll i=a;i<b;++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 MOD (ll)(1e9+7)
#define ADD(a,b) a=((a)+(b))%MOD
#define FIX(a) ((a)%MOD+MOD)%MOD

ll calc(string s,int p){
  ll n = s.size();
  while((n=s.size())<5)s="0"+s;

  vvvvl dp(2,           // pos
      vvvl(3,           // mod3
       vvl(2,           // found 3 flag
        vl(2,0))));     // less flag
  int now = 0;
  int nxt = 1;
  dp[0][0][0][0] = 1;
  REP(pos,n-5){
    int d = s[pos]-'0';
    REP(md,3){
      REP(t,10){
        ll nmd = (md*10+t)%3;
        if(t==3){
          // necessarily aho
          dp[nxt][nmd][1][1] += dp[now][md][1][1] + dp[now][md][0][1];
          if(t < d) dp[nxt][nmd][1][1] += dp[now][md][1][0] + dp[now][md][0][0];
          else if(t == d) dp[nxt][nmd][1][0] += dp[now][md][1][0] + dp[now][md][0][0];
        }else{
          // not sure aho
          dp[nxt][nmd][1][1] += dp[now][md][1][1];
          dp[nxt][nmd][0][1] += dp[now][md][0][1];
          if(t < d){
            dp[nxt][nmd][1][1] += dp[now][md][1][0];
            dp[nxt][nmd][0][1] += dp[now][md][0][0];
          }else if(t==d){
            dp[nxt][nmd][1][0] += dp[now][md][1][0];
            dp[nxt][nmd][0][0] += dp[now][md][0][0];
          }
        }
        dp[nxt][nmd][0][0] %= MOD;
        dp[nxt][nmd][0][1] %= MOD;
        dp[nxt][nmd][1][0] %= MOD;
        dp[nxt][nmd][1][1] %= MOD;
      }
      dp[now][md][0][0]=dp[now][md][0][1]=dp[now][md][1][0]=dp[now][md][1][1]=0;
    }
    now = 1-now;
    nxt = 1-nxt;
  }
  // last 5 digit
  vvvvvl ddp(2,           // pos
       vvvvl(3,           // mod3
        vvvl(2,           // found 3 flag
         vvl(2,           // less flag
          vl(1e5,0)))));   // last 5 digit
  REP(i,3)REP(j,2)REP(k,2)ddp[now][i][j][k][0] = dp[now][i][j][k];
  REP(_pos,5){
    ll pos = n-5+_pos;
    int d = s[pos]-'0';
    REP(ld5,1e5)REP(md,3){
      // if(_pos==0 && ld5!=0) continue;
      REP(t,10){
        ll nmd = (md*10+t)%3;
        ll nld = (ld5*10+t);
        if(nld>1e5)continue;
        if(t==3){
          // necessarily aho
          ddp[nxt][nmd][1][1][nld] += ddp[now][md][1][1][ld5] + ddp[now][md][0][1][ld5];
          if(t < d) ddp[nxt][nmd][1][1][nld] += ddp[now][md][1][0][ld5] + ddp[now][md][0][0][ld5];
          else if(t == d) ddp[nxt][nmd][1][0][nld] += ddp[now][md][1][0][ld5] + ddp[now][md][0][0][ld5];
        }else{
          // not sure aho
          ddp[nxt][nmd][1][1][nld] += ddp[now][md][1][1][ld5];
          ddp[nxt][nmd][0][1][nld] += ddp[now][md][0][1][ld5];
          if(t < d){
            ddp[nxt][nmd][1][1][nld] += ddp[now][md][1][0][ld5];
            ddp[nxt][nmd][0][1][nld] += ddp[now][md][0][0][ld5];
          }else if(t==d){
            ddp[nxt][nmd][1][0][nld] += ddp[now][md][1][0][ld5];
            ddp[nxt][nmd][0][0][nld] += ddp[now][md][0][0][ld5];
          }
        }
        ddp[nxt][nmd][0][0][nld] %= MOD;
        ddp[nxt][nmd][0][1][nld] %= MOD;
        ddp[nxt][nmd][1][0][nld] %= MOD;
        ddp[nxt][nmd][1][1][nld] %= MOD;
      }
      ddp[now][md][0][0][ld5]=ddp[now][md][0][1][ld5]=ddp[now][md][1][0][ld5]=ddp[now][md][1][1][ld5]=0;
    }
    now = 1-now;
    nxt = 1-nxt;
  }
  ll result = 0;
  REP(is3,2)REP(less,2)REP(md3,3)REP(ckp,1e5){
    if(ckp%p==0)continue;
    if(is3==1 || md3%3==0)result = (result+ddp[now][md3][is3][less][ckp])%MOD;
  }
  return result;
}

int main(){
  string a,b;
  int p;
  cin>>a>>b>>p;
  ll iter = a.size()-1;
  while(iter>=0){
    if(a[iter]!='0'){
      a[iter]--;
      break;
    }
    a[iter]='9';
    iter--;
  }
  ll result = calc(b,p)-calc(a,p);
  result = FIX(result);
  cout << result << endl;
  return 0;
}
0