题意:给定一个大小n<=2e5的滚动数组a,对于贡献值如此定义:k<n&&s<n,从s开始每次向前移动k,移动n次,经过下标的ai的值的累和的最大值
思路:首先应该想到的是a的非约数k都会将a数组遍历一遍,所得到的贡献值是一样的,而a的约数则会选取a数组中n/k个数*k进行贡献值换算,由于s可以任选,我们注意到,n的非约数的遍历实际是约数k遍历不同起点s的1/k的和,因此可以排除非约数k的选取,同时,我们应注意到,约数之间也有这样的关系,对于k1%k2==0 k2与k1的关系实际上等同于非约数与约数的关系
因此,我们只需要考虑n的极大约数即可,即n/它自身的质因子
我们将决策枚举降到了logn级别,但如何处理单点修改中,相隔k的数值累和的最大值问题呢?
我们可以先这样暴力出所有结果,对于每一个极大约数k,我们将从1到k的所有起点的值全部算出来总共nlogn的时间复杂度,注意算的时候遍历下标%记录
用可删除堆就ok了,就是用一个堆记录删除了的值,一个堆记录所有的值,输出的时候如果两个堆顶相等则都pop,否则返回记录ans的堆顶
也可以用multiset(但是常数大一点,而这题卡常数。。。)
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string>
#include<bitset>
#include<cmath>
#include<array>
#include<atomic>
#include<sstream>
#include<stack>
#include<iomanip>
//#include<bits/stdc++.h>
//#define int ll
#define pb push_back
#define endl '\n'
#define x first
#define y second
#define Endl endl
#define pre(i,a,b) for(int i=a;i<=b;i++)
#define rep(i,b,a) for(int i=b;i>=a;i--)
#define si(x) scanf("%d", &x);
#define sl(x) scanf("%lld", &x);
#define ss(x) scanf("%s", x);
#define YES {puts("YES");return;}
#define NO {puts("NO"); return;}
#define all(x) x.begin(),x.end()
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
typedef pair<int, PII> PIII;
typedef pair<char, int> PCI;
typedef pair<int, char> PIC;
typedef pair<double, double> PDD;
typedef pair<ll, ll> PLL;
const int N = 200010, M = 2 * N, B = N, MOD = 998244353;
const int INF = 0x3f3f3f3f;
const ll LLINF = 0x3f3f3f3f3f3f3f3f;
int dx[4] = { -1,0,1,0 }, dy[4] = { 0,1,0,-1 };
int n, m, k;
int a[N];
vector<int> fac;
vector<ll> sum[N];
multiset<ll> heap;
priority_queue<ll> del, val;
ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }
ll lowbit(ll x) { return x & -x; }
ll qmi(ll a, ll b, ll mod) {
ll res = 1;
while (b) {
if (b & 1) res = res * a % mod;
a = a * a % mod;
b >>= 1;
}
return res;
}
inline void init() {}
void getprimes(int n)
{
fac.clear();
pre(i, 2, n)
{
if (n % i == 0)
{
while (n % i == 0)n /= i;
fac.pb(i);
}
}
}
ll get_ans()
{
while (del.size() && val.size() && del.top() == val.top())del.pop(), val.pop();
return val.top();
}
void slove()
{
cin >> n >> m;
getprimes(n);
pre(i, 1, n) cin >> a[i];
heap.clear();
if (fac.empty()) {
ll sum = 0;
pre(i, 1, n) sum += a[i];
cout << sum << endl;
while (m--) {
int p, x;
cin >> p >> x;
cout << sum - a[p] + x << endl;
a[p] = x;
}
return;
}
for (int& p : fac)
{
int k = n / p;
sum[k].resize(k);
pre(i, 0, k - 1)sum[k][i] = 0;
pre(i, 1, n) {
if (i <= k)
sum[k][i % k] = a[i];
else sum[k][i % k] += a[i];
}
//pre(i, 0, k - 1) heap.insert(sum[k][i] * k);
pre(i, 0, k - 1)val.push(sum[k][i] * k);
}
//cout << *heap.rbegin() << endl;
cout << get_ans() << endl;
while (m--)
{
int p, x;
cin >> p >> x;
for (int& fact : fac)
{
int k = n / fact;
//heap.erase(lower_bound(heap.begin(), heap.end(), sum[k][p % k]*k));
del.push(sum[k][p % k] * k);
sum[k][p % k] += x - a[p];
//heap.insert(sum[k][p % k] * k);
val.push(sum[k][p % k] * k);
}
//cout << *heap.rbegin() << endl;
cout << get_ans() << endl;
a[p] = x;
}
while (val.size())val.pop();
while (del.size())del.pop();
return;
}
signed main()
{
int _;
si(_);
//_ = 1;
init();
while (_--)
{
slove();
}
return 0;
}