結果
| 問題 |
No.868 ハイパー部分和問題
|
| コンテスト | |
| ユーザー |
convexineq
|
| 提出日時 | 2020-11-22 04:41:27 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 4,748 ms / 7,000 ms |
| コード長 | 1,044 bytes |
| コンパイル時間 | 525 ms |
| コンパイル使用メモリ | 82,232 KB |
| 実行使用メモリ | 78,180 KB |
| 最終ジャッジ日時 | 2024-07-23 16:14:50 |
| 合計ジャッジ時間 | 65,791 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 38 |
ソースコード
# coding: utf-8
# Your code here!
import sys
readline = sys.stdin.readline
read = sys.stdin.read
n,k = map(int,readline().split())
*a, = map(int,readline().split())
def mul(c,m,dp): #multiply (1+x^c) up to x^m
for i in range(m,c-1,-1):
dp[i] += dp[i-c]
dp[i] %= MOD
def div(c,m,dp): #divide by (1+x^c) up to x^m
for i in range(m-c+1):
dp[i+c] -= dp[i]
dp[i+c] %= MOD
def mul2(c,m,dp): #multiply (1+x^c) up to x^m
for i in range(m,c-1,-1):
dp[i] += dp[i-c]
dp[i] %= MOD2
def div2(c,m,dp): #divide by (1+x^c) up to x^m
for i in range(m-c+1):
dp[i+c] -= dp[i]
dp[i+c] %= MOD2
MOD = 10**9+7
MOD2 = 10**9+9
dp = [0]*(k+1)
dp[0] = 1
dp2 = [0]*(k+1)
dp2[0] = 1
for i in a:
mul(i,k,dp)
mul2(i,k,dp2)
q = int(readline())
for _ in range(q):
x,v = map(int,readline().split())
if a[x-1]:
div(a[x-1],k,dp)
div2(a[x-1],k,dp2)
if v:
mul(v,k,dp)
mul2(v,k,dp2)
a[x-1] = v
print(1 if dp[-1] or dp2[-1] else 0)
convexineq