牛牛的凑数游戏

题目

传送门 to nowcoder

思路

为什么我这么弱啊 😢

我想了半个小时才明白，为什么 $a$ 的范围不好用。原来可以直接将 $a$ 排序，然后从小到大取。如果有一个地方出现了断层，也就是 $1+{\sum }_{i=1}^{k-1}{a}_i<a_k$ 时，那么答案就是 $1+{\sum }_{i=1}^{k-1}{a}_i$ 了。

但是这个没法推广到正解！这太臭了吧！因为正解根本不排序！

正解是另一种暴力思路：当前答案（即所谓的 ${\sum }_{i=1}^{k-1}{a}_i$）已经是 $x$ 的时候，我们再把所有不超过 $x+1$ 的数字加入。因为排序只是用来保证，我们只需要做一次。

这玩意儿要做多少次呢？最坏情况下，每次只加入一个上一次的 $x+2$ （因为它在上一次尚未被加入）。那么可以看出，下一次一定大于上一次的两倍。即：操作两次，$x$ 至少翻一倍。而 $x\ge A$ 时就一定会加入所有数字，结束了。所以只需要做 $\mathcal{O}(\log A)$ 次。

我们再借用整体二分的思想。我们进行整体迭代。将所有没有得到答案的询问不停的进行该操作，每次都利用树状数组，$\mathcal{O}[(n+q)\log n]$ 的求出新的答案。显然最多 $\mathcal{O}(\log A)$ 次。足矣。

代码

#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long int_;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline int qkpow(int_ b,int_ q,int Mod){
	int ans = 1;
	for(; q; q>>=1,b=b*b%Mod)
		if(q&1) ans = ans*b%Mod;
	return ans;
}

const int MaxN = 100005;

int_ bit[MaxN], n;
void modify(int id,int v){
	for(int i=id; i<=n; i+=(i&-i))
		bit[i] += v;
}
int_ query(int id){
	int_ res = 0;
	for(int i=id; i; i-=(i&-i))
		res += bit[i];
	return res;
}

int a[MaxN]; // 已经经过 hash
int tmp[MaxN], *End;

struct Query{
	int id, l, r; int_ ans;
	operator int() const {
		return id;
	}
};
Query q[MaxN]; // 所有询问
int_ now[MaxN]; // 新的答案
vector< int > L[MaxN], R[MaxN];
void updata(int len){
	for(int i=1; i<=n; ++i){
		bit[i] = 0;
		L[i].clear();
		R[i].clear();
	}
	for(int i=1; i<=len; ++i){
		L[q[i].l].push_back(i);
		R[q[i].r].push_back(i);
	}
	for(int i=1,x,len; i<=n; ++i){
		len = L[i].size();
		for(int j=0; j<len; ++j){
			x = lower_bound(tmp+1,End,
				q[L[i][j]].ans+2)-tmp-1;
			now[L[i][j]] = -query(x);
		}
		modify(a[i],tmp[a[i]]);
		len = R[i].size();
		for(int j=0; j<len; ++j){
			x = lower_bound(tmp+1,End,
				q[R[i][j]].ans+2)-tmp-1;
			now[R[i][j]] += query(x);
		}
	}
	for(int i=1; i<=len; ++i)
		if(now[i] == q[i].ans) // 结束
			swap(q[i --],q[len --]);
		else q[i].ans = now[i];
	if(len) updata(len); // 继续迭代
}

int main(){
	n = readint();
	int m = readint();
	for(int i=1; i<=n; ++i)
		a[i] = tmp[i] = readint();
	sort(tmp+1,tmp+n+1);
	End = unique(tmp+1,tmp+n+1);
	for(int i=1; i<=n; ++i)
		a[i] = lower_bound(tmp+1,
			End,a[i])-tmp;
	for(int i=1; i<=m; ++i){
		q[i].l = readint();
		q[i].r = readint();
		q[i].ans = 0, q[i].id = i;
	}
	updata(m); sort(q+1,q+m+1);
	for(int i=1; i<=m; ++i)
		printf("%lld\n",q[i].ans+1);
	return 0;
}