SGU132 状态压缩~最优值
2013-06-06 18:38
295 查看
看了一天,然后由于各种事情一直没写,今天搞了一下午才搞出来,每次存两行的状态,保证now状态之前的行满足条件。
ACcode:
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int NS=80;
const int MS=10;
const int MAX=7000;
int n,m,c,sj,sx;
char str[NS];
int dp[2][MAX];
int p3[MS],g[NS];
int bit[MAX][MS];
void prepare()
{
p3[0]=1;
for (int i=1;i<MS;i++)
p3[i]=p3[i-1]*3;
for (int i=0;i<MAX;i++)
{
for (int t=i,j=0;t>0&&j<MS;j++)
bit[i][j]=t%3,t/=3;
}
}
void init(int x)
{
for (int i=0;i<p3[m+1];i++)
dp[x][i]=MAX;
}
int Min(int a1,int b1)
{
return a1<b1?a1:b1;
}
void dfs(int pos,int old,int now,int sum)
{
if (pos>0&&bit[old][pos-1]==2) return ;
if (pos>1&&bit[old][pos-1]==0&&bit[old][pos-2]==0) return ;
if (pos==m)
{
dp[sx][now]=Min(dp[sx^1][sj]+sum,dp[sx][now]);
return ;
}
if ( bit[old][pos]!=1&&bit[now][pos]!=1 )
dfs(pos+1,old-bit[old][pos]*p3[pos]+p3[pos],now+p3[pos],sum+1);
if ( pos+1<m&&bit[now][pos]!=1&&bit[now][pos+1]!=1 )
dfs(pos+1,old,now+p3[pos]+p3[pos+1],sum+1);
if (bit[old][pos]==0&&bit[now][pos]==0)
now=now+p3[pos]+p3[pos];
dfs(pos+1,old,now,sum);
}
int main()
{
prepare();
scanf("%d%d",&n,&m);
{
for (int i=0;i<n;i++)
{
g[i]=0;
scanf("%s",str);
for (int j=0;j<m;j++)
g[i]+=str[j]=='*'?p3[j]:0;
}
int pre=1,cur=0;
init(cur),c=0;
for (int i=0;i<m;i++)
c+=p3[i];
dp[cur][c]=0,g
=c;
for (int i=0;i<=n;i++)
{
pre^=1,cur^=1;
init(cur);
for (int j=0;j<p3[m+1];j++)
if (dp[pre][j]<MAX)
{
sj=j,sx=cur;
dfs(0,j,g[i],0);
}
}
printf("%d\n",dp[cur][c]);
}
return 0;
}
ACcode:
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int NS=80;
const int MS=10;
const int MAX=7000;
int n,m,c,sj,sx;
char str[NS];
int dp[2][MAX];
int p3[MS],g[NS];
int bit[MAX][MS];
void prepare()
{
p3[0]=1;
for (int i=1;i<MS;i++)
p3[i]=p3[i-1]*3;
for (int i=0;i<MAX;i++)
{
for (int t=i,j=0;t>0&&j<MS;j++)
bit[i][j]=t%3,t/=3;
}
}
void init(int x)
{
for (int i=0;i<p3[m+1];i++)
dp[x][i]=MAX;
}
int Min(int a1,int b1)
{
return a1<b1?a1:b1;
}
void dfs(int pos,int old,int now,int sum)
{
if (pos>0&&bit[old][pos-1]==2) return ;
if (pos>1&&bit[old][pos-1]==0&&bit[old][pos-2]==0) return ;
if (pos==m)
{
dp[sx][now]=Min(dp[sx^1][sj]+sum,dp[sx][now]);
return ;
}
if ( bit[old][pos]!=1&&bit[now][pos]!=1 )
dfs(pos+1,old-bit[old][pos]*p3[pos]+p3[pos],now+p3[pos],sum+1);
if ( pos+1<m&&bit[now][pos]!=1&&bit[now][pos+1]!=1 )
dfs(pos+1,old,now+p3[pos]+p3[pos+1],sum+1);
if (bit[old][pos]==0&&bit[now][pos]==0)
now=now+p3[pos]+p3[pos];
dfs(pos+1,old,now,sum);
}
int main()
{
prepare();
scanf("%d%d",&n,&m);
{
for (int i=0;i<n;i++)
{
g[i]=0;
scanf("%s",str);
for (int j=0;j<m;j++)
g[i]+=str[j]=='*'?p3[j]:0;
}
int pre=1,cur=0;
init(cur),c=0;
for (int i=0;i<m;i++)
c+=p3[i];
dp[cur][c]=0,g
=c;
for (int i=0;i<=n;i++)
{
pre^=1,cur^=1;
init(cur);
for (int j=0;j<p3[m+1];j++)
if (dp[pre][j]<MAX)
{
sj=j,sx=cur;
dfs(0,j,g[i],0);
}
}
printf("%d\n",dp[cur][c]);
}
return 0;
}
相关文章推荐
- SGU 132 一道很恶心的状态压缩dp
- SGU132 - Another Chocolate Maniac(状态压缩DP)
- SGU - 132 Another Chocolate Maniac (状态压缩)
- SGU 132 状态压缩
- SGU 132 Another Chocolate Maniac 状态压缩DP
- SGU 131 贴地砖类型 状态压缩DP
- hdu4284 状态压缩~最优值
- SGU 536 Berland Chess(状态压缩 + bfs)
- hdu 1074 状态压缩~最优值
- hdu2167 状态压缩~最优值
- sgu 223 Little Kings 状态压缩dp
- sgu 222 - Little Rooks 简单的状态压缩DP
- SGU 131.Hardwood floor 状态压缩DP
- SGU223 -状态压缩DP
- SGU 197 状态压缩DP+矩阵乘法
- hdu1565 状态压缩~最优值
- hdu4281 状态压缩~MTSP最优值
- SGU - 131 Hardwood floor (状态压缩)
- SGU 327 Yet Another Palindrome(状态压缩DP)
- 状态压缩dp/sgu 131 Hardwood floor