結果
| 問題 |
No.1172 Add Recursive Sequence
|
| コンテスト | |
| ユーザー |
Shun_PI
|
| 提出日時 | 2020-08-14 23:00:28 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 446 ms / 4,000 ms |
| コード長 | 1,717 bytes |
| コンパイル時間 | 1,624 ms |
| コンパイル使用メモリ | 175,640 KB |
| 実行使用メモリ | 9,088 KB |
| 最終ジャッジ日時 | 2024-10-10 16:20:26 |
| 合計ジャッジ時間 | 3,639 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 16 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using lint = long long int;
using P = pair<int, int>;
using PL = pair<lint, lint>;
#define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i<i##_end_;i++)
#define IFOR(i, begin, end) for(int i=(end)-1,i##_begin_=(begin);i>=i##_begin_;i--)
#define REP(i, n) FOR(i,0,n)
#define IREP(i, n) IFOR(i,0,n)
#define ALL(a) (a).begin(),(a).end()
constexpr int MOD = 1000000007;
constexpr lint B1 = 1532834020;
constexpr lint M1 = 2147482409;
constexpr lint B2 = 1388622299;
constexpr lint M2 = 2147478017;
constexpr int INF = 2147483647;
void yes(bool expr) {cout << (expr ? "Yes" : "No") << "\n";}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
lint K, N, M;
cin >> K >> N >> M;
vector<lint> a(N);
REP(i, K) cin >> a[i];
vector<lint> c(K);
REP(i, K) cin >> c[i];
FOR(i, K, N) REP(j, K) a[i] = (a[i] + c[j]*a[i-1-j] % MOD) % MOD;
vector<lint> eventl(N);
vector<vector<lint>> eventr(N);
REP(i, M) {
lint l, r;
cin >> l >> r;
r--;
eventl[l]++;
eventr[r].push_back(r-l);
}
vector<lint> cnt(K);
vector<lint> as(K);
REP(i, N) {
REP(j, K) as[j] = (as[j] + cnt[K-1]*a[j] % MOD) % MOD;
lint asnext = 0;
REP(j, K) asnext = (asnext + c[j]*as[K-1-j] % MOD) % MOD;
REP(j, K-1) as[j] = as[j+1];
as[K-1] = asnext;
IREP(j, K-1) cnt[j+1] = cnt[j];
cnt[0] = eventl[i];
lint ans = as[K-1];
REP(j, K) ans = (ans + a[j]*cnt[j] % MOD) % MOD;
cout << ans << "\n";
REP(j, eventr[i].size()) {
if(eventr[i][j] < K) {
cnt[eventr[i][j]]--;
} else {
REP(k, K) {
as[K-1-k] = (as[K-1-k] - a[eventr[i][j]-k] + MOD) % MOD;
}
}
}
}
}
Shun_PI