USACO算法系列三十三——buylow

    题目：http://www.nocow.cn/index.php/Translate:USACO/buylow

    求一个最长的一个递减系列的长度，和序列的种类。

    第一个问题可以用动态规划解决。map[i] = max( map[i+1 ..... n] ) +1; if( map[i] > map[i+1 ... n]) )；

    第二个问题在第一个问题基础之上，通过第一个问题，我们可以得到一个map的数组，然后通过map数组，使用递归我们可以得到所有的系列，但是好像种类数很多，一下子就超出了long int的范围，因为我们必须自己写一个大整数的处理方法。

    代码如下：

#include <iostream> #include <fstream> #include <string.h> #define SIZE 5001 #define LARGER 0X7FFFFFFF #define max(a,b) (a)>(b)?(a):(b) using namespace std; int n; //价格、每个节点以后的数量 long price[SIZE] = {0}; int map[SIZE]; int buytype[SIZE][20] = {0}; bool caltype[SIZE]; bool isuse[SIZE]; ifstream fin("buylow.in"); ofstream fout("buylow.out"); //动态规划计算最长购买系列 int buylow(int index) { int large = -1; for (int i=index+1; i <=n; i ++) { if (price[i] < price[index]) { if (map[i] == -1) { map[i] = buylow(i); } large = max(large, map[i]); } }//end for return (large+1); } void add(int index, int pos) { for (int i=0; i < 19; i ++) { buytype[index][i] += buytype[pos][i]; if (buytype[index][i] > 99999 ) { buytype[index][i+1] ++; buytype[index][i] %= 100000; } } } void print() { fout << map[0] << " "; int i; for (i=19; i >=0; i --) { if (buytype[0][i] > 0) { break; } } fout << buytype[0][i]; for (i=i-1; i >=0; i --) { int bas = 10000; while (bas > 0) { fout << buytype[0][i] / bas; buytype[0][i] = buytype[0][i] % bas; bas /= 10; } } fout << endl; } //深度遍历计算种类数 void buylowtype(int index, int number) { if (number == -1) { caltype[index] = true; buytype[index][0] = 1; } else { for (int i=index+1; i<=n; i ++) { if (map[i]==number && price[i] < price[index] ) { //判断是否重复 bool flag = false; for (int j=i-1; j > index; j --) { if (map[j] == map[i] && price[j] == price[i] && isuse[j] == true) { flag = true; break; }//end if }//end for j if (flag) //重复出现，跳过 { continue; }//end if flag if (!caltype[i]) { buylowtype(i, number-1); isuse[i] = true; } add(index, i); }//end if }//end for } caltype[index] = true; } int main() { //输入数据 fin >> n; price[0] = LARGER; for (int i = 1; i <= n; i ++) { fin >> price[i]; } //初始化数据 memset(map, -1, sizeof(int) * SIZE); //动态规划求最长购买系列 map[0] = buylow(0); buylowtype(0, map[0]-1); print(); return 0; }

    运行结果如下：

Executing... Test 1: TEST OK [0.000 secs, 3456 KB] Test 2: TEST OK [0.000 secs, 3456 KB] Test 3: TEST OK [0.011 secs, 3456 KB] Test 4: TEST OK [0.011 secs, 3456 KB] Test 5: TEST OK [0.011 secs, 3456 KB] Test 6: TEST OK [0.022 secs, 3456 KB] Test 7: TEST OK [0.011 secs, 3456 KB] Test 8: TEST OK [0.000 secs, 3456 KB] Test 9: TEST OK [0.054 secs, 3456 KB] Test 10: TEST OK [0.270 secs, 3456 KB] All tests OK.