<span>洛谷 题解 P3161 【[CQOI2012]模拟工厂】</span>

本蒟蒻又双叒叕被爆踩辣！

题目链接

Solution:

这题又是一道贪心。。

数据范围：

n<=15 ti<=100,000 gi<=10^9 mi<=10^9

这里就可以看到几个小细节，也可以想出一些思路。

要开long long！

首先，n <= 15， 这个数据可以说是很小辣，它代表有n个任务，又因为

如果接受订单i，则必须恰好在时刻ti交易，不能早也不能晚。同一时刻可以接受多个订单，但每个订单只能被接受一次。

窝们本不知道选择哪些任务才会是利润最高，而且这很明显也无法通过一些规律来知道选择哪个会最优。

所以窝们应该怎么办？？？

枚举！！！

因为只有数据范围只有15，所以窝们可以枚举1到(1 << n) - 1用二进制表示当前这个任务选不选，那我们已经解决了这个任务辣！

接下来就是，然后处理每个任务方案。

窝们针对每个任务方案，依次判断它的可行性。

然后在可行的方案中找出利润最大的（就是保证任务完成的情况下，保证生产力更大。

假设窝们当前有那么对于订单i和j，我们都会得到方程：make为完成订单i时的生产力，time为距离j订单的时间，time_create为用来提升生产力的时间，need是订单j需求量

(make + time_create) × (time − time_create) = need

化简就是：

time_create * time_create + (make - time) * time_create + need - make * time = 0;

只要保证起“根的判别式” >= 0就行。

所以复杂度应该是O（2 的 n次方 * n 的平方）；

具体见代码，这里只有一些思路，个人介意对着代码读，更容易理解

Code:

#include<bits/stdc++.h>

using namespace std;

#define maxn 25
#define maxm
#define ll long long//mi数据1e9，你不开longlong看看 
#define int long long
#define Rep(x, a, b) for(int x = a; x <= b; ++ x)
#define Dep(x, a, b) for(int x = a; x >= b; -- x)
#define Next(x, u) for(int i = head[u]; i ; i = e[i].nxt)

//int read(){窝不喜欢用快读，但缺省源里面放了，所以也没删 
//	int f = 1, x = 0;
//	char c = getchar();
//	while(c < '0' || c > '9'){
//	if(c == '-'){
//	f = -1;
//	}
//	c = getchar();
//	}
//	while(c >= '0' && c <= '9'){
//	x = (x << 1) + (x << 3) + (c ^ 48);
//	c = getchar();
//	}
//	return x * f;
//}

struct node{
int t, g, m;//都如题，t是在什么时刻交易，g是需要多少货，m是可以得到的利润 
}e[maxn], s[maxn];

//e是记录所有的任务
//s是记录当前情况窝们已经选择的任务 

int n, ans;//不多说

bool cmp(node x, node y){
return x.t < y.t;
}//要先按照t排序，要不不好处理

int discriminant(int b, int c){//这个单词是“判别式”，百度说的
if(b * b - 4 * c < 0){//都知道二次函数根的判别式吧（b * b - 4 * a * c） 
return -1;//返回-1，即根的判别式小于0 
}
return floor((-b + sqrt(b * b - 4 * c)) / 2);//否则就返回它的根:((-b + sqrt(b * b - 4 * a * c) / 2 / a) 
}

void dfs(int x){
int now_have = 0, make_force = 1, ans_now = 0, tot = 0;
//now_have是窝们当前拥有的货物
//make_force是当前的生产力
//ans_now是这种方案的利润 
//tot是当前这种选择方案任务的个数 
Rep(i, 1, n){//这里就是寻找当前方案的任务 
if(x & (1 << i - 1)){//用二进制 
s[++ tot] = e[i];//进入s 
ans_now += e[i].m;//就是当前方案的总利润 
}
}
Rep(i, 1, tot){
int t = s[i].t - s[i - 1].t, sum = 0;//sum是所需要的货物 
//t是时间，s[i].t - s[i- 1].t是它的最大值，因为当它前面的货物足够购买这两个货物时，它的t就是此值 
Rep(j, i, tot){//算t的值 
sum += s[j].g;//累加所需要的货物 
if(now_have >= sum){//如果货物还够，那就下一层循环 
continue;
}//否则，就更新t，为t，和它的根的最小值 
t = min(t, discriminant(make_force - s[j].t + s[i - 1].t, sum - now_have - make_force * (s[j].t - s[i - 1].t)));
}
if(t < 0){//如果无法完成次方案 
return;//return 
}
make_force += t;//生产力就要加，（尽可能加就行 
now_have += make_force * (s[i].t - s[i - 1].t - t) - s[i].g;//其余就现在有的产量 
}
ans = max(ans, ans_now);//更新值 
}

signed main(){
scanf("%lld", &n);//输入 
Rep(i, 1, n){//输入每个任务 
scanf("%lld%lld%lld", &e[i].t, &e[i].g, &e[i].m);
}
sort(e + 1, e + n + 1, cmp);//按时刻排序 
Rep(i, 0, (1 << n) - 1){
dfs(i);//每种方案跑一遍 
}
printf("%lld", ans);//输出 
return 0;
}

下面还给两组数据

Inout:
5
5 5 5
7 7 7
9 9 9
11 1 15
16 18 778

Output:
814

Input:
10
1 1 1
2 2 2
3 3 3
4 4 4
5 5 5
6 6 6
7 7 7
8 8 8
9 9 9
10 10 10

Output:
30

Ps:请看懂再抄