hdu 状态压缩dp小集(中低难度)

源自---------------------------------------------------- Ice_Crazy

hdu1074
    状态压缩dp入门题,关键是输出顺序。

#include"stdio.h"
#include"string.h"
#include"stdlib.h"

int n,base[20];
struct node
{
    char name[111];
    int dead,cost;
}MEM[20];
struct node1
{
    int days,score;
    int pre,now;
}DP[1<<15];

int main()
{
    int T;
    int i,l;
    int tot,t,reduce;
    for(i=0;i<15;i++)    base[i]=1<<i;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(i=0;i<n;i++)    scanf("%s%d%d",MEM[i].name,&MEM[i].dead,&MEM[i].cost);

        memset(DP,127,sizeof(DP));
        DP[0].days=0;
        DP[0].score=0;
        DP[0].pre=-1;
        tot=1<<n;
        for(i=0;i<tot;i++)
        {
            for(l=n-1;l>=0;l--)        //从大到小,并且在后面的第七行用了“<”:保证顺序
            {
                if(i & base[l])
                {
                    t=i-base[l];
                    reduce=DP[t].days+MEM[l].cost-MEM[l].dead;
                    if(reduce<0)    reduce=0;
                    if(DP[t].score+reduce<DP[i].score)
                    {
                        DP[i].pre=t;
                        DP[i].now=l;
                        DP[i].days=DP[t].days+MEM[l].cost;
                        DP[i].score=DP[t].score+reduce;
                    }
                }
            }
        }

        int K=0,temp=(1<<n)-1;
        char ans[20][111];
        while(DP[temp].pre!=-1)
        {
            strcpy(ans[K++],MEM[DP[temp].now].name);
            temp=DP[temp].pre;
        }
        printf("%d\n",DP[(1<<n)-1].score);
        for(i=K-1;i>=0;i--)    printf("%s\n",ans[i]);
    }
    return 0;
}




hdu3182
    状态压缩dp水题。




hdu4568
    最短路+状态压缩dp、不难,两部分都较果。




hdu4336
    概率dp+状态压缩 求期望。
    没啥难的、主要是关于期望的运算。。期望。。高中的数学知识
啊、我推了半张演草纸才想起来啊。。。

#include"iostream"  
#include"cstdio"  
#include"cstring"  
using namespace std;  
const int N=1<<20;  
  
int n,base[20];  
double p[20],dp[N];  
int main()  
{  
    int i,l;  
    int limit;  
    double t,sum;  
    for(i=0;i<20;i++)    base[i]=1<<i;  
    while(scanf("%d",&n)!=-1)  
    {  
        for(i=0;i<n;i++) scanf("%lf",&p[i]);  
        limit=1<<n;  
        dp[0]=0;  
        for(i=1;i<limit;i++)  
        {  
            t=sum=0;  
            for(l=0;l<n;l++)  
            {  
                if(!(i & base[l]))  continue;  
                t+=p[l]*dp[i-base[l]];  
                sum+=p[l];  
            }  
            dp[i]=(t+1)/sum;  
        }  
        printf("%.5lf\n",dp[(1<<n)-1]);  
    }  
    return 0;  
} 




hdu1565
    很果的二进制dp。
    不过这题的数据小,所以可以这样暴力,好的方法是网络流,
可以做下hdu1569。




hdu2167
    果果的状态压缩dp。
    可以先做一下hdu1565,那个比这个更果。。
    这个题的话可以先算一下,15位的格子只有1597个可用来
递推的状态,所以,尽情的暴力吧~!= =II。




hdu4539
    很水的状态压缩dp,数据也不强,4层for循环写的pl的可以迅速过之,当然不pl也能过= =。




hdu3001
    三进制dp,思路同二进制。

#include"iostream"
#include"cstdio"
#include"cstring"
using namespace std;
const int N=60000;
const int M=1000;
const int inf=123456789;

int n,m,map[20][20],dp[N][10],sp[N][10],base[10];
struct Edge{
    int v,len,next;
}edge[M];
int tot,head[10];
void add(int a,int b,int c){
    edge[tot].v=b;edge[tot].len=c;edge[tot].next=head[a];head[a]=tot++;
}

void init()
{
    int i,k,temp,limit;
    base[0]=1;
    for(i=1;i<10;i++)    base[i]=base[i-1]*3;
    limit=base[9]*3;
    memset(sp,0,sizeof(sp));
    for(i=0;i<limit;i++)
    {
        k=0;temp=i;
        while(temp)    {sp[i][k++]=temp%3;temp/=3;}
    }
}
void build_map()
{
    int i,l;
    int a,b,c;
    for(i=1;i<=n;i++)
    {
        for(l=1;l<=n;l++)    map[i][l]=inf;
        map[i][i]=0;
    }
    while(m--)
    {
        scanf("%d%d%d",&a,&b,&c);
        map[a][b]=map[b][a]=map[a][b]>c?c:map[a][b];
    }
    tot=0;
    memset(head,-1,sizeof(head));
    for(i=1;i<=n;i++)
    for(l=1;l<=n;l++)
    {
        if(i==l || map[i][l]==inf)    continue;
        add(i-1,l-1,map[i][l]);
    }
}
int DP()
{
    int i,l,j,v;
    int t,next,flag,limit,ans;
    limit=1;
    for(i=0;i<n;i++)    limit*=3;
    for(i=0;i<limit;i++)
    for(l=0;l<n;l++)
        dp[i][l]=inf;
    t=1;
    for(i=0;i<n;i++)    {dp[t][i]=0;t*=3;}

    ans=inf;
    for(i=0;i<limit;i++)
    {
        flag=1;
        for(l=0;l<n;l++)
        {
            if(!sp[i][l])    {flag=0;continue;}
            if(dp[i][l]==inf)    continue;
            for(j=head[l];j!=-1;j=edge[j].next)
            {
                v=edge[j].v;
                if(sp[i][v]>=2)    continue;
                next=i+base[v];
                if(dp[i][l]+edge[j].len<dp[next][v])
                    dp[next][v]=dp[i][l]+edge[j].len;
            }
        }
        if(!flag)    continue;
        for(l=0;l<n;l++)if(dp[i][l]<ans)ans=dp[i][l];
    }
    if(ans==inf)    ans=-1;
    return ans;
}
int main()
{
    init();
    while(scanf("%d%d",&n,&m)!=-1)
    {
        build_map();
        printf("%d\n",DP());
    }
    return 0;
}


  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
### 回答1: hdu 2829 Lawrence 斜率优化dp 这道题是一道经典的斜率优化dp题目,需要用到单调队列的思想。 题目大意是给定一个序列a,求出一个序列b,使得b[i]表示a[1]~a[i]中的最小值,且满足b[i] = min{b[j] + (i-j)*k},其中k为给定的常数。 我们可以将上式拆开,得到b[i] = min{b[j] - j*k} + i*k,即b[i] = i*k + min{b[j] - j*k},这个式子就是斜率优化dp的形式。 我们可以用单调队列来维护min{b[j] - j*k},具体思路如下: 1. 首先将第一个元素加入队列中。 2. 从第二个元素开始,我们需要将当前元素加入队列中,并且需要维护队列的单调性。 3. 维护单调性的方法是,我们从队列的末尾开始,将队列中所有大于当前元素的元素弹出,直到队列为空或者队列中最后一个元素小于当前元素为止。 4. 弹出元素的同时,我们需要计算它们对应的斜率,即(b[j]-j*k)/(j-i),并将这些斜率与当前元素的斜率比较,如果当前元素的斜率更小,则将当前元素加入队列中。 5. 最后队列中的第一个元素就是min{b[j] - j*k},我们将它加上i*k就得到了b[i]的值。 6. 重复以上步骤直到处理完所有元素。 具体实现可以参考下面的代码: ### 回答2: HDU 2829 Lawrence 斜率优化 DP 是一道经典的斜率优化 DP 题目,其思想是通过维护一个下凸包来优化 DP 算法。下面我们来具体分析一下这道题目。 首先,让我们看一下该题目的描述。题目给定一些木棒,要求我们将这些木棒割成一些给定长度,且要求每种长度的木棒的数量都是一样的,求最小的割枝次数。这是一个典型的背包问题,而且在此基础上还要求每种长度的木棒的数量相同,这就需要我们在状态设计上走一些弯路。 我们来看一下状态的定义。定义 $dp[i][j]$ 表示前 $i$ 个木棒中正好能割出 $j$ 根长度为 $c_i$ 的木棒的最小割枝次数。对于每个 $dp[i][j]$,我们可以分类讨论: 1. 不选当前的木棒,即 $dp[i][j]=dp[i-1][j]$; 2. 选当前的木棒,即 $dp[i][j-k]=dp[i-1][j-k]+k$,其中 $k$ 是 $j/c_i$ 的整数部分。 现在问题再次转化为我们需要在满足等量限制的情况下,求最小的割枝次数。可以看出,这是一个依赖于 $c_i$ 的限制。于是,我们可以通过斜率优化 DP 来解决这个问题。 我们来具体分析一下斜率优化 DP 算法的思路。我们首先来看一下动态规划的状态转移方程 $dp[i][j]=\min\{dp[i-1][k]+x_k(i,j)\}$。可以发现,$dp[i][j]$ 的最小值只与 $dp[i-1][k]$ 和 $x_k(i,j)$ 有关。其中,$x_k(i,j)$ 表示斜率,其值为 $dp[i-1][k]-k\times c_i+j\times c_i$。 接下来,我们需要维护一个下凸包,并通过斜率进行优化。我们具体分析一下该过程。假设我们当前要计算 $dp[i][j]$。首先,我们需要找到当前点 $(i,j)$ 在凸包上的位置,即斜率最小值的位置。然后,我们根据该位置的斜率计算 $dp[i][j]$ 的值。接下来,我们需要将当前点 $(i,j)$ 加入到下凸包上。 我们在加入点的时候需要注意几点。首先,我们需要将凸包中所有斜率比当前点小的点移除,直到该点能够加入到凸包中为止。其次,我们需要判断该点是否能够加入到凸包中。如果不能加入到凸包中,则直接舍弃。最后,我们需要保证凸包中斜率是单调递增的,这就需要在加入新的点之后进行上一步操作。 以上就是该题目的解题思路。需要注意的是,斜率优化 DP 算法并不是万能的,其使用情况需要根据具体的问题情况来确定。同时,该算法中需要维护一个下凸包,可能会增加一些算法的复杂度,建议和常规 DP 算法进行对比,选择最优的算法进行解题。 ### 回答3: 斜率优化DP是一种动态规划优化算法,其主要思路是通过对状态转移方程进行变形,提高算法的时间复杂度。HDU2829 Lawrence问题可以用斜率优化DP解决。 首先,我们需要了解原问题的含义。问题描述如下:有$n$个人在数轴上,第$i$个人的位置为$A_i$,每个人可以携带一定大小的行李,第$i$个人的行李重量为$B_i$,但是每个人只能帮助没有他们重量大的人搬行李。若第$i$个人搬运了第$j$个人的行李,那么第$i$个人会累加$C_{i,j}=\left|A_i-A_j\right|\cdot B_j$的体力消耗。求$m$个人帮助每个人搬运行李的最小体力消耗。 我们可以通过斜率优化DP解决这个问题。记$f_i$为到前$i$个人的最小体力消耗,那么状态转移方程为: $$f_i=\min_{j<i}\{f_j+abs(A_i-A_j)\cdot B_i\}$$ 如果直接使用该方程,时间复杂度为$O(n^2)$,如果$n=10^4$,则需要计算$10^8$次,运算时间极长。斜率优化DP通过一些数学推导将方程变形,将时间复杂度降低到$O(n)$,大大缩短了计算时间。 通过斜率优化DP的推导式子,我们可以得到转移方程为: $$f_i=\min_{j<i}\{f_j+slope(j,i)\}$$ 其中,$slope(j,i)$表示直线$j-i$的斜率。我们可以通过如下方式来求解$slope(j,i)$: $$slope(j,i)=\frac{f_i-f_j}{A_i-A_j}-B_i-B_j$$ 如果$slope(j,i)\leq slope(j,k)$,那么$j$一定不是最优,可以直接舍去,降低计算时间。该算法的时间复杂度为$O(n)$。 综上所述,斜率优化DP是一种动态规划优化算法,可以大大缩短计算时间。在处理类似HDU2829 Lawrence问题的时候,斜率优化DP可以很好地解决问题。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值