这道题以为是简单的模板题(确实简单啊。。。),卡了好久
1.内存限制,之前题放的很开,这里要卡一下最多的状态数,但没卡滚动数组(也卡不了吧)
2.注意要统计的是最值还是数量。
3.二进制状态
2n和2n+1
2
n
和
2
n
+
1
是分开的 可以预处理
220
2
20
再对不同n取前缀,但效果不明显。
#include<iostream>
#include<stdio.h>
#include<string>
#include<cstring>
#include<vector>
#include<stack>
#include<algorithm>
using namespace std;
#define clr(a, x) memset(a, x, sizeof(a))
#define mp(x, y) make_pair(x, y)
#define pb(x) push_back(x)
#define X first
#define Y second
#define fastin \
ios_base::sync_with_stdio(0); \
cin.tie(0);
typedef long long LL;
typedef pair<int, int> PII;
typedef vector<int> VI;
const int INF=0x3f3f3f3f;
LL G[30][30];
LL state[20000];
LL sum[20][20000];
LL dp[20][20000];
int n;
LL dfs(int r,int id){
if(sum[r][id]) return sum[r][id];
int ans=0;
for(int i=0;i<n;i++)
{
if(state[id]&(1<<i))
ans+=G[r][i];
}
return ans;
}
int main()
{
int num=0;
for(int i=0;i<(1<<20);i++)
if(!(i&(i<<1))) state[num++] = i;
while(cin>>n){
clr(sum,0);
clr(dp,0);
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
scanf("%d",&G[i][j]);
int total=(1<<n);
int id=0;
for(;id<num;id++) if(state[id]>total) break;
LL ans=0;
for(int i=0;i<id;i++)
{
dp[0][i] = dfs(0,i);
ans=max(ans,dp[0][i]);
}
for(int r=0;r<n-1;r++)
for(int i=0;i<id;i++)
if(dp[r][i])
for(int j=0;j<id;j++)
{
if(state[i]&state[j]) continue;
dp[r+1][j] = max(dp[r][i] + dfs(r+1,j),dp[r+1][j]);
ans=max(ans,dp[r+1][j]);
}
cout<<ans<<endl;
}
return 0;
}
发现一个轮廓线dp做法 有意思
初始值(初始值的选取也有讲究)
0 | 0 | 0 | 0 | 0 |
---|
需要更新的值
x | x | 0 | 1 | 0 |
---|---|---|---|---|
1 | 0 | ? | x | x |
状态
n2−1
n
2
−
1
就是形似上图所示的意义。
对于一个
n×n
n
×
n
的矩阵,从行,列一一遍历,同时更新每个可到达的状态的最大值。其他和入门题一致,比较简单。
#include<iostream>
#include<string>
#include<cstring>
#include<stdio.h>
#include<algorithm>
using namespace std;
#define mp(x, y) make_pair(x, y)
#define pb(x) push_back(x)
#define X first
#define Y second
#define fastin \
ios_base::sync_with_stdio(0); \
cin.tie(0);
#define clr(a, x) memset(a, x, sizeof(a))
typedef long long LL;
typedef pair<int, int> PII;
const int INF=0x3f3f3f3f;
int n,cur;
LL G[25][25];
LL dp[2][1<<20];
int main()
{
while(~scanf("%d",&n)){
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
scanf("%lld",&G[i][j]);
clr(dp,0);
cur=0;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
{
cur=1^cur;
clr(dp[cur],0);
for(int st=0;st<(1<<n);st++){
//取
if( (j==0||(st&1)==0) && (st&(1<<(n-1)))==0) //已经保证首位不是1
dp[cur][(st<<1)|1] = max(dp[cur][(st<<1)|1],dp[1^cur][st]+G[i][j]);
//不取
dp[cur][(st<<1)&(1<<n)-1] = max(dp[cur][(st<<1)&(1<<n)-1],dp[1^cur][st]);
}
}
LL ans=0;
for(int i=0;i<(1<<n);i++)
ans=max(ans,dp[cur][i]);
printf("%lld\n",ans);
}
return 0;
}