結果

問題 No.315 世界のなんとか3.5
ユーザー rickythetarickytheta
提出日時 2015-12-08 17:28:13
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 347 ms / 2,000 ms
コード長 4,171 bytes
コンパイル時間 2,444 ms
コンパイル使用メモリ 167,696 KB
実行使用メモリ 37,672 KB
最終ジャッジ日時 2023-10-12 22:07:45
合計ジャッジ時間 13,226 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 198 ms
36,900 KB
testcase_01 AC 197 ms
36,964 KB
testcase_02 AC 199 ms
36,952 KB
testcase_03 AC 199 ms
36,768 KB
testcase_04 AC 202 ms
36,952 KB
testcase_05 AC 202 ms
36,840 KB
testcase_06 AC 197 ms
36,832 KB
testcase_07 AC 198 ms
36,960 KB
testcase_08 AC 203 ms
37,116 KB
testcase_09 AC 202 ms
36,828 KB
testcase_10 AC 203 ms
36,888 KB
testcase_11 AC 199 ms
36,960 KB
testcase_12 AC 235 ms
37,216 KB
testcase_13 AC 233 ms
36,960 KB
testcase_14 AC 267 ms
37,052 KB
testcase_15 AC 266 ms
37,144 KB
testcase_16 AC 270 ms
37,148 KB
testcase_17 AC 270 ms
37,084 KB
testcase_18 AC 234 ms
36,892 KB
testcase_19 AC 234 ms
36,884 KB
testcase_20 AC 238 ms
37,216 KB
testcase_21 AC 237 ms
36,956 KB
testcase_22 AC 270 ms
37,088 KB
testcase_23 AC 274 ms
37,100 KB
testcase_24 AC 268 ms
37,148 KB
testcase_25 AC 274 ms
37,100 KB
testcase_26 AC 267 ms
37,088 KB
testcase_27 AC 268 ms
37,220 KB
testcase_28 AC 265 ms
37,124 KB
testcase_29 AC 335 ms
37,532 KB
testcase_30 AC 333 ms
37,460 KB
testcase_31 AC 328 ms
37,388 KB
testcase_32 AC 337 ms
37,384 KB
testcase_33 AC 273 ms
37,188 KB
testcase_34 AC 334 ms
37,672 KB
testcase_35 AC 347 ms
37,536 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