HDU 6274 Master of sequence （公式处理 + 二分 + 树状数组）

题意：给定你两个序列a和b，然后有三种不同的操作，分别是修改a数组，修改b数组，给定一个k值，查询一个满足条件的最小的t值，条件是S(t) >= k。

思路：操作3是找一个满足条件的最小的t值，然后很明显t值越大s（t）这个式子也是越大的，满足单调性，很明显我们需要二分这个t来解决。然后我们在观察操作3让我们询问的式子，因为答案是整数，所以这样的分式求和没法直接变成一个式子解决，我们知道分式向下取整，会分为两个部分，一个是bi/ai的整数部分，另一个余数的分数部分即bi%ai，两个这样的分数求差就可以变为整数部分相减之后存下这个值x，然后分数部分的差值如果大于等于0的话，x就不用变，否则的话，x是需要减一的。

我们需要求得S(t) 可以变成 t/ai - bi/ai 整数部分相减得到一个差值，然后对于每个t%ai,有多少个bi%ai比它大那我们就要把前面得到的差值减去这个个数，我们相求出bi/ai的和值，然后枚举t的时候先减去一个这个和，然后就是解决大于这个差值的个数。

题目说ai的值小于1000，也就是说bi%ai的值也最多就只有一千个，所以我们完全可以用数据结构去维护小于当前t%ai的bi%ai有多少个，每次二分t的时候，对于1~1000的范围去查询。时间复杂度就是O(1e31e3logt*log1e3)在时限为10s的情况下是可以跑过去的。

至于check函数，因为我们的树状数组存的值为ai的个数，所以我们需要也用一个数组保存当前的a数组中，每个数的出现次数为多少。然后每次答案累加的就是：

for(int i = 1;i <= 1000;i ++){
	sum += (t/i) * cnt[i] - (cnt[i]-bit[i].query(t%i));
}

t/i*cnt[i]代表的是我当前枚举分母为i时的整数部分答案贡献，分数部分cnt[i]最多减cnt[i]次一，然后我们把小的找出来相减就得了大的次数，然后答案要减去这部分。
考虑如何用数据结构去优化这个查询求和的过程。

代码：

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int MAXN = 1e5 + 7;
const int N = 1e3 + 7;
int a[MAXN],b[MAXN],n,m,cnt[N];

int lowbit(int x){ return x&(-x); }

struct Binary_Index_Tree{
	int c[N];//c[i]维护当前小于t%i的个数是多少
	void add(int p,int val){
		p++;//这个地方要把p先修改一下因为树状数组 维护从1开始 而我们的余数是从0开始的 !!!
		while(p < N){
			c[p] += val;
			p += lowbit(p);
		}
	}
	int query(int p){
		int ans = 0;p++;
		while(p){
			ans += c[p];
			p -= lowbit(p);
		}
		return ans;
	}
}bit[N];

bool check(ll x,ll y){
	ll res = 0;
	for(int i = 1;i <= 1000;i ++){
		res += (x/i)*cnt[i] - (cnt[i]-bit[i].query(x%i));
		if(res >= y + 1000 - i)
			return true;
	}
	return res >= y;
}

int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		memset(cnt,0,sizeof(cnt));
		memset(bit,0,sizeof(bit));
		ll sum = 0;
		scanf("%d%d",&n,&m);
		for(int i = 1;i <= n;i ++) scanf("%d",&a[i]),cnt[a[i]]++;
		for(int i = 1;i <= n;i ++) scanf("%d",&b[i]);
		for(int i = 1;i <= n;i ++){
			sum += b[i]/a[i];
			bit[a[i]].add(b[i]%a[i],1);
		}
		int op,x,y;
		while(m--){
			scanf("%d",&op);
			if(op == 1){
				scanf("%d%d",&x,&y);
				sum += b[x]/y - b[x]/a[x];
				cnt[a[x]]--,cnt[y]++;
				bit[a[x]].add(b[x]%a[x],-1);
				bit[y].add(b[x]%y,1);
				a[x] = y;
			}
			else if(op == 2){
				scanf("%d%d",&x,&y);
				sum += y/a[x] - b[x]/a[x];
				bit[a[x]].add(b[x]%a[x],-1);
				bit[a[x]].add(y%a[x],1);
				b[x] = y;
			}
			else if(op == 3){
				scanf("%d",&x);
				ll ans,xx = x + sum;
				ll l = 0,r = 2e12,mid;
				while(l <= r){
					mid = (l+r)>>1;
					if(check(mid,xx)){
						ans = mid;
						r = mid - 1;
					}
					else
						l = mid + 1;
				}
				printf("%lld\n",ans);
			}
		}
	}
	return 0;
}