結果
| 問題 | 
                            No.685 Logical Operations
                             | 
                    
| コンテスト | |
| ユーザー | 
                             maspy
                         | 
                    
| 提出日時 | 2020-02-26 22:35:53 | 
| 言語 | Python3  (3.13.1 + numpy 2.2.1 + scipy 1.14.1)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 34 ms / 2,000 ms | 
| コード長 | 851 bytes | 
| コンパイル時間 | 104 ms | 
| コンパイル使用メモリ | 12,416 KB | 
| 実行使用メモリ | 10,752 KB | 
| 最終ジャッジ日時 | 2024-10-13 15:32:57 | 
| 合計ジャッジ時間 | 1,791 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge2 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 27 | 
ソースコード
#!/usr/bin/env python3
# %%
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
from functools import lru_cache
# %%
N = int(read())
MOD = 10**9 + 7
# %%
@lru_cache(None)
def solve(N, a, b, c):
    """count x,y such that
    0 <= x + a <= y <= N, x&y + b <= x^y, x^y+c <= x|y"""
    ret = 0
    if N < 2:
        for x in range(N + 1):
            for y in range(N + 1):
                ret += all((x + a <= y, (x & y) + b <= x ^ y, (x ^ y) + c <= x | y))
        return ret
    # x = 2x, y = 2y
    ret += solve(N // 2, a, b, c)
    # x = 2x+1, y = 2y
    ret += solve(N // 2, 1, 0, c)
    # x = 2x, y = 2y+1
    ret += solve((N - 1) // 2, 0, 0, c)
    # x = 2x+1, y = 2y+1
    ret += solve((N - 1) // 2, a, 1, 0)
    return ret % MOD
# %%
answer = solve(N, 0, 1, 1)
print(answer)
            
            
            
        
            
maspy