結果

問題 No.315 世界のなんとか3.5
ユーザー koyumeishi
提出日時 2015-12-08 23:56:45
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 47 ms / 2,000 ms
コード長 2,130 bytes
コンパイル時間 997 ms
コンパイル使用メモリ 74,492 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-14 20:31:30
合計ジャッジ時間 3,018 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:82:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   82 |         scanf("%s %s %d", a_,b_, &p);
      |         ~~~~~^~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff #
プレゼンテーションモードにする

#include <vector>
#include <iostream>
#include <string>
#include <cstdio>
using namespace std;
constexpr long long MOD = 1000000007;
vector<vector<vector<pair<int,int>>>> memo(10, vector<vector<pair<int,int>>>(8));
void make_memo(){
    for(int k=0; k<10; k++){
        for(int a=0; a<2; ++a)
        for(int b=0; b<4; ++b){
            int s = (a<<2)|b;
            vector<int> next(8, 0);
            for(int d=0; d<10; ++d){
                if(a==0 && d>k) break;
                int a_ = a | (d<k);
                int b_ = (b==3||d==3)?3:(b+d)%3;
                int s_ = (a_<<2)|b_;
                next[s_]++;
            }
            for(int i=0; i<8; i++){
                if(next[i]>0) memo[k][s].push_back({i, next[i]});
            }
        }
    }
}
long long calc(string& x, int p){
    vector<vector<long long>> dp(p, vector<long long>(8, 0));
    vector<vector<long long>> dp_(p, vector<long long>(8, 0));
    dp[0][0] = 1;
    int sz = x.size();
    for(int i=0; i<sz; ++i){
        int k = x[i] - '0';
        for(int j=0; j<((i<sz-4)?1:p); ++j){
            for(int a=0; a<2; ++a)
            for(int b=0; b<4; ++b){
                int s=(a<<2)|b;
                long long val = dp[j][s];
                dp[j][s] = 0;
                if(i<sz-5){
                    for(auto& p: memo[k][s]){
                        dp_[0][p.first] += val * p.second;
                        if(dp_[0][p.first] >= MOD){
                            dp_[0][p.first] %= MOD;
                        }
                    }
                }else for(int d=0; d<10; ++d){
                    if(a==0 && d>k) break;
                    int a_ = a | (d<k);
                    int b_ = (b==3||d==3)?3:(b+d)%3;
                    int s_ = (a_<<2)|b_;
                    int j_ = (j*10+d)%p;
                    dp_[j_][s_] += val;
                    if(dp_[j_][s_]>=MOD) dp_[j_][s_]-=MOD;
                }
            }
        }
        swap(dp,dp_);
    }
    long long ret = 0;
    for(int j=1; j<p; ++j)
    for(int a=0; a<2; ++a)
    for(int b=0; b<4; ++b){
        if(b!=0 && b!=3) continue;
        int s=(a<<2)|b;
        ret += dp[j][s];
        if(ret>=MOD) ret-=MOD;
    }
    return ret;
}
int main(){
    char a_[200010];
    char b_[200010];
    int p;
    //cin >> a >> b >> p;
    scanf("%s %s %d", a_,b_, &p);
    string a(a_);
    string b(b_);
    
    make_memo();
    long long ans = calc(b, p);
    ans = (ans - calc(a,p) + MOD) % MOD;
    int x=0;
    int y=0;
    int z=0;
    for(auto& c:a){
        int k = c-'0';
        y |= k==3;
        if(y==0) x = (x+k)%3;
        z = (z*10+k)%p;
    }
    if((x==0 || y==1) && z!=0){
        ans++;
        ans %= MOD;
    }
    printf("%lld\n", ans);
    //cout << ans << endl;
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0