在输入 $r_i$ 的时候应当立即取模，否则更新时

dp[i][j] += (dp[i - 1][j] + dp[i - 1][(j - r[i] + f) % f]) % mod;

$j-r_i + f$ 有可能仍然为负数。