TopCoder SRM 667 Div.2题解
概览：
T1 枚举 T2 状压DP T3 DP
TopCoder SRM 667 Div.2 T1
解题思路
由于数据范围很小，所以直接枚举所有点，判断是否可行。时间复杂度O(δX × δY)，空间复杂度O(1)。
参考程序段
#include <bits/stdc++.h>
using namespace std;
class PointDistance {
public:
vector <int> findPoint( int x1, int y1, int x2, int y2 );
};
int sqr(int x){ return x * x; }
vector<int> ans;
vector <int> PointDistance::findPoint(int x1, int y1, int x2, int y2) {
ans.clear();
for(int x = -100; x <= 100; x++)
for(int y = -100; y <= 100; y++)
if(sqr(x1 - x) + sqr(y1 - y) > sqr(x2 - x) + sqr(y2 - y)){
ans.push_back(x);
ans.push_back(y);
return ans;
}
return ans;
}
TopCoder SRM 667 Div.2 T2
解题思路
可能大家的第一反应都是每一步选花费最小来贪心，但是发现，某一步有多个花费最小时，会出问题。
我们看一组数据
0000
1100
0011
0111
如果贪心， 可能顺序是1,2,3,4， 这样答案为8， 而事实上按顺序1,3,4,2,最优，为6。
于是我们考虑DP。
首先，我们令dp[i]表示在i的二进制表示下为1的那些指令处理完毕所需的最小花费。则可以从任何二进制表示下比i少一个1的j转移过来，我们取其最小。开始状态为dp[0] = 0，结束状态为dp[(1 << n) - 1]。
参考程序段
程序中，并非枚举dp[j]向dp[i]转移，而是对于dp[j]向每个可能的dp[i]转移，并使用了类似于spfa的写法。
同时，程序中q数组保存的数指向dp数组的下标，q_used对应花费。
#include <bits/stdc++.h>
using namespace std;
class OrderOfOperationsDiv2 {
public:
int minTime( vector <string> s );
};
int n;
int rec[20];
int l, r, q[1 << 23], q_used[1 << 23];
int dp[1 << 21];
int len;
int OrderOfOperationsDiv2::minTime(vector <string> s) {
n = s.size();
len = s[0].size();
for(int i = 0; i < n; i++){
rec[i] = 0;
for(int j = 0; j < len; j++)
if(s[i][j] == '1') rec[i] |= (1 << j);//转成二进制方便处理
}
memset(dp, 255, sizeof(dp));
dp[0] = 0;
l = 0; r = 1; q[1] = 0; q_used[1] = 0;
while(l < r){
l++;
int rec_vis = q[l];
int rec_used = q_used[l];
for(int i = 0; i < n; i++){//枚举可能被更新的状态
if((rec_vis >> i) & 1) continue;
int x = rec_vis | (1 << i);
int y = rec_used | rec[i];
int tt = 0;
for(int j = 0; j < len; j++)
if(((y >> j) & 1) == 1 && ((rec_used >> j) & 1) == 0) tt++;
tt = tt * tt;
if(dp[x] == -1 || dp[x] > dp[rec_vis] + tt){
dp[x] = dp[rec_vis] + tt;
r++;
q[r] = x;
q_used[r] = y;
}
}
}
return dp[(1 << n) - 1];
}
TopCoder SRM 667 Div.2 T3
解题思路
乍一看好像很麻烦呐。看起来有后效性？仿佛不能用dp？唉？然而数据只有30？所以考虑以下乱搞？
emmm
我们还是dp，大不了扩大维数，增加状态数。嗯。我们令dp[i][j][k]表示前i幢楼，在第i幢楼这里开了j家店，同时计划在下一幢楼开k家能获得的利益。那么我们可以得到状态方程
\[dp[i][j][k] = max(dp[i][j][k], dp[i - 1][l][j] + j * c[i * 3 * m + (j + k + l) - 1])(0 <= l <= m)
\]
特判一下第一栋楼和最后一栋楼。
同时注意下小细节
1、题目里并没有给出c[-1]的定义，所以任意相邻三栋楼中一定要开一家：）
2、只有一栋楼的情况：）
时间复杂度O(n^4), 空间复杂度O(n^3)【这个我觉得还是叫暴力比较好：捂脸：】
参考代码段
#include <bits/stdc++.h>
using namespace std;
class ShopPositions {
public:
int maxProfit( int n, int m, vector <int> c );
};
int f[40][40][40];
int ShopPositions::maxProfit(int n, int m, vector <int> c) {
memset(f, 0, sizeof(f));
int ans = 0;
for(int j = 0; j <= m; j++)//特判第一栋楼
for(int k = 0; k <= m; k++){
if(j + k == 0) continue;
f[0][j][k] = max(f[0][j][k], j * c[0 * 3 * m + (j + k) - 1]);
ans = max(ans, f[0][j][k]);
}
for(int i = 1; i < n - 1; i++)
for(int j = 0; j <= m; j++)
for(int k = 0; k <= m; k++)//O(n^4)DP
for(int l = 0; l <= m; l++){
if(j + k + l == 0) continue;
f[i][j][k] = max(f[i][j][k], f[i - 1][l][j] + j * c[i * 3 * m + (j + k + l) - 1]);
ans = max(ans, f[i][j][k]);
}
if(n > 1)//如果不止一栋楼，就特判最后一栋楼
for(int j = 0; j <= m; j++)
for(int l = 0; l <= m; l++){
if(j + l == 0) continue;
f[n - 1][j][0] = max(f[n - 1][j][0], f[n - 2][l][j] + j * c[(n - 1) * 3 * m + (j + l) - 1]);
ans = max(ans, f[n - 1][j][0]);
}
return ans;
}
