結果

問題 No.3207 Digital Font
ユーザー lif4635
提出日時 2025-07-18 23:05:08
言語 PyPy3
(7.3.15)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 5,663 bytes
コンパイル時間 339 ms
コンパイル使用メモリ 82,584 KB
実行使用メモリ 192,520 KB
最終ジャッジ日時 2025-07-19 01:40:18
合計ジャッジ時間 40,051 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 37 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

# input
import sys
input = sys.stdin.readline
II = lambda : int(input())
MI = lambda : map(int, input().split())
LI = lambda : [int(a) for a in input().split()]
SI = lambda : input().rstrip()
LLI = lambda n : [[int(a) for a in input().split()] for _ in range(n)]
LSI = lambda n : [input().rstrip() for _ in range(n)]
MI_1 = lambda : map(lambda x:int(x)-1, input().split())
LI_1 = lambda : [int(a)-1 for a in input().split()]

def graph(n:int, m:int, dir:bool=False, index:int=-1) -> list[set[int]]:
    edge = [set() for i in range(n+1+index)]
    for _ in range(m):
        a,b = map(int, input().split())
        a += index
        b += index
        edge[a].add(b)
        if not dir:
            edge[b].add(a)
    return edge

def graph_w(n:int, m:int, dir:bool=False, index:int=-1) -> list[set[tuple]]:
    edge = [set() for i in range(n+1+index)]
    for _ in range(m):
        a,b,c = map(int, input().split())
        a += index
        b += index
        edge[a].add((b,c))
        if not dir:
            edge[b].add((a,c))
    return edge

inf = 1001001001001001001
ordalp = lambda s : ord(s)-65 if s.isupper() else ord(s)-97
ordallalp = lambda s : ord(s)-39 if s.isupper() else ord(s)-97
yes = lambda : print("Yes")
no = lambda : print("No")
yn = lambda flag : print("Yes" if flag else "No")
def acc(a:list[int]):
    sa = [0]*(len(a)+1)
    for i in range(len(a)):
        sa[i+1] = a[i] + sa[i]
    return sa

prinf = lambda ans : print(ans if ans < 1000001001001001001 else -1)
alplow = "abcdefghijklmnopqrstuvwxyz"
alpup = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
alpall = "abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ"
URDL = {'U':(-1,0), 'R':(0,1), 'D':(1,0), 'L':(0,-1)}
DIR_4 = [[-1,0],[0,1],[1,0],[0,-1]]
DIR_8 = [[-1,0],[-1,1],[0,1],[1,1],[1,0],[1,-1],[0,-1],[-1,-1]]
DIR_BISHOP = [[-1,1],[1,1],[1,-1],[-1,-1]]
prime60 = [2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59]
sys.set_int_max_str_digits(0)
# sys.setrecursionlimit(10**6)
# import pypyjit
# pypyjit.set_param('max_unroll_recursion=-1')

from collections import defaultdict,deque
from heapq import heappop,heappush
from bisect import bisect_left,bisect_right
DD = defaultdict
BSL = bisect_left
BSR = bisect_right

mod = (1 << 61) - 1

class BIT:
    __slots__ = ["n", "data"]
    
    def __init__(self, n):
        self.n = n
        self.data = [0]*(n+1)
    
    def build(self, arr):
        for i,a in enumerate(arr):
            self.data[i+1] = a
        for i in range(1, self.n+1):
            if i + (i&-i) <= self.n:
                self.data[i + (i&-i)] += self.data[i] % mod
    
    def add(self, p, x):
        p += 1
        while p <= self.n:
            self.data[p] += x % mod
            p += p& -p
    
    def sum0(self, r):
        s = 0
        while r:
            s += self.data[r] % mod
            r -= r& -r
        return s % mod
    
    def sum(self, l, r):
        s = 0
        while r:
            s += self.data[r] % mod
            r -= r& -r
        while l:
            s -= self.data[l] % mod
            l -= l& -l
        return s % mod
    
    def get(self, p):
        return self.sum0(p+1) - self.sum0(p)

    def bisect_left(self, w):
        """
        not velify
        sum_{0 <= i < x} >= w となる 最小のx
        """
        if w <= 0: return 0
        x = 0
        r = 1 << self.n.bit_length()
        while r:
            if x + r <= self.n and self.data[x + r] < w:
                w -= self.data[x + r]
                x += r
            r >>= 1
        return x + 1
    
    def bisect_right(self, w):
        """
        sum_{0 <= i < x} <= w となる 最大のx
        """
        assert w >= 0
        x = 0
        r = 1 << self.n.bit_length()
        while r:
            if x + r <= self.n and self.data[x + r] <= w:
                w -= self.data[x + r]
                x += r
            r >>= 1
        return x

    def __str__(self):
        return str([self.get(i) for i in range(self.n)])

h, w = MI()

n = II()

from random import randint

b1 = 37 # 横
b2 = 61 # 縦

num = [randint(1, 36) for i in range(10)]

ib1 = pow(b1, -1, mod)
ib2 = pow(b2, -1, mod) 

hb = [1] * (h + 10)
ihb = [1] * (h + 10)
for i in range(h + 9):
    hb[i+1] = hb[i] * b1 % mod
    ihb[i+1] = ihb[i] * ib1 % mod

wb = [1] * (w + 10)
iwb = [1] * (w + 10)
for i in range(w + 9):
    wb[i+1] = wb[i] * b2 % mod
    iwb[i+1] = iwb[i] * ib2 % mod
    
# hash
# 区間のhashをつくってそれが一致するかどうか？
# 矩形hash？

a = [[] for i in range(h + 1)]
b = [[] for i in range(h + 1)]
for i in range(n):
    i, j, x = MI()
    i -= 1
    j -= 1
    if x == 0: continue
    a[i].append((j, x))
    if x == 6 or x == 9:
        b[h-1-i].append((w-1-j, 15-x))
    else:
        b[h-1-i].append((w-1-j, x))

q = II()
qa = [[] for i in range(h + 1)]
qb = [[] for i in range(h + 1)]
oa = [1] * q
ob = [1] * q
for i in range(q):
    l, d, r, u = MI()
    l -= 1
    d -= 1
    qa[l].append((i, -1, d, u))
    qa[r].append((i, 1, d, u))
    qb[h-r].append((i, -1, w-u, w-d))
    qb[h-l].append((i, 1, w-u, w-d))
    oa[i] = ihb[l] * iwb[d] % mod
    ob[i] = ihb[h-r] * iwb[w-u] % mod


def calc(a, qa, oa):
    bit = BIT(w + 1)
    ans = [0] * q
    for i in range(h + 1):
        for idx, f, d, u in qa[i]:
            p = bit.sum(d, u) # offsetがかかっている
            ans[idx] += f * p
        
        for j, x in a[i]: # この行に追加
            bit.add(j, num[x] * hb[i+1] % mod * wb[j+1] % mod)
        # print(bit)
    for i in range(q):
        ans[i] *= oa[i]
        ans[i] %= mod
    return ans


a1 = calc(a, qa, oa)
a2 = calc(b, qb, ob)

for i in range(q):
    yn(a1[i] % mod == a2[i] % mod)
0