我乱写了一个算法，太慢，过不了，请问是我的思路有问题还是有优化方法

#include<bits/stdc++.h>
using namespace std;
#define N 200001
#define ll long long
int T,n,q,s,num,a[N];
struct dat{
	ll *a;
	int ys;
	multiset<ll>s;
}d[500];
void get_ys(int x){
	for(int i=1;i<=x;++i){
		if(!(x%i)){
			if(i!=n)d[++num].ys=i;
			if(i*i!=x&&x/i!=n)d[++num].ys=x/i;
		}
	}
}
ll max(ll x,ll y){
	return x>y?x:y;
}
int main(){
	scanf("%d",&T);
	while(T--){
		ll ans=0;
		scanf("%d%d",&n,&q);num=0;
		get_ys(n);
		for(int i=1;i<=n;++i)scanf("%d",a+i);
		for(int i=1;i<=num;++i){
			d[i].a=new ll[d[i].ys];d[i].s.clear();
			for(int j=0;j<d[i].ys;++j){
				d[i].a[j]=0;
				for(int k=j;k<=n;k+=d[i].ys){
					d[i].a[j]+=a[k];
				}
				(d[i].s).insert(d[i].a[j]);
			}
			ans=max(ans,d[i].ys*(*(--(d[i].s).end())));
		}
		printf("%lld\n",ans);
		while(q--){
			int p,x;scanf("%d%d",&x,&p);
			ans=0;
			for(int i=1;i<=num;++i){
				(d[i].s).erase((d[i].s).find(d[i].a[x%d[i].ys]));
				d[i].a[x%d[i].ys]+=p-a[x];
				(d[i].s).insert(d[i].a[x%d[i].ys]);
				ans=max(ans,d[i].ys*(*(--(d[i].s).end())));
			}
			printf("%lld\n",ans);
		}
		for(int i=1;i<=num;++i)delete d[i].a;
	}
	return 0;
}