把过去和现在的笔记记录也搬了过来,也算是给以后留个念想吧,想想一开始打acm就是图一乐,后来发现这游戏还挺上头的,也算是度过了一段电竞生涯(xs)
早些时候的笔记写的好中二,连我自己看着都羞耻。
不过,就喜欢这种羞耻的感觉。
收录的题目大部分是个人认为质量不错的题目,以DP为主,非DP的题目都用※进行了标识。
当然,有些题解思路本身也是源自其他人的,不过除非特殊标注,否则都是用的自己的代码。

题目大意:

HDU 1059 Dividing
题目大意:有六种价值的大理石,每种Ni个(合不超过2e4),问能否分成价值相同的两份。

解:

很容易想到的是二进制优化背包,但是有很多细节,不如是hdu特有的坑吧……
主要的坑点:①只给了32mb内存,导致必须用滚动数组来存,②没有给有多少组样例(深坑)
代码中我保留了之前写二进制的坏习惯,把所有背包平分成了两组,后来发现其实根本没必要,因为二进制分包已经包含了所有可能的分配情况,包括偶数平分成两份的情况。
不过再后来发现,其实还是有必要的,为了简化整个背包的价值范围。
最坏情况的价值范围:全是6块,并且2e4个,这个时候需要2e5的存储范围,每次循环的时候由于滚动数组的原因还得去遍历多次,而且在没有给具体样例数量的情况下很有可能t,因此存储范围的值必须通过题目数据确定,不然很容易被大量小数据的样例卡掉……
最后t了很多发,靠交题二分(?)的方法找到了一个稳定的常数(最大的那个背包重量*4和100的较大值)来作为背包范围,终于过了。
只要没有不给组数的多组样例,你永远可以相信hdu。
不过这题提供了很不错的背包优化思想(不是二分,不是二分),因为求差值的多重背包肯定是存在差值的极限的,而二进制优化背包虽然说优化了物品个数,但是反而拉高了对于差值极限的要求,本题中我之前的代码(分成两部分二进制化)虽然在正常题目中没有必要,但在本题中反而发挥了奇效:仅仅是背包中的物品加倍换来了差值极限的减半,从而导致空间不会爆,对滚动数组的运算量也减少,背包的复杂度两个乘区的相乘,而背包中物品的加倍倍率是比差值极限减小倍率小很多的(表面上是两个乘数一个乘2,一个除2,但因为误差的原因,两个乘数越是接近,乘积就会明显减少……)
总之就是非常神奇,不过主要目的还是过空间复杂度,时间上其实也没有省下太多。

代码

代码
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89

#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <algorithm>
#include <utility>
#include <vector>
#include <istream>
#include <map>
#include <cmath>
#include <stack>
#include <set>
#include <cstring>
#include <string>
#define ll long long
#define maxn 1000005
#define mdl 1000000007
#define clr(a,n) for(int i=0;i<n;i++)a[i]=0
#define cfast std::ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define pll pair<ll,ll>
#define inc(i,a,n) for(ll i=a;i<n;i++)
#define vset(a,n,m) for(ll i=0;i<n;i++)a[i]=m;
#define endl '\n'
#define cst 130000using namespace std;
ll num[6];
ll dp[2][cst*2];
ll val[200];
int main() {
cfast;
int tim = 0;
while (1) {
tim++;
ll sum = 0;
inc(i, 0, 6) {
cin >> num[i];
sum += num[i];
}
if (sum == 0)break;
ll cs =100;
int ts = 0;
inc(i, 0, 6) {
ll tmp = 1;
ll tp = num[i] / 2;
if (num[i] % 2 == 1)val[ts++] = i + 1;
while (tp >= tmp) {
val[ts++] = (i + 1)*tmp;
val[ts++] = (i + 1)*tmp;
tp -= tmp;
tmp *= 2;
cs = max(cs,2*val[ts - 1]);
}
if (tp > 0) {
val[ts++] = (i + 1)*tp;
val[ts++] = (i + 1)*tp;
}
}
//cout << cs << endl; inc(i, 0, 2) {
inc(j, 0, cs * 2) {
dp[i][j] = 0;
}
}
inc(i, 0, ts) {
inc(j, 0, cs * 2)dp[1][j] = 0;
inc(j, 0, cs*2) {
if (i == 0) {
dp[1][cs + val[i]] = 1;
dp[1][cs - val[i]] = 1;
}
else {
if (dp[0][j] == 1) {
if (j + val[i] < cs*2)dp[1][j + val[i]] = 1;
if (j - val[i] >= 0)dp[1][j - val[i]] = 1;
}

}
}
inc(j, 0, cs * 2)dp[0][j] = dp[1][j];

}
cout << "Collection #" << tim << ":" << endl;
if (dp[1][cs] == 1) {
cout << "Can be divided." << endl;
}
else cout << "Can't be divided." << endl;
cout << endl;
}
}