牛客挑战赛39 - 牛牛的等差数列
牛牛的等差数列
https://ac.nowcoder.com/acm/contest/5157/C
题意
一个数列,两个区间操作:加上一个等差数列;求和,对 取模。
算法(
)
线段树应用,等差数列的处理要注意。
会发现直接储存会爆 long long,而模数又不固定。观察 的范围,会发现
,而
,所以在操作过程中,可以直接对
取模。
代码
#include <cstdio>
#include <algorithm>
using namespace std;
#define ll long long
const int MOD = 3 * 5 * 7 * 11 * 13 * 17 * 19 * 23;
const int INV2 = MOD / 2 + 1;
struct TreeNode {
int val, d, s;
};
TreeNode t[200010 * 4];
int n, q, op, l, r, x, y, a[200010];
void build(int node, int l, int r) {
if (l == r) {
t[node] = {0, 0, a[l]};
return;
}
int mid = (l + r) >> 1;
build(node << 1, l, mid);
build(node << 1 | 1, mid + 1, r);
t[node] = {0, 0, ((ll)t[node << 1].s + t[node << 1 | 1].s) % MOD};
}
void node_upd(int node, ll val, ll d, int n) {
t[node].val = (t[node].val + val) % MOD;
t[node].d = (t[node].d + d) % MOD;
t[node].s = (t[node].s + (val + val + d * (n - 1)) % MOD *
n % MOD * INV2) % MOD;
}
void update(int node, int l, int r, int ul, int ur, ll val, ll d) {
if (r < ul || ur < l) return;
if (ul <= l && r <= ur) {
node_upd(node, (val + d * (l - ul)) % MOD, d, r - l + 1);
return;
}
int mid = (l + r) >> 1;
node_upd(node << 1, t[node].val, t[node].d, mid - l + 1);
node_upd(node << 1 | 1, (t[node].val + (mid - l + 1) * (ll)t[node].d) % MOD,
t[node].d, r - mid);
t[node].val = t[node].d = 0;
update(node << 1, l, mid, ul, ur, val, d);
update(node << 1 | 1, mid + 1, r, ul, ur, val, d);
t[node].s = ((ll)t[node << 1].s + t[node << 1 | 1].s) % MOD;
}
int query(int node, int l, int r, int ql, int qr) {
if (r < ql || qr < l) return 0;
if (ql <= l && r <= qr) return t[node].s;
int mid = (l + r) >> 1;
node_upd(node << 1, t[node].val, t[node].d, mid - l + 1);
node_upd(node << 1 | 1, (t[node].val + (mid - l + 1) * (ll)t[node].d) % MOD,
t[node].d, r - mid);
t[node].val = t[node].d = 0;
return ((ll)query(node << 1, l, mid, ql, qr) +
query(node << 1 | 1, mid + 1, r, ql, qr)) % MOD;
}
int main() {
scanf("%d", &n);
for (int i = 1; i <= n; ++i) scanf("%d", &a[i]);
build(1, 1, n);
scanf("%d", &q);
while (q--) {
scanf("%d%d%d%d", &op, &l, &r, &x);
if (op == 1) {
scanf("%d", &y);
update(1, 1, n, l, r, x, y);
} else printf("%d\n", query(1, 1, n, l, r) % x);
}
return 0;
}
查看2道真题和解析