Codeforces div1 C. Tonya and Burenka-179

题意：给定一个大小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;
}