[2017CCPC杭州] HDU 6275 树状数组+二分

给出$n\le 1e6$的两个序列$A,B$。有$m\le 1e5$次操作：1.把$a_x$修改为$y$。2.把$b_x$修改为$y$。3.给出$k$，求最小的$t$，满足$S(t)={\sum }_{i=1}^n\lfloor \frac{t-b_i}{a_i}\rfloor$。其中$a_i\le 1e3$。
直接计算难以处理。$S(t)={\sum }_{i=1}^n(\lfloor \frac{t}{a_i}\rfloor -\lfloor \frac{b_i}{a_i}\rfloor -[t\%a_i-b_i\%a_i<0])$
${\sum }_{i=1}^n\lfloor \frac{b_i}{a_i}\rfloor$可以$O(n)$计算，然后$O(1)$修改。然后要算的就是剩下的部分。注意到$a_i$的值很小，考虑基于值的做法。
统计出$a_i$的个数，这个也同样支持$O(1)$修改。这样只要能快速求出对于相同的$a_i$有多少个$b_i\%a_i$是大于现在$t\%i$的。这个可以用树状数组实现。
由于查询次数不会超过$1e3$次。所以复杂度最高为$1e3\ast 1e3\ast log(S)\ast log(1e3)$。
过的略为勉强。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int N=1e5+7;
ll a[N],b[N];
int cnt[1001];
struct BIT {
	ll t[1001];
	void empty() {
		memset(t,0,sizeof(t));
	}
	int lowbit(int x) {
		return x&(-x);
	}
	void add(int x,ll v) {
		for(int i=x;i<=1000;i+=lowbit(i)) 
			t[i]+=v;
	}
	ll query(int x) {
		ll sum=0;
		for(int i=x;i>0;i-=lowbit(i)) 
			sum+=t[i];
		return sum;
	}
}B[1001];
void empty() {
	for(int i=1;i<=1000;i++) {
		B[i].empty(); 
		cnt[i]=0;
	}
}
ll cal(ll x) {
	ll sum=0;
	for(int i=1;i<=1000;i++) {
		sum+=(x/i)*cnt[i];
		sum-=(cnt[i]-B[i].query(x%i+1));
	}
	return sum; 
}
int main() {
	int T;
	scanf("%d",&T);
	while(T--) {
		empty();
		int n,m;
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++)
			scanf("%lld",&a[i]);
		for(int i=1;i<=n;i++)
			scanf("%lld",&b[i]);
		ll sum=0;
		for(int i=1;i<=n;i++) {
			sum+=b[i]/a[i];
			cnt[a[i]]++;
			B[a[i]].add(b[i]%a[i]+1,1);
		}
		while(m--) {
			int opt;
			scanf("%d",&opt);
			if(opt==1) {
				ll x,y;
				scanf("%lld%lld",&x,&y); // ax->y 
				cnt[a[x]]--;cnt[y]++;
				sum-=b[x]/a[x];sum+=b[x]/y;
				B[a[x]].add(b[x]%a[x]+1,-1);
				B[y].add(b[x]%y+1,1);
				a[x]=y;
			}
			else if(opt==2) {
				ll x,y;
				scanf("%lld%lld",&x,&y); // bx->y 
				sum-=b[x]/a[x];sum+=y/a[x];
				B[a[x]].add(b[x]%a[x]+1,-1);
				B[a[x]].add(y%a[x]+1,1);
				b[x]=y;
			}
			else if(opt==3) {
				ll k;
				scanf("%lld",&k);
				k+=sum;
				ll l=1,r=1e13,ans=0;
				while(l<=r) {
					ll mid=(l+r)>>1;
					if(cal(mid)>=k) ans=mid,r=mid-1;
					else l=mid+1;
				}
				printf("%lld\n",ans);
			}
		}
	} 
	return 0;
}