Bin Packing Problem (线段树+multiset）

原题链接

题目：

输入输出样例：

题目大意：

装箱问题。给你含有$n$个数字的一个数组，每个数字代表物品的重量，给你每个箱子的最大容量$c$。有两种装箱子的策略：
1：首次适应算法：每次从前往后遍历每个已经开的箱子，如果可以放得下，那就放进去；如果都放不下，那就新开一个箱子。
2：最佳适应算法：找到所有已经打开的箱子中可以放下这个物品的箱子中剩余容量最小的那个箱子。如果都放不下，那就新开一个箱子。

求：分别运用两种算法，按顺序将所有物品装箱，至少需要几个箱子。

解题思路：

1：首次适应算法：
最先想到的应该是朴素暴力，每次从第一个箱子开始找，如果箱子的剩余容量大于该物品的重量，那就放进去，修改这个箱子的剩余容量，结束寻找；如果找完所有的箱子，都放不进，那就新开一个箱子放这件物品。但是 $10^6$ 的数据太大了，所以得想$O(nlogn)$的办法查找，可以尝试用线段树处理。
设置一个 $cnt$ 变量记录当前开了几个箱子了，初始值为0；线段树中维护的是区间最大值。先开 $n$ 个容量为 $c$ 的箱子，由于需要找最先大于等于 $a_i$ （物体重量）的箱子，那么当该区间的最大值大于等于 $a_i$ 时，我们找该点的左子树比较，直到左子树的区间最大值不满足大于等于 $a_i$ 时，我们再找右子树。当找到位置后，与 $cnt$ 比较，如果小于 $cnt$ ，那么就直接将该位置箱子的容量减去 $a_i$，更新线段树；如果大于 $cnt$ ，代表要重新开一个箱子了，那就把该位置箱子的容量修改为 $c-a_i$，更新线段树。

2：最佳适应算法：
由于不需要记录每个物体放在哪个箱子的序列号，只需要记录最终开的箱子个数，所以，我们想用一个可以自动排序的数据结构，那么查找就方便很多了，那么肯定能想到 $set$ ，但是 $set$ 中不支持重复的元素，所以可以使用它的升级版 $multiset$ 存储已开的箱子的剩余容量，使用 lower_bound() 函数查找可以装入的箱子的位置。如果没找到，那就新开一个箱子，将剩余容量加入到 $multiset$ 中。

代码：

#include <bits/stdc++.h>
#include <iostream>
#include <string.h>
#include <algorithm>
#define qc ios::sync_with_stdio(false); cin.tie(0);cout.tie(0)
#define ll long long
using namespace std;
const int MAXN = 1e6 + 7;
ll tree[MAXN<<2],a[MAXN],b[MAXN],n,c,cnt;
multiset<ll> st;//可以存放有重复的元素 
multiset<ll>::iterator it;

//线段树维护区间最大值 
void pushup(int rt){ tree[rt] = max(tree[rt<<1],tree[rt<<1|1]); }

void build(int rt,int l,int r){
	if(l == r){
		tree[rt] = c;
		return ;
	}
	int mid = (l+r)>>1;
	build(rt<<1,l,mid);build(rt<<1|1,mid+1,r);
	pushup(rt);
}

//将p位置的大小减去val 
void modify(int rt,int l,int r,int p,ll val){
	if(l == r){
		tree[rt] -= val;
		return ;
	}
	int mid = (l+r)>>1;
	if(p <= mid) modify(rt<<1,l,mid,p,val);
	else if(p > mid) modify(rt<<1|1,mid+1,r,p,val);
	pushup(rt);
}

//询问>=val 的第一个箱子的位置 
int query(int rt,int l,int r,ll val){
	if(l == r) return l;
	int mid = (l+r)>>1;
	int res;
	if(tree[rt<<1] >= val)//首次适应算法，第一个>=val的箱子即可装入，所以，尽可能往左子树走，不满足时，再往右子树走 
		res = query(rt<<1,l,mid,val);
	else if(tree[rt<<1|1] >= val)
		res = query(rt<<1|1,mid+1,r,val);
	return res;
}

int main()
{
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%lld%lld",&n,&c);
		for(int i = 1;i <= n;i ++) scanf("%lld",&a[i]);
		build(1,1,n);//假设开了n个箱子，每个箱子的体积为c 
		cnt = 0;//已经开了几个箱子 
		for(int i = 1;i <= n;i ++){
			int p = query(1,1,n,a[i]);
			if(p > cnt) cnt++;//查到的当前合法位置 大于cnt了 证明需要申请新的空间了。
			modify(1,1,n,p,a[i]);
		}
		ll ans = 0;
		st.insert(c - a[1]);
		for (int i = 2; i <= n; ++i) {
			it = st.lower_bound(a[i]);
			if (it == st.end()) {
				st.insert(c - a[i]);//开新的箱子 
				continue;
			}
			int v = *it - a[i];
			st.erase(it);
			st.insert(v);
		}
        printf("%lld %lld\n",cnt,st.size());
        st.clear();
  	}
	return 0;
}

碎碎念：

这个首次适应算法和最佳适应算法就是操作系统里面的分配存储空间等使用算法，既熟悉又陌生的感觉，hhhhhh