結果

問題 No.315 世界のなんとか3.5
ユーザー rickytheta
提出日時 2015-12-08 17:28:13
言語 C++11（廃止可能性あり）
(gcc 13.3.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0