結果

問題 No.315 世界のなんとか3.5
ユーザー rpy3cpprpy3cpp
提出日時 2015-12-09 00:43:27
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 958 ms / 2,000 ms
コード長 1,640 bytes
コンパイル時間 162 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 11,648 KB
最終ジャッジ日時 2024-09-14 20:37:06
合計ジャッジ時間 5,959 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,752 KB
testcase_01 AC 30 ms
10,880 KB
testcase_02 AC 194 ms
10,880 KB
testcase_03 AC 32 ms
10,752 KB
testcase_04 AC 33 ms
10,624 KB
testcase_05 AC 57 ms
10,752 KB
testcase_06 AC 70 ms
10,880 KB
testcase_07 AC 58 ms
10,752 KB
testcase_08 AC 54 ms
10,752 KB
testcase_09 AC 34 ms
10,752 KB
testcase_10 AC 32 ms
10,624 KB
testcase_11 AC 31 ms
10,752 KB
testcase_12 AC 63 ms
11,136 KB
testcase_13 AC 63 ms
11,008 KB
testcase_14 AC 107 ms
11,136 KB
testcase_15 AC 108 ms
11,136 KB
testcase_16 AC 108 ms
11,264 KB
testcase_17 AC 110 ms
11,264 KB
testcase_18 AC 62 ms
11,008 KB
testcase_19 AC 62 ms
11,136 KB
testcase_20 AC 62 ms
10,880 KB
testcase_21 AC 62 ms
11,136 KB
testcase_22 AC 108 ms
11,136 KB
testcase_23 AC 108 ms
11,264 KB
testcase_24 AC 108 ms
11,264 KB
testcase_25 AC 108 ms
11,264 KB
testcase_26 AC 107 ms
11,136 KB
testcase_27 AC 108 ms
11,008 KB
testcase_28 AC 173 ms
11,264 KB
testcase_29 AC 260 ms
11,520 KB
testcase_30 AC 258 ms
11,520 KB
testcase_31 AC 262 ms
11,648 KB
testcase_32 AC 183 ms
11,648 KB
testcase_33 AC 93 ms
11,136 KB
testcase_34 AC 336 ms
11,520 KB
testcase_35 AC 958 ms
11,648 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def solve(A, B, P):
    return (f(B, P) - f(A, P) + is_valid(A, P)) % (10**9 + 7)


def is_valid(A, P):
    if len(A) < 6:
        a = int(A)
        return (('3' in A) or (a % 3 == 0)) and (a % P != 0)
    sumA = sum(map(int, A))
    lowerA = int(A[-5:])
    return (('3' in A) or (sumA % 3 == 0)) and (lowerA % P != 0)


def f(A, P):
    if P == 8:
        return fcore(A, P, 4, 8750, 3750)
    elif P == 80:
        return fcore(A, P, 5, 98750, 38742)
    elif P == 800:
        return fcore(A, P, 6, 998750, 353670)

def f_naive(n, P, cum=0):
    cum %= 3
    count = 0
    for i in range(n + 1):
        if i % P == 0:
            continue
        if ('3' in str(i)) or ((i + cum) % 3 == 0):
            count += 1
    return count


def fcore(A, P, d, large, small):
    mod = 10**9 + 7
    inv9 = 111111112
    inv10 = 700000005
    if len(A) <= d:
        return f_naive(int(A), P)
    head = A[:-d]
    tail = A[-d:]
    n = len(head) - 1
    pow10 = pow(10, n, mod)
    pow9 = pow(9, n, mod)
    cum0 = 0
    cum1 = 0
    has_no_3 = True
    for a in head:
        m = ord(a) - 48  # ord('0') = 48
        cum0 = (cum0 + pow10 * m) % mod
        pow10 = (pow10 * inv10) % mod
        if has_no_3:
            if m == 3:
                has_no_3 = False
            m -= m > 3
            cum1 = (cum1 + pow9 * m) % mod
            pow9 = (pow9 * inv9) % mod
    if has_no_3:
        cum = sum(map(int, head))
        return (large * cum0 - small * cum1 + f_naive(int(tail), P, cum)) % mod
    else:
        return (large * cum0 - small * cum1 + int(tail) - int(tail)//P) % mod

A, B, P = input().split()
print(solve(A, B, int(P)))
0