Submission #6997771


Source Code Expand

import sys
input = sys.stdin.readline
sys.setrecursionlimit(10 ** 7)

from functools import lru_cache

MOD = 10 ** 9 + 7

"""
桁数が異なると矛盾
したがって商が1であることが必要
すると、xがyのsubsetであることが必要
一番上の桁を含むsubsetがx
"""

@lru_cache()
def F_naive(L,R):
    # 
    ret = 0
    ret += F((L+1)//2,(R+1)//2) # 2x,2y
    answer = 0
    for x in range(L,R+1):
        for y in range(x,R+1):
            if x.bit_length() == y.bit_length() and (x^y)+x == y:
                answer += 1
    return answer

@lru_cache()
def F(L,R):
    # x subset y かつ 2x > y
    if L < 0:
        L = 0
    if R < L:
        return 0
    if R == 0:
        return 0
    ret = 0
    # 下一桁が0,0ととる場合
    ret += F((L+1)//2, R//2)
    # 下一桁が0,1ととる場合
    # 2(2x) > (2y+1) iff 2x > y
    ret += F((L+1)//2,(R-1)//2)
    # 下一桁が1,0ととる場合
    # これはsubsetにならない

    # 下一桁が1,1ととる場合
    # 2(2x+1) > 2y+1 iff 2x >= y
    ret += G(L//2,(R-1)//2)
    return ret

@lru_cache()
def G(L,R):
    # x subset y かつ 2x >= y
    if L < 0:
        L = 0
    if R < L:
        return 0
    if R == 0:
        return 1
    ret = 0
    # 下一桁が0,0ととる場合
    ret += G((L+1)//2, R//2)
    # 下一桁が0,1ととる場合
    # 2(2x) >= (2y+1) iff 2x > y
    ret += F((L+1)//2,(R-1)//2)
    # 下一桁が1,0ととる場合
    # これはsubsetにならない

    # 下一桁が1,1ととる場合
    # 2(2x+1) >= 2y+1 iff 2x >= y
    ret += G(L//2,(R-1)//2)
    return ret

L,R = map(int,input().split())
answer = F(L,R) % MOD
print(answer)

Submission Info

Submission Time
Task F - Coincidence
User maspy
Language Python (3.4.3)
Score 600
Code Size 1742 Byte
Status AC
Exec Time 27 ms
Memory 3796 KiB

Judge Result

Set Name Sample All
Score / Max Score 0 / 0 600 / 600
Status
AC × 3
AC × 23
Set Name Test Cases
Sample a01, a02, a03
All a01, a02, a03, b04, b05, b06, b07, b08, b09, b10, b11, b12, b13, b14, b15, b16, b17, b18, b19, b20, b21, b22, b23
Case Name Status Exec Time Memory
a01 AC 22 ms 3572 KiB
a02 AC 23 ms 3572 KiB
a03 AC 26 ms 3680 KiB
b04 AC 22 ms 3572 KiB
b05 AC 23 ms 3572 KiB
b06 AC 26 ms 3680 KiB
b07 AC 22 ms 3572 KiB
b08 AC 27 ms 3796 KiB
b09 AC 27 ms 3680 KiB
b10 AC 25 ms 3680 KiB
b11 AC 25 ms 3680 KiB
b12 AC 24 ms 3680 KiB
b13 AC 25 ms 3680 KiB
b14 AC 25 ms 3680 KiB
b15 AC 24 ms 3680 KiB
b16 AC 23 ms 3680 KiB
b17 AC 26 ms 3680 KiB
b18 AC 23 ms 3572 KiB
b19 AC 23 ms 3572 KiB
b20 AC 24 ms 3572 KiB
b21 AC 25 ms 3680 KiB
b22 AC 25 ms 3680 KiB
b23 AC 26 ms 3680 KiB