BZOJ 4152 状压DP 解题报告

来源:互联网 发布:linux 设置ip 编辑:程序博客网 时间:2024/05/21 12:40

4152: [AMPPZ2014]The Captain

Time Limit:

给定平面上的n个点,定义(x1,y1)到(x2,y2)的费用为min(|x1-x2|,|y1-y2|),求从1号点走到n号点的最小费用。

Input

第一行包含一个正整数n(2<=n<=200000),表示点数。
接下来n行,每行包含两个整数x[i],yi,依次表示每个点的坐标。

Output

一个整数,即最小费用。

Sample Input

5
2 2
1 1
4 5
7 1
6 7

Sample Output

2

【解题报告】

一看题感觉是状压,一看数据范围更加确定是状压,当然打表也可以算作下下策。

代码如下:

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include<cmath>#include<vector> using namespace std;int n,m,k,stay[101],cnt[101];long long ans,dp[10][101][101];//行,王,状态 bool map[101][101];void dfs(int p,int put,int num){    stay[++m]=num;    cnt[m]=p;    if(p>=k||p>=(n+1)/2)return;    for(int i=put+2;i<=n;i++)    dfs(p+1,i,num+(1<<(i-1)));}//枚举出在一行上放马的所有合法状态int main(){    memset(dp,0,sizeof(dp));    memset(stay,0,sizeof(stay));    memset(cnt,0,sizeof(cnt));    memset(map,0,sizeof(map));    scanf("%d%d",&n,&k);    m=0;    dfs(0,-1,0);    for(int i=1;i<=m;i++)    for(int j=1;j<=m;j++)    map[i][j]=map[j][i]=((stay[i]&stay[j])||((stay[i]<<1)&stay[j])||((stay[i]>>1)&stay[j]))?0:1;    for(int i=1;i<=m;i++)    dp[1][cnt[i]][i]=1ll;    for(int i=2;i<=n;i++)    for(int j=0;j<=k;j++)    for(int now=1;now<=m;now++)    {        if(cnt[now]>j)continue;        for(int l=1;l<=m;l++)        if(map[now][l]&&cnt[l]+cnt[now]<=j) dp[i][j][now]+=dp[i-1][j-cnt[now]][l];    }    ans=0;    for(int i=1;i<=m;i++)ans+=dp[n][k][i];    printf("%lld\n",ans);       return 0;}