Problem:


Score : 700 points

Problem Statement

We have a rooted binary tree with N vertices, where the vertices are numbered 1 to N.
Vertex 1 is the root, and the parent of Vertex i (i \geq 2) is Vertex \left[ \frac{i}{2} \right].

Each vertex has one item in it. The item in Vertex i has a value of V_i and a weight of W_i.
Now, process the following query Q times:

Given are a vertex v of the tree and a positive integer L.
  Let us choose some (possibly none) of the items in v and the ancestors of v so that their total weight is at most L.
  Find the maximum possible total value of the chosen items.

Here, Vertex u is said to be an ancestor of Vertex v when u is an indirect parent of v, that is, there exists a sequence of vertices w_1,w_2,\ldots,w_k (k\geq 2) where w_1=v, w_k=u, and w_{i+1} is the parent of w_i for each i.

Constraints

All values in input are integers.

1 \leq N < 2^{18}

1 \leq Q \leq 10^5

1 \leq V_i \leq 10^5

1 \leq W_i \leq 10^5

For the values v and L given in each query, 1 \leq v \leq N and 1 \leq L \leq 10^5.

Input

Let v_i and L_i be the values v and L given in the i-th query.
Then, Input is given from Standard Input in the following format:

N
V_1 W_1
:
V_N W_N
Q
v_1 L_1
:
v_Q L_Q

Output

For each integer i from 1 through Q,
the i-th line should contain the response to the i-th query.

Sample Input 1

3
1 2
2 3
3 4
3
1 1
2 5
3 5

Sample Output 1

0
3
3

In the first query, we are given only one choice: the item with (V, W)=(1,2). Since L = 1, we cannot actually choose it, so our response should be 0.

In the second query, we are given two choices: the items with (V, W)=(1,2) and (V, W)=(2,3). Since L = 5, we can choose both of them, so our response should be 3.

Sample Input 2

15
123 119
129 120
132 112
126 109
118 103
115 109
102 100
130 120
105 105
132 115
104 102
107 107
127 116
121 104
121 115
8
8 234
9 244
10 226
11 227
12 240
13 237
14 206
15 227

Sample Output 2

256
255
250
247
255
259
223
253



Sample input:
3
1 2
2 3
3 4
3
1 1
2 5
3 5


Sample output:
0
3
3

Submission status: Runtime Error

Solution code:
import math
from collections import deque
import time

def filter(w1,w2,v):
    if w2<=w1:
        return v
    else:
        return 0

def main():
    N=int(input())
    V=[]
    W=[]
    for i in range(N):
        v,w=map(int,input().split())
        V.append(v)
        W.append(w)
    Q=int(input())
    v=[]
    L=[]
    for i in range(Q):
        v_,l=map(int,input().split())
        v.append(v_)
        L.append(l)

    #start=time.time()

    M=pow(10,5)
    D= int(math.log(N,2))+1#bound of tree
    if D>=2:
        D=D*2//3

    #knapsack DP
    # dp=[[filter(w,W[v],V[v]) for w in range(M+1)] for v in range(N)]
    # for i in range(1, pow(2, D) - 1):
    #     for w in range(1,M+1):
    #         if (w>=W[i]):
    #             dp[i][w]=max(dp[(i-1)//2][w-W[i]]+V[i],dp[(i-1)//2][w])
    #         else:
    #             dp[i][w]=dp[(i-1)//2][w]

    dp = [[filter(w, W[0], V[0]) for w in range(M+1)]]
    for i in range(1, pow(2, D)):
        dp.append([])
        for w in range(1, M + 1):
            if (w >= W[i]):
                dp[i].append(max(dp[(i - 1) // 2][w - W[i]] + V[i], dp[(i - 1) // 2][w]))
            else:
                dp[i].append(dp[(i - 1) // 2][w])
    BoundIndex=pow(2,D)-2

    #start2=time.time()


    for q in range(Q):
        v_q=v[q]
        if v_q-1>BoundIndex:
            d=int(math.log(v_q,2))+1
            k=d-D
            ALL=pow(2,k)
            rest_v=[v_q-1]
            v_=v_q-1
            for _ in range(k-1):
                v_=(v_-1)//2
                rest_v.append(v_)

            parent_node=(v_-1)//2
            res=0
            parent_dp=dp[parent_node]
            for a in range(ALL):
                restW=0
                restV=0
                #bit search
                for flag_index in range(k):
                    if a & 1==1:
                        ind=rest_v[flag_index]
                        restW+=W[ind]
                        restV+=V[ind]
                    a=a>>1
                if L[q]-restW>=0:
                    res=max(res,restV+parent_dp[L[q]-restW])
            print(res)
        else:
            print(dp[v[q]-1][L[q]])

    # exeTime=time.time()-start
    # exeTime2 = time.time() - start2
    #
    # print("time:{}[ms]".format(exeTime*1000))
    # print("bit time:{}[ms]".format(exeTime2*1000))

    # for q in range(Q):
    #     print(dp[v[q]-1][L[q]])





if __name__ == "__main__":
    main()

Submission status: Runtime Error

Solution code:
import math
from collections import deque
import time

def filter(w1,w2,v):
    if w2<=w1:
        return v
    else:
        return 0

def main():
    N=int(input())
    V=[]
    W=[]
    for i in range(N):
        v,w=map(int,input().split())
        V.append(v)
        W.append(w)
    Q=int(input())
    v=[]
    L=[]
    for i in range(Q):
        v_,l=map(int,input().split())
        v.append(v_)
        L.append(l)

    #start=time.time()

    M=pow(10,5)
    D= int(math.log(N,2))+1#bound of tree
    if D>=2:
        D=D*2//3

    #knapsack DP
    # dp=[[filter(w,W[v],V[v]) for w in range(M+1)] for v in range(N)]
    # for i in range(1, pow(2, D) - 1):
    #     for w in range(1,M+1):
    #         if (w>=W[i]):
    #             dp[i][w]=max(dp[(i-1)//2][w-W[i]]+V[i],dp[(i-1)//2][w])
    #         else:
    #             dp[i][w]=dp[(i-1)//2][w]

    dp = [[filter(w, W[0], V[0]) for w in range(M+1)]]
    for i in range(1, pow(2, D) - 1):
        dp.append([])
        for w in range(1, M + 1):
            if (w >= W[i]):
                dp[i].append(max(dp[(i - 1) // 2][w - W[i]] + V[i], dp[(i - 1) // 2][w]))
            else:
                dp[i].append(dp[(i - 1) // 2][w])
    BoundIndex=pow(2,D)-2

    #start2=time.time()


    for q in range(Q):
        v_q=v[q]
        if v_q-1>BoundIndex:
            d=int(math.log(v_q,2))+1
            k=d-D
            ALL=pow(2,k)
            rest_v=[v_q-1]
            v_=v_q-1
            for _ in range(k-1):
                v_=(v_-1)//2
                rest_v.append(v_)

            parent_node=(v_-1)//2
            res=0
            parent_dp=dp[parent_node]
            for a in range(ALL):
                restW=0
                restV=0
                #bit search
                for flag_index in range(k):
                    if a & 1==1:
                        ind=rest_v[flag_index]
                        restW+=W[ind]
                        restV+=V[ind]
                    a=a>>1
                if L[q]-restW>=0:
                    res=max(res,restV+parent_dp[L[q]-restW])
            print(res)
        else:
            print(dp[v[q]-1][L[q]])

    # exeTime=time.time()-start
    # exeTime2 = time.time() - start2
    #
    # print("time:{}[ms]".format(exeTime*1000))
    # print("bit time:{}[ms]".format(exeTime2*1000))

    # for q in range(Q):
    #     print(dp[v[q]-1][L[q]])





if __name__ == "__main__":
    main()

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3
from collections import defaultdict
import sys
input = sys.stdin.buffer.readline
def LI(): return map(int, input().split())
def II(): return int(input())

#solve
def solve():
    n = II()
    v,w = LI()
    nbit = 1 << n.bit_length()
    dp = [[]] * (nbit)
    dp[1] = [(v, w)]
    for i in range(1, nbit >> 1):
        dpi = dp[i]
        if 2 * i + 1 <= n:
            v0, w0 = LI()
            v1, w1 = LI()
            dp21 = dpi[:]
            dp20 = dpi[:]
            dp21.append((v1, w1))
            dp20.append((v0, w0))
            dp[2 * i + 1] = dp21
            dp[2 * i] = dp20
        elif 2 * i <= n:
            v0, w0 = LI()
            dp20 = dpi[:]
            dp20.appned((v0, w0))
            dp[2 * i] = dp20
        else:
            break
    return
    M = 10 ** 5 + 1
    M0 = [0] * M
    d = [M0]
    for i in range(1, 1 << n.bit_length() // 2):
        d1 = defaultdict(int)
        d1[0] = 0
        di = M0[:]
        for v, w in dp[i]:
            for key, value in list(d1.items()):
                keyw = key + w
                valuev = value + v
                if keyw >= M: continue
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    di[keyw] = valuev
        for j in range(M - 1):
            if di[j + 1] < di[j]:
                di[j + 1] = di[j]
        d.append(di)
    for _ in range(II()):
        V, L = LI()
        t = dp[V]
        lent = V.bit_length()
        n1 = lent // 2
        n2 = lent - n1
        vw2 = t[n1:]
        V = V >> n2
        lw = [0] * (1 << n2)
        lv = [0] * (1 << n2)
        dv = d[V]
        ans = dv[L]
        for mask in range(1, 1 << n2):
            a = mask & -mask
            b = a.bit_length() - 1
            v, w = vw2[b]
            ww = lw[mask - a] + w
            lw[mask] = ww
            if ww > L: continue
            vv = lv[mask - a] + v
            lv[mask] = vv
            if dv[L - ww] + vv > ans:
                ans = dv[L - ww] + vv
        print(ans)
    return


solve()

Submission status: Runtime Error

Solution code:
import sys
import numpy as np
from numba import njit


@njit('i8[:](i8,i8[:],i8[:],i8,i8[:],i8[:])')
def solve(n, vvv, www, q, uuu, lll):
    weight_limit = 10 ** 5
    precalc_limit = min(1 << 9, n + 1)
    precalc = [np.zeros(weight_limit + 1, dtype=np.int64)]

    for u in range(1, precalc_limit):
        v = vvv[u]
        w = www[u]
        parent_dp = precalc[u >> 1]
        dp = parent_dp.copy()
        dp[w:] = np.maximum(parent_dp[w:], parent_dp[:-w] + v)
        precalc.append(dp)

    buf = np.zeros(q, dtype=np.int64)
    for i in range(q):
        u = uuu[i]
        l = lll[i]
        if u < precalc_limit:
            buf[i] = precalc[u][l]
            continue
        dp = {0: 0}
        while u >= precalc_limit:
            v = vvv[u]
            w = www[u]
            for cw, cv in list(dp.items()):
                nw = cw + w
                nv = cv + v
                if nw > weight_limit:
                    continue
                if nw not in dp or dp[nw] < nv:
                    dp[nw] = nv
            u >>= 1
        ans = 0
        for w, v in dp.items():
            ans = max(ans, v + precalc[u][l - w])
        buf[i] = ans

    return buf


n, *inp = map(int, sys.stdin.buffer.read().split())
vvv = np.array([0] + inp[0:n * 2:2], dtype=np.int64)
www = np.array([0] + inp[1:n * 2:2], dtype=np.int64)
q = inp[n * 2]
uuu = np.array(inp[n * 2 + 1::2], dtype=np.int64)
lll = np.array(inp[n * 2 + 2::2], dtype=np.int64)
ans = solve(n, vvv, www, q, uuu, lll)
print('\n'.join(map(str, ans)))

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3
from collections import defaultdict
import sys
input = sys.stdin.buffer.readline
def LI(): return map(int, input().split())
def II(): return int(input())

#solve
def solve():
    n = II()
    v,w = LI()
    nbit = 1 << n.bit_length()
    dp = [[]] * (nbit)
    dp[1] = [(v, w)]
    for i in range(1, nbit >> 1):
        dpi = dp[i]
        if 2 * i + 1 <= n:
            v0, w0 = LI()
            v1, w1 = LI()
            dp21 = dpi[:]
            dp20 = dpi[:]
            dp21.append((v1, w1))
            dp20.append((v0, w0))
            dp[2 * i + 1] = dp21
            dp[2 * i] = dp20
        elif 2 * i <= n:
            v0, w0 = LI()
            dp20 = dpi[:]
            dp20.appned((v0, w0))
            dp[2 * i] = dp20
        else:
            break
    M = 10 ** 5 + 1
    M0 = [0] * M
    d = [M0]
    for i in range(1, 1 << n.bit_length() // 2):
        d1 = defaultdict(int)
        d1[0] = 0
        di = M0[:]
        for v, w in dp[i]:
            for key, value in list(d1.items()):
                keyw = key + w
                valuev = value + v
                if keyw >= M: continue
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    di[keyw] = valuev
        for j in range(M - 1):
            if di[j + 1] < di[j]:
                di[j + 1] = di[j]
        d.append(di)
    for _ in range(II()):
        V, L = LI()
        t = dp[V]
        lent = V.bit_length()
        n1 = lent // 2
        n2 = lent - n1
        vw2 = t[n1:]
        V = V >> n2
        lw = [0] * (1 << n2)
        lv = [0] * (1 << n2)
        dv = d[V]
        ans = dv[L]
        for mask in range(1, 1 << n2):
            a = mask & -mask
            b = a.bit_length() - 1
            v, w = vw2[b]
            ww = lw[mask - a] + w
            lw[mask] = ww
            if ww > L: continue
            vv = lv[mask - a] + v
            lv[mask] = vv
            if dv[L - ww] + vv > ans:
                ans = dv[L - ww] + vv
        print(ans)
    return


solve()

Submission status: Runtime Error

Solution code:
N=int(input())
VW=[[0,0]]+[list(map(int,input().split())) for _ in range(N)]
Q=int(input())
QU=[list(map(int,input().split())) for _ in range(Q)]
W=10**5

from collections import deque
dp=[[0]*(W+1) for _ in range(2**9)]
q=deque([1])
#深さ優先探索的に調べる
while q:
  a=q.pop()
  for i in range(W+1):
    if i>=VW[a][1]:
      dp[a][i]=max([dp[a//2][i],dp[a//2][i-VW[a][1]]+VW[a][0]])
    else:
      dp[a][i]=dp[a//2][i]
  if a<2**8:
    if 2*a<=N:
      q.append(2*a)
    if 2*a+1<=N:
      q.append(2*a+1)

for i in range(Q):
  if QU[i][0]<2**9:
    print(dp[QU[i][0]][QU[i][1]])
  else:
    b=QU[i][0]
    ls=[b]
    while b>=2**9:
      b=b//2
      if b>=2**9:
        ls.insert(0,b)  
    A=[]  
    for k in range(2**(len(ls))):
      v=0
      w=0
      for j in range(len(ls)):
        if (k>>j)&1==1:
          v+=VW[ls[j]][0]
          w+=VW[ls[j]][1]
      A.append([v,w])
    B=[]  
    for j in len(A):
      B.append(VW[b][QU[j][1]-w]+QU[j][0])
    print(max(B))

Submission status: Runtime Error

Solution code:
N=int(input())
VW=[[0,0]]+[list(map(int,input().split())) for _ in range(N)]
Q=int(input())
QU=[list(map(int,input().split())) for _ in range(Q)]
W=10**5

from collections import deque
dp=[[0]*(W+1) for _ in range(2**9)]
q=deque([1])
#深さ優先探索的に調べる
while q:
  a=q.pop()
  for i in range(W+1):
    if i>=VW[a][1]:
      dp[a][i]=max([dp[a//2][i],dp[a//2][i-VW[a][1]]+VW[a][0]])
    else:
      dp[a][i]=dp[a//2][i]
  if a<2**8:
    if 2*a<=N:
      q.append(2*a)
    if 2*a+1<=N:
      q.append(2*a+1)

for i in range(Q):
  if QU[i][0]<2**9:
    print(dp[QU[i][0]][QU[i][1]])
  else:
    b=QU[i][0]
    ls=[b]
    while b>=2**9:
      b=b//2
      ls.insert(0,b)  
    A=[]  
    for k in range(2**len(ls)):
      v=0
      w=0
      for j in range(len(ls)):
        if (k>>j)&1==1:
          v+=VW[ls[j]][0]
          w+=VW[ls[j]][1]
      A.append([v,w])
    B=[]  
    for j in len(A):
      B.append(VW[b//2][QU[j][1]-w]+QU[j][0])
    print(max(B))

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3
from collections import defaultdict
import sys
input = sys.stdin.buffer.readline
def LI(): return map(int, input().split())
def II(): return int(input())

#solve
def solve():
    n = II()
    v,w = LI()
    nbit = 1 << n.bit_length()
    dp = [[]] * (nbit)
    dp[1] = [(v, w)]
    for i in range(1, nbit >> 1):
        dpi = dp[i]
        if 2 * i + 1 <= n:
            v0, w0 = LI()
            v1, w1 = LI()
            dp21 = dpi[:]
            dp20 = dpi[:]
            dp21.append((v1, w1))
            dp20.append((v0, w0))
            dp[2 * i + 1] = dp21
            dp[2 * i] = dp20
        elif 2 * i <= n:
            v0, w0 = LI()
            dp20 = dpi[:]
            dp20.appned((v0, w0))
            dp[2 * i] = dp20
        else:
            break
    M = 10 ** 5 + 1
    M0 = [0] * M
    d = [M0]
    for i in range(1, 1 << n.bit_length() // 2):
        d1 = defaultdict(int)
        d1[0] = 0
        di = M0[:]
        for v, w in dp[i]:
            for key, value in list(d1.items()):
                keyw = key + w
                valuev = value + v
                if keyw >= M: continue
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    di[keyw] = valuev
        for j in range(M - 1):
            if di[j + 1] < di[j]:
                di[j + 1] = di[j]
        d.append(di)
    for _ in range(II()):
        V, L = LI()
        t = dp[V]
        lent = V.bit_length()
        n1 = lent // 2
        n2 = lent - n1
        vw2 = t[n1:]
        V = V >> n2
        d1 = defaultdict(int)
        d1[0] = 0
        dv = d[V]
        ans = dv[L]
        for v, w in vw2:
            for key, value in list(d1.items()):
                keyw = key + w
                if keyw > L: continue
                valuev = value + v
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    if dv[L - keyw] + valuev > ans:
                        ans = dv[L - keyw] + valuev
        print(ans)
    return


solve()

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3
from collections import defaultdict
import sys
input = sys.stdin.buffer.readline
def LI(): return map(int, input().split())
def II(): return int(input())

#solve
def solve():
    n = II()
    v,w = LI()
    nbit = 1 << n.bit_length()
    dp = [[]] * (nbit)
    dp[1] = [(v, w)]
    for i in range(1, nbit >> 1):
        dpi = dp[i]
        if 2 * i + 1 <= n:
            v0, w0 = LI()
            v1, w1 = LI()
            dp21 = dpi[:]
            dp20 = dpi[:]
            dp21.append((v1, w1))
            dp20.append((v0, w0))
            dp[2 * i + 1] = dp21
            dp[2 * i] = dp20
        elif 2 * i <= n:
            v0, w0 = LI()
            dp20 = dpi[:]
            dp20.appned((v0, w0))
            dp[2 * i] = dp20
        else:
            break
    M = 10 ** 5 + 1
    d = [[]]
    M0 = [0] * M
    for i in range(1, 1 << n.bit_length() // 2):
        d1 = defaultdict(int)
        d1[0] = 0
        di = M0[:]
        for v, w in dp[i]:
            for key, value in list(d1.items()):
                keyw = key + w
                valuev = value + v
                if keyw >= M: continue
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    di[keyw] = valuev
        for j in range(M - 1):
            if di[j + 1] < di[j]:
                di[j + 1] = di[j]
        d.append(di)
    for _ in range(II()):
        V,L = LI()
        t = dp[V]
        lent = V.bit_length()
        n1 = lent // 2
        n2 = lent - n1
        vw2 = t[n1:]
        V = V >> n2
        d1 = defaultdict(int)
        d1[0] = 0
        dv = d[V]
        ans = dv[L]
        for v, w in vw2:
            for key, value in list(d1.items()):
                keyw = key + w
                if keyw > L: continue
                valuev = value + v
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    if dv[L - keyw] + valuev > ans:
                        ans = dv[L - keyw] + valuev
        print(ans)
    return


solve()

Submission status: Runtime Error

Solution code:
import math
from collections import deque
import time

def filter(w1,w2,v):
    if w2<=w1:
        return v
    else:
        return 0

def main():
    N=int(input())
    V=[]
    W=[]
    for i in range(N):
        v,w=map(int,input().split())
        V.append(v)
        W.append(w)
    Q=int(input())
    v=[]
    L=[]
    for i in range(Q):
        v_,l=map(int,input().split())
        v.append(v_)
        L.append(l)

    #start=time.time()

    M=pow(10,5)
    D= int(math.sqrt(N+1))#bound of tree
    if D>=2:
        D=D*2//3

    #knapsack DP
    dp=[[filter(w,W[v],V[v]) for w in range(M+1)] for v in range(pow(2,D))]
    for i in range(1, pow(2, D) - 1):
        for w in range(1,M+1):
            if (w>=W[i]):
                dp[i][w]=max(dp[(i-1)//2][w-W[i]]+V[i],dp[(i-1)//2][w])
            else:
                dp[i][w]=dp[(i-1)//2][w]
    BoundIndex=pow(2,D)-2

    #start2=time.time()


    for q in range(Q):
        v_q=v[q]
        if v_q-1>BoundIndex:
            d=int(math.sqrt(v_q))+1
            k=d-D
            ALL=pow(2,k)
            rest_v=[v_q-1]
            v_=v_q-1
            for _ in range(k-1):
                v_=(v_-1)//2
                rest_v.append(v_)

            parent_node=(v_-1)//2
            res=0
            parent_dp=dp[parent_node]
            for a in range(ALL):
                restW=0
                restV=0
                #bit search
                for flag_index in range(k):
                    if a & 1==1:
                        ind=rest_v[flag_index]
                        restW+=W[ind]
                        restV+=V[ind]
                    a=a>>1
                if L[q]-restW>=0:
                    res=max(res,restV+parent_dp[L[q]-restW])
            print(res)
        else:
            print(dp[v[q]-1][L[q]])

    # exeTime=time.time()-start
    # exeTime2 = time.time() - start2
    # 
    # print("time:{}[ms]".format(exeTime*1000))
    # print("bit time:{}[ms]".format(exeTime2*1000))

    # for q in range(Q):
    #     print(dp[v[q]-1][L[q]])





if __name__ == "__main__":
    main()

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3
from collections import defaultdict
import sys, itertools
input = sys.stdin.buffer.read
def LI(): return list(map(int, input().split()))
def II(): return int(input())
def LIR(n): return [LI() for _ in range(n)]

#solve
def solve():
    n = II()
    vw = LIR(n)
    nbit = 1 << n.bit_length()
    b = [i.bit_length() for i in range(nbit)]
    dp = [[None] * b[i] for i in range(nbit)]
    dp[1] = [(vw[0][0], vw[0][1])]
    for i in range(1 , nbit >> 1):
        dpi = dp[i]
        if 2 * i + 1 <= n:
            v1, w1 = vw[2 * i]
            v0, w0 = vw[2 * i - 1]
            dp21 = dp[2 * i + 1]
            dp20 = dp[2 * i]
            for i, j in enumerate(dpi):
                dp21[i] = j
                dp20[i] = j
            dp21[i + 1] = (v1, w1)
            dp20[i + 1] = (v0, w0)
        elif 2 * i <= n:
            v0, w0 = vw[2 * i - 1]
            dp20 = dp[2 * i]
            for i, j in enumerate(dpi):
                dp20[i] = j
            dp20[i + 1] = (v0, w0)
    M = 10 ** 5 + 1
    d = [[0] * M for _ in range(1 << n.bit_length() // 2)]
    for i in range(1, 1 << n.bit_length() // 2):
        d1 = defaultdict(int)
        d1[0] = 0
        di = d[i]
        for v, w in dp[i]:
            for key, value in list(d1.items()):
                keyw = key + w
                valuev = value + v
                if keyw >= M: continue
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    di[keyw] = valuev
        for j in range(M - 1):
            if di[j + 1] < di[j]:
                di[j + 1] = di[j]
    for _ in range(II()):
        V, L = LI()
        t = dp[V]
        lent = V.bit_length()
        n1 = lent // 2
        n2 = lent - n1
        vw2 = t[n1:]
        V = V >> n2
        d1 = defaultdict(int)
        d1[0] = 0
        dv = d[V]
        ans = dv[L]
        for v, w in vw2:
            for key, value in list(d1.items()):
                keyw = key + w
                if keyw > L: continue
                valuev = value + v
                if d1[keyw] < valuev:
                    d1[keyw] = valuev
                    if dv[L - keyw] + valuev > ans:
                        ans = dv[L - keyw] + valuev
        print(ans)
    return


#main
if __name__ == '__main__':
    solve()

Submission status: Runtime Error

Solution code:
# coding: utf-8
import sys
import numpy as np
sys.setrecursionlimit(10 ** 7)

sr = lambda: sys.stdin.readline().rstrip()
ir = lambda: int(sr())
lr = lambda: list(map(int, sr().split()))

N = ir()
VW = [lr() for _ in range(N)]
X = [None] * (N+1)

def dfs(se, i):
    v, w = VW[i-1]
    se2 = set()
    for value, weight in se:
        x1 = value + v
        x2 = weight + w
        if x2 <= 10 ** 5 + 1:
            se2.add((x1, x2))
    se |= se2
    X[i] = se
    if i * 2 <= N // 2:
        dfs(se, i*2)
    if i * 2 + 1 <= N // 2:
        dfs(se, i*2+1)

dfs({(0, 0)}, 1)
Q = ir()
Z = []

def dfs_rev(se, i):
    assert i > N // 2
    v, w = VW[i-1]
    se2 = set()
    for value, weight in se:
        x1 = value + v
        x2 = weight + w
        if x2 <= 10 ** 5 + 1:
            se2.add((x1, x2))
    se |= se2
    if i // 2 > N // 2:
        return dfs(se, i // 2)
    else:
        return (se, i // 2)

for _ in range(Q):
    v, L = lr()
    if v <= N // 2:
        answer = 0
        for value, weight in X[v]:
            if weight <= L and value > answer:
                answer = value
        print(answer)
    else:
        down, parent = dfs_rev({(0, 0)}, v)
        up = X[parent]
        up_list = [0] * (L+1)
        for value, weight in up:
            if weight <= L and up_list[weight] < value:
                up_list[weight] = value
        down_list = [0] * (L+1)
        for value, weight in down:
            if weight <= L and down_list[weight] < value:
                down_list[weight] = value
        up_list = np.array(up_list); up_list = np.maximum.accumulate(up_list)
        down_list = np.array(down_list); down_list = np.maximum.accumulate(down_list)
        answer = np.max(up_list + down_list[::-1])
        print(answer)
        

#print('\n'.join(map(str, Z)))

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3

import sys
import math
from bisect import bisect_right as br
from bisect import bisect_left as bl
sys.setrecursionlimit(2147483647)
from heapq import heappush, heappop,heappushpop
from collections import defaultdict
from itertools import accumulate
from collections import Counter
from collections import deque
from operator import itemgetter
from itertools import permutations
mod = 10**9 + 7
inf = float('inf')
def I(): return int(sys.stdin.readline())
def LI(): return list(map(int,sys.stdin.readline().split()))

n = I()
V = []
W = []
for _ in range(n):
    v, w = LI()
    V.append(v)
    W.append(w)

q = I()
vl = []
maxL = 0
for _ in range(q):
    v, l = LI()
    maxL = max(maxL, l)
    vl.append([v,l])

m = min(n, 1 << 10)
dp = [[0]*(maxL + 1) for _ in range(m + 1)]

for i in range(1, m+1):
    v, w = V[i-1], W[i-1]
    for j in range(maxL+1):
        if j < w:
            dp[i][j] = dp[i>>1][j]
        else:
            dp[i][j] = max(dp[i>>1][j], dp[i>>1][j - w] + v)


for i in range(q):
    v, L = vl[i]
    if v <= m:
        print(dp[v][L])
        continue
    cnt = 0
    tmp = v
    v_ = []
    w_ = []
    while tmp > m:
        v_.append(v[tmp-1])
        w_.append(w[tmp-1])
        tmp >>= 1
        cnt += 1
    ret = dp[tmp][L]
    for j in range(1, 1 << cnt):
        weight_sum = 0
        value_sum = 0
        for v, w in zip(v_, w_):
            if j & 1:
                weight_sum += w
                value_sum += v
                if weight_sum > L:
                    break
            j >>= 1
        else:
            value_sum += dp[tmp][L - weight_sum]
            ret = max(ret, value_sum)
    print(ret)

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3

import sys
import math
from bisect import bisect_right as br
from bisect import bisect_left as bl
sys.setrecursionlimit(2147483647)
from heapq import heappush, heappop,heappushpop
from collections import defaultdict
from itertools import accumulate
from collections import Counter
from collections import deque
from operator import itemgetter
from itertools import permutations
mod = 10**9 + 7
inf = float('inf')
def I(): return int(sys.stdin.readline())
def LI(): return list(map(int,sys.stdin.readline().split()))

n = I()
V = []
W = []
for _ in range(n):
    v, w = LI()
    V.append(v)
    W.append(w)

q = I()
vl = []
maxL = 0
for _ in range(q):
    v, l = LI()
    maxL = max(maxL, l)
    vl.append([v,l])

m = min(n, 1 << 10)
dp = [[0]*(maxL + 1) for _ in range(m + 1)]

for i in range(1, m+1):
    v, w = V[i-1], W[i-1]
    for j in range(maxL+1):
        if j < w:
            dp[i][j] = dp[i>>1][j]
        else:
            dp[i][j] = max(dp[i>>1][j], dp[i>>1][j - w] + v)


for i in range(q):
    v, L = vl[i]
    if v <= m:
        print(dp[v][L])
        continue
    cnt = 0
    tmp = v
    v_ = []
    w_ = []
    while tmp > m:
        v_.append(v[tmp-1])
        w_.append(w[tmp-1])
        tmp >>= 1
        cnt += 1
    ret = 0
    for j in range(1 << cnt):
        weight_sum = 0
        value_sum = 0
        idx = j
        for v, w in zip(v_, w_):
            if idx & 1:
                weight_sum += w
                value_sum += v
                if weight_sum > L:
                    break
            idx >>= 1
        else:
            value_sum += dp[tmp][L - weight_sum]
            ret = max(ret, value_sum)
    print(ret)

Submission status: Runtime Error

Solution code:
#!/usr/bin/env python3

import sys
import math
from bisect import bisect_right as br
from bisect import bisect_left as bl
sys.setrecursionlimit(2147483647)
from heapq import heappush, heappop,heappushpop
from collections import defaultdict
from itertools import accumulate
from collections import Counter
from collections import deque
from operator import itemgetter
from itertools import permutations
mod = 10**9 + 7
inf = float('inf')
def I(): return int(sys.stdin.readline())
def LI(): return list(map(int,sys.stdin.readline().split()))

n = I()
V = []
W = []
for _ in range(n):
    v, w = LI()
    V.append(v)
    W.append(w)

q = I()
vl = []
maxL = 0
for _ in range(q):
    v, l = LI()
    maxL = max(maxL, l)
    vl.append([v,l])

m = min(n, 1 << 10)
dp = [[0]*(maxL + 1) for _ in range(m + 1)]

for i in range(1, m+1):
    v, w = V[i-1], W[i-1]
    for j in range(maxL+1):
        if j < w:
            dp[i][j] = dp[i>>1][j]
        else:
            dp[i][j] = max(dp[i>>1][j], dp[i>>1][j - w] + v)


for i in range(q):
    v, L = vl[i]
    if v <= m:
        print(dp[v][L])
        continue
    cnt = 0
    tmp = v
    v_ = []
    w_ = []
    while tmp > m:
        v_.append(v[tmp-1])
        w_.append(w[tmp-1])
        tmp >>= 1
        cnt += 1
    ret = 0
    for j in range(1 << cnt):
        weight_sum = 0
        value_sum = 0
        for v, w in zip(v_, w_):
            if j & 1:
                weight_sum += w
                value_sum += v
                if weight_sum > L:
                    break
            j >>= 1
        else:
            value_sum += dp[tmp][L - weight_sum]
            ret = max(ret, value_sum)
    print(ret)

Submission status: Runtime Error

Solution code:
N = int(input())
V = []
W = []

for _ in range(N):
    v, w = map(int,input().split())
    V.append(v)
    W.append(w)

L = [[-1 for j in range(10 ** 5 + 1)] for j in range(2 ** 9)]
L[0][0] = 0

for i in range(min(N, 2 ** 9 - 1)):
    for j in range(10 ** 5 + 1):
        L[i+1][j] = max(L[i+1][j], L[(i+1)//2][j])
        if j + W[i] <= 10 ** 5:
            L[i+1][j+W[i]] = max(L[(i+1)//2][j] + V[i], L[i+1][j+W[i]])
        if j != 0:
            L[i+1][j] = max(L[i+1][j], L[i+1][j-1])

Q = int(input())
for _ in range(Q):
    v, l = map(int,input().split())
    if v < 2 ** 9:
        print(L[v][l])
    else:
        ans = 0
        K = []
        while v >= 2 ** 9:
            K.append(v-1)
            v //= 2
        for i in range(len(K) ** 2):
            tmp_v = 0
            tmp_w = 0
            for j in range(K):
                if i & 1:
                    tmp_v += V[K[j]]
                    tmp_w += W[K[j]]
                i >>= 1
            if l - tmp_w >= 0:
                ans = max(ans, L[v][l-tmp_w] + tmp_v)
        print(ans)