CF103Dカワバンを襲撃する時間题解

この質問は線分ツリーを混乱させているようですが、実際に線分ツリーを使用している場合は、役に立たないことがわかります。

したがって、ルート番号アルゴリズムを検討します。

ブロックに分割されていますか？違う！根本的なアルゴリズムもあります-根本的な分割統治。

ルートサインで分割統治するという考えは、しきい値SSを設定することです$S$、ブルートフォース計算がしきい値より大きく、高速計算がしきい値より小さく、S = n S = \ sqrtnを取る$S=\sqrt{n}$ 時間は最適です（それを証明する方法を私に聞かないでください）。

したがって、この質問では、kik_iを目指し$k_{私}$しきい値を手動で500に設定します。500より大きい場合は直接計算され、500以下です。プレフィックスの合計が計算されますが、この質問カードのスペースがあるため、クエリをオフラインにする必要があります。同じkik_iを処理します$k_{私}$ それでおしまい。

コード：

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;
const int MAXN = 320000 + 10;
int n, m;
LL a[MAXN], t[MAXN], k[MAXN], ans[MAXN], f[MAXN];
vector <LL> v[MAXN];

LL read()
{
    
    
	LL sum = 0, fh = 1; char ch = getchar();
	while (ch < '0' || ch > '9') {
    
    if (ch == '-') fh = -1; ch = getchar();}
	while (ch >= '0' && ch <= '9') {
    
    sum = (sum << 3) + (sum << 1) + (ch ^ 48); ch = getchar();}
	return sum * fh;
}

int main()
{
    
    
	n = read();
	for (int i = 1; i <= n; ++i) a[i] = read();
	m = read();
	for (int i = 1; i <= m; ++i) t[i] = read(), k[i] = read();
	for (int i = 1; i <= m; ++i)
	{
    
    
		if (k[i] > 500)
		{
    
    
			for (int j = t[i]; j <= n; j += k[i]) ans[i] += a[j];
			continue;
		}
		v[k[i]].push_back(i);
	}
	for (int i = 1; i <= 500; ++i)
	{
    
    
		if (v[i].empty()) continue;
		memset(f, 0, sizeof(f));
		for (int j = n; j >= 1; --j) f[j] = a[j] + f[j + i];
		for (int j = 0; j < v[i].size(); ++j) ans[v[i][j]] = f[t[v[i][j]]];
	}
	for (int i = 1; i <= m; ++i) printf("%lld\n", ans[i]);
	return 0;
}