状态DP

于 2014-07-15 10:35:01 发布
阅读量699 收藏
点赞数
文章标签： 状态DP TSP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011495472/article/details/37809813
版权

由于队友XJP在做状态DP，我也跟着做了一发，作为数学系的学生，记得两年前在运筹学中学到了一点点动态规划的皮毛，当初对动态规划不屑一顾，后来接触了背包算法，以及各种编程中用到了不少动态规划的知识才发现动态规划之博大精深，不是动态规划不强大，而是自己太弱了，学习动态规划是学习了这里的文章。



做的第一道状态DP是  zoj 3777Problem Arrangement   题意是随机安排题目使得兴趣值大于给定M<=300的期望，题目数量1<=N<=12

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=13;
const int M=502;

int dp[1<<N][M];
int a[N][N],f[N];

void init(){
    f[0]=1;
    for(int i=1;i!=N;++i){
        f[i]=f[i-1]*i;
    }
}

int gcd(int a,int b){
    return b?gcd(b,a%b):a;
}

int main(){
    int T,n,m;
    init();
//    freopen("in.txt","r",stdin);
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=0;i!=n;++i){
            for(int j=0;j!=n;++j){
                scanf("%d",&a[i][j]);
            }
        }
        memset(dp,0,sizeof(dp));
        dp[0][0]=1;
        for(int i=0;i!=(1<<n);++i){
            int tmp=0;
            for(int j=0;j!=n;++j){
                if(i&(1<<j))    ++tmp;
            }
            for(int j=0;j!=n;++j){
                if(i&(1<<j))    continue;
                for(int k=0;k<=m;++k){
                    if(k+a[tmp][j]>=m)      dp[i|(1<<j)][m]+=dp[i][k];
                    else                    dp[i|(1<<j)][k+a[tmp][j]]+=dp[i][k];
                }
            }
        }
        if(dp[(1<<n)-1][m]){
            int tmp=gcd(dp[(1<<n)-1][m],f[n]);
            printf("%d/%d\n",f[n]/tmp,dp[(1<<n)-1][m]/tmp);
        }else printf("No solution\n");
    }
    return 0;
}
第二题  poj3254Corn Fields 题意在可以种草的地方 种草且草不相邻的种数mod1e9 

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=13;
const int MOD=100000000;

int a[N],dp[1<<N][N],next[1<<N];

int main(){
    int n,m,tmp;
//    freopen("in.txt","r",stdin);
    while(scanf("%d%d",&n,&m)!=EOF){
        int mm=1<<m;
        memset(dp,0,sizeof(dp));
        memset(a,0,sizeof(a));
        tmp=0;
        for(int i=0;i!=mm;++i){
            if((i&(i>>1))==0){
                next[tmp]=i;tmp=i;
            }
        }
        next[tmp]=mm;

        for(int i=0;i!=n;++i){
            for(int j=0;j!=m;++j){
                scanf("%d",&tmp);
                if(tmp==0)  a[i]|=(1<<j);
            }
        }
        for(int i=0;i!=mm;i=next[i]){
            if((i&a[0])==0)   dp[i][0]=1;
        }
        for(int j=1;j!=n;++j){
            for(int i=0;i!=mm;i=next[i]){
                for(int k=0;k!=mm;k=next[k]){
                    if((k&a[j])==0 && (i&k)==0)   dp[k][j]=(dp[k][j]+dp[i][j-1])%MOD;
                }
            }
        }
        int ans=0;
        for(int i=0;i!=mm;i=next[i]){
            ans=(ans+dp[i][n-1])%MOD;
        }
        printf("%d\n",ans);
    }
    return 0;
}

第三题   poj1185 炮兵阵地  题意就是在平原放能放置的最大大炮数且满足同行同列的大炮之间距离大于2 

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
typedef __int64 LL;
using namespace std;
const int N=11;

int a[102],next[1<<N],p[70];
int dp[102][70][70],num[1<<N];
char s[102];

int bcount(int x){
   int cnt=0;
   while(x){
       ++cnt;
       x&=(x-1);
   }
   return cnt;
}

inline int max(int a,int b){
    return a>b?a:b;
}

int main(){
    int n,m,tmp;
    //freopen("in.txt","r",stdin);
    while(scanf("%d%d",&n,&m)!=EOF){
        int mm=1<<m;
        memset(dp,0,sizeof(dp));
        memset(a,0,sizeof(a));
        tmp=0;
        for(int i=0;i!=mm;++i){
            if((i&(i>>1))==0 && (i&(i>>2))==0){
                next[tmp]=i;tmp=i;
            }
        }
        next[tmp]=mm;
        int len=0;
        for(int i=0;i!=mm;i=next[i])    p[len++]=i;
        for(int i=0;i!=len;++i) num[p[i]]=bcount(p[i]);
        for(int i=0;i!=n;++i){
            scanf("%s",s);
            for(int j=0;j!=m;++j){
                if(s[j]=='H')  a[i]|=(1<<j);
            }
        }

        for(int i=0;i!=len;++i){
            if((p[i]&a[0])==0)   dp[0][0][i]=num[p[i]];
        }
        if(n>1) for(int i=0;i!=len;++i){
            for(int j=0;j!=len;++j){
                if((p[j]&a[1])==0 && (p[i]&p[j])==0)   dp[1][i][j]=dp[0][0][i]+num[p[j]];
            }
        }
        for(int i=2;i<n;++i){
            for(int j=0;j!=len;++j){
                for(int k=0;k!=len;++k){
                    if(p[j]&p[k])   continue;
                    for(int l=0;l!=len;++l){
                        if(p[j]&p[l])   continue;
                        if(p[k]&p[l])   continue;
                        if(a[i]&p[l])   continue;
                        dp[i][k][l]=max(dp[i][k][l],dp[i-1][j][k]+num[p[l]]);
                    }
                }
            }
        }
        int ans=0;
        for(int i=0;i!=len;++i){
            for(int j=0;j!=len;++j){
                ans=max(ans,dp[n-1][i][j]);
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}
第四题  TSP旅行商问题，不过有点小变化   poj3311Hide with the pie 
#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
typedef __int64 LL;
using namespace std;

const int N=12;
const int INF=1e8;

int a[N][N],dp[1<<N][N];

inline int min(int x,int y){
    return x<y?x:y;
}

int main(){
    int n;
    while(scanf("%d",&n)!=EOF&&n){
        ++n;
        int nn=1<<n;
        for(int i=0;i!=n;++i){
            for(int j=0;j!=n;++j){
                scanf("%d",&a[i][j]);
            }
        }
        for(int k=0;k!=n;++k){  //Floyd
            for(int i=0;i!=n;++i){
                for(int j=0;j!=n;++j){
                    a[i][j]=min(a[i][j],a[i][k]+a[k][j]);
                }
            }
        }
        for(int i=0;i!=nn;++i){
            for(int j=0;j!=n;++j){
                dp[i][j]=INF;
            }
        }
        dp[1][0]=0;
        for(int i=0;i!=nn;++i){
            for(int j=0;j!=n;++j){
                if(i&(1<<j))    continue;
                for(int k=0;k!=n;++k){
                    if(i&(1<<k)) dp[i|(1<<j)][j]=min(dp[i|(1<<j)][j],dp[i][k]+a[k][j]);
                }
            }
        }
        int ans=INF;
        for(int i=0;i!=n;++i){
            ans=min(ans,dp[nn-1][i]+a[i][0]);
        }
        printf("%d\n",ans);
    }
    return 0;
}
第四题 hdu3681Prison Break本质上与上题没啥区别，只是需要转化一下把距离求出来，并且这里需要二分查找答案 

#include <iostream>
#include <cstdio>
#include <vector>
#include <queue>
#include <cstring>
#include <algorithm>
typedef __int64 LL;
using namespace std;

const int N=16;
struct Node{
    int x,y;
    Node(int _x=0,int _y=0):x(_x),y(_y){}
}a[N];

char str[N];
int map[N][N];
int dt[N][N][N][N],n,m,len,ans;
int dp[1<<N][N],dist[N][N];
int dir[4][2]={{-1,0},{1,0},{0,1},{0,-1}};

void bfs(int s){
    queue<Node> q;
    while(!q.empty())   q.pop();
    int sx=a[s].x,sy=a[s].y;
    dt[sx][sy][sx][sy]=0;
    q.push(a[s]);
    while(!q.empty()){
        Node cur=q.front();
        q.pop();
        int x=cur.x,y=cur.y;
        for(int i=0;i!=4;++i){
            int tx=x+dir[i][0];
            int ty=y+dir[i][1];
            if(tx<0||tx>=m||ty<0||ty>=m||map[tx][ty]<0)  continue;
            if(dt[sx][sy][tx][ty]==-1){
                dt[sx][sy][tx][ty]=dt[sx][sy][x][y]+1;
                q.push(Node(tx,ty));
            }
        }
    }
}
bool check(int num){
    int top=1<<len;
    memset(dp,-1,sizeof(dp));
    dp[1][0]=num;
    for(int i=0;i!=top;++i){
        for(int j=0;j!=len;++j){
            if((i&(1<<j))==0||dp[i][j]==-1)    continue;
            if((i&ans)==ans)    return true;
            for(int k=0;k!=len;++k){
                if((i&(1<<k))||dist[j][k]==-1||i==j) continue;
                dp[i|(1<<k)][k]=max(dp[i|(1<<k)][k],dp[i][j]-dist[j][k]);
                if(dp[i|(1<<k)][k]>=0&&map[a[k].x][a[k].y]==2)  dp[i|(1<<k)][k]=num;
            }
        }
    }
    return false;
}
int main(){
//    freopen("in.txt","r",stdin);
    while(scanf("%d%d",&n,&m)!=EOF&&(n+m)){
        len=1;ans=0;
        for(int i=0;i!=n;++i){
            scanf("%s",str);
            for(int j=0;j!=m;++j){
                if(str[j]=='S')         map[i][j]=0;
                else if(str[j]=='D')    map[i][j]=-1;
                else if(str[j]=='F'){
                    map[i][j]=0;
                    a[0]=Node(i,j);
                }else{
                    a[len++]=Node(i,j);
                    if(str[j]=='Y'){
                        map[i][j]=1;
                        ans|=1<<(len-1);
                    }
                    else                map[i][j]=2;
                }
            }
        }
        memset(dt,-1,sizeof(dt));
        for(int i=0;i!=len;++i) bfs(i);
        for(int i=0;i!=len;++i){
            for(int j=0;j!=i;++j){
                dist[j][i]=dist[i][j]=dt[a[i].x][a[i].y][a[j].x][a[j].y];
            }
            dist[i][i]=0;
        }
        int low=0,high=303;
        while(low<=high){
            int mid=(low+high)>>1;
            if(check(mid))  high=mid-1;
            else            low=mid+1;
        }
        if(low==304)    printf("-1\n");
        else            printf("%d\n",low);
    }
    return 0;
}

总结可以认为第一个是无限状态DP，第二个是1状态DP，第三个是2状态，类似于背包前两个比较好处理，而最后一个处理起来就相对较麻烦，另外状态dp适合于求解n!以及2^n的复杂度问题，它并没有降低复杂度，只是减少了枚举，并能很方便的解决小规模问题



dna049
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
状态压缩DP
04-01
例如，在0/1背包问题中，状态转移方程为`dp[i][j] = max(dp[i-1][j], dp[i-1][j-w[i]] + v[i])`，其中`w[i]`是第i个物品的重量，`v[i]`是第i个物品的价值。 #### 四、状态压缩动态规划的应用场景 **1. 棋盘问题** ...
状态压缩dp
03-15
状态图压缩入门的资料，如果谁需要，谁就下吧！
dp状态压缩资料
08-08
状态压缩是一种在动态规划（Dynamic Programming, 简称DP）中优化空间复杂度的技巧。在处理一些具有大量状态的DP问题时，如果每个状态之间有明显的位运算关系，或者状态之间的信息重叠，我们可以将多个状态压缩到一...
动态规划：状态压缩DP - 愤怒的小鸟
07-09
### 动态规划：状态压缩DP - 愤怒的小鸟 #### 题目介绍 本题源自AcWing平台编号为524的题目《愤怒的小鸟》。题目要求寻找最少数量的抛物线，使得每条抛物线都能覆盖给定的一组点中的若干个点。具体而言，给定一组...
DP 背包问题详解pdf
05-19
- **滚动数组**：由于dp[i][j]只依赖于上一行的状态dp[i-1][j]和dp[i-1][j-v[i]]，因此可以使用滚动数组来减少空间复杂度，只需维护两行状态数组即可。 - **内层循环方向**：为了防止重复计算或数据污染，内层循环...
vasp6.3.2安装包
08-19
VASP（Vienna Ab-initio Simulation Package）是一个用于电子结构计算和量子力学-分子动力学模拟的软件包，由维也纳大学的Hafner小组开发。 它基于密度泛函理论（DFT）和赝势平面波基组，广泛应用于材料模拟和计算物质科学研究领域。VASP能够处理周期性和非周期性体系，包括金属、半导体、分子、晶体、掺杂体系、复合材料、纳米材料等，提供精确的电子性质预测和理化性质分析。由于其高效的大规模并行计算能力，VASP在国际上享有很高的声誉，是第一性原理计算领域最权威的软件之一。每年，以VASP作为科研工具发表的论文数量众多，包括在顶级科学期刊上发表的研究成果。物理与电子科学学院等多名教师使用VASP进行科研工作，并取得了多项研究成果。‌
计算机毕业设计SSM员工人事工资管理系统论文
08-19
计算机毕业设计SSM员工人事工资管理系统论文
根据文档框架生成内容草稿版
08-19
根据文档框架生成内容
毕业设计，基于SpringBoot+Vue+MySQL开发的企业客户管理系统，源码+论文答辩+毕业论文+视频演示
08-19
毕业设计，基于SpringBoot+Vue+MySQL开发的企业客户管理系统，源码+论文答辩+毕业论文+视频演示 本论文主要论述了如何使用JAVA语言开发一个企业客户管理系统，本系统将严格按照软件开发流程进行各个阶段的工作，采用B/S架构，面向对象编程思想进行项目开发。在引言中，作者将论述企业客户管理系统的当前背景以及系统开发的目的，后续章节将严格按照软件开发流程，对系统进行各个阶段分析设计。 企业客户管理系统的主要使用者分为管理员和员工，实现功能包括管理员：首页、个人中心、员工管理、客户信息管理、行业类型管理、项目信息管理、项目类型管理、项目收益管理，员工：首页、个人中心、客户信息管理、项目信息管理、项目收益管理等功能。由于本网站的功能模块设计比较全面，所以使得整个企业客户管理系统信息管理的过程得以实现。 本系统的使用可以实现本企业客户管理的信息化，可以方便管理员进行更加方便快捷的管理，可以提高管理人员的工作效率。 关键词：企业客户管理系统 JAVA语言；MySQL数据库；Spring Boot框架
探索图的世界：一键检测连通性的奥秘
最新发布
08-19
数据结构是计算机科学中的一个基本概念，它指的是数据的组织、管理和存储方式，以及对数据的操作。数据结构使得数据的访问和修改更加高效和有序。常见的数据结构包括： 1. **数组**（Array）：一种线性数据结构，可以存储相同类型的元素，并通过索引访问。 2. **链表**（Linked List）：一种线性数据结构，由一系列节点组成，每个节点包含数据部分和指向下一个节点的指针。 3. **栈**（Stack）：一种后进先出（LIFO, Last In First Out）的数据结构，只能在一端进行添加或删除操作。 4. **队列**（Queue）：一种先进先出（FIFO, First In First Out）的数据结构，允许在一端添加元素，在另一端删除元素。 5. **哈希表**（Hash Table）：通过键值对存储数据的数据结构，可以快速地通过键来访问数据。 6. **树**（Tree）：一种层次结构的数据结构，每个节点有零个或多个子节点，通常用于表示具有层次关系的数据。 7. **图**（Graph）：由顶点（节点）和边组成，可以表示复杂的关系和网络结构。 每种数据结构都有其
weixin030英语学习交流平台小程序+ssm--pf.zip
08-19
随着互联网技术的快速发展，网络时代的到来，网络信息也将会改变当今社会。各行各业在日常企业经营管理等方面也在慢慢的向规范化和网络化趋势汇合。英语学习交流平台小程序的信息化程度体现在将互联网与信息技术应用于经营与管理，以现代化工具代替传统手工作业。无疑，使用网络信息化管理使信息管理更先进、更高效、更科学，信息交流更迅速。 对于之前英语学习交流平台小程序的管理，大部分都是使用传统的人工方式去管理，这样导致了管理效率低下、出错频率高。而且，时间一长的话，积累下来的数据信息不容易保存，对于查询、更新还有维护会带来不少问题。对于数据交接也存在很大的隐患。如果采用电子化的存储方式就会带来很大的改善，而且给用户的查询带来了很大便利，因此设计一个英语学习交流平台小程序刻不容缓，能够提高信息的管理水平。
Dynamic Test Loading
08-19
Dynamic Test Loading
基于 Spring Boot+AIGC的智能数据分析平台 区别于传统 BI，用户只需要导入原始数据集、 并输入分析诉求，就能自动
08-19
基于 Spring Boot+AIGC的智能数据分析平台。区别于传统 BI，用户只需要导入原始数据集、 并输入分析诉求，就能自动生成可视化图表及分析结论，
基于微信小程序的新生报到系统--.zip
08-19
本网站所实现的是一个新生报到系统，该系统严格按照需求分析制作相关模块，并利用所学知识尽力完成，但是本人由于学识浅薄，无法真正做到让该程序可以投入市场使用，仅仅简单实现部分功能，希望日后还能改善。 本系统具有以下优点： 该系统具有较高的适用性，选用B/S结构，可以在绝大部分个人平台上使用该系统。 系统将用户权限进行划分，管理员和学生能看到及操作的信息不一样，两者具备不同的操作权限。 该系统操作界面简单明了，大部分人都可以正常使用。
#-ssm-016-mysql-图书仓储管理系统-(1实用1).zip
08-19
随着社会经济的迅速发展和科学技术的全面进步，计算机事业的飞速发展，以计算机与通信技术为基础的信息系统正处于蓬勃发展的时期，随着经济文化水平的显著提高，人们对生活质量及工作环境的要求也越来越高，书籍做为人类的精神食粮，在现代社会中越来越受到重视，大量的书籍出现在市场上，人们有了各种各样不同的选择，与此同时，为了管理大量的图书，仓储的管理问题也就提上了日程，随着图书的大量增加，其管理难度也越来越大，如何优化仓储的日常管理也就成为了一个大众化的课题。因此，我选择开发这款图书仓储管理系统。 根据需求，确定系统采用JSP技术，SSM框架，JAVA作为编程语言，MySQL作为数据库。整个系统要操作方便、易于维护、灵活实用。主要实现了人员管理、库位管理、图书管理、图书报废管理、图书退回管理等功能。 关键词：图书仓储管理系统； JSP; MYSQL
JAVA WEB jstl资源jar包 包括jstl.jar和standard.jar
08-19
当然可以帮助您撰写关于`jstl.jar`和`standard.jar`这两个JAR包的介绍。以下是一段示例文本： --- ### JAR包资源介绍 #### 1. jstl.jar `jstl.jar`（JavaServer Pages Standard Tag Library，简称 JSTL）是一个用于简化 JavaServer Pages (JSP) 开发的标准标签库。它提供了一组功能强大的标签，这些标签可以用来控制流程、格式化输出、国际化处理、SQL操作等。使用JSTL可以显著减少页面中的Java代码量，提高代码的可读性和可维护性。 - **主要特性**： - **核心标签库**：提供了条件判断、循环、URL重定向等功能。 - **函数标签库**：支持日期、数字、字符串等数据类型的格式化。 - **国际化标签库**：支持多语言和本地化的页面展示。 - **XML标签库**：处理XML文档。 - **SQL标签库**：执行数据库操作。 - **适用场景**： - 适用于任何需要使用JSP技术进行Web开发的应用程序。 - **兼容性**：
基于Android学生课程表软件APP设计与实现.docx
08-19
基于Android学生课程表软件APP设计与实现.docx
计算机毕业设计仿慕课视频系统hsgAHB论文
08-19
计算机毕业设计仿慕课视频系统hsgAHB论文
物联网开发 物联网开发新手教程 包含源码和说明
08-19
物联网开发 物联网开发新手教程 包含源码和说明 【项目资源】：包含前端、后端、移动开发、操作系统、人工智能、物联网、信息化管理、数据库、硬件开发、大数据、课程资源、音视频、网站开发等各种技术项目的源码。包括STM32、ESP8266、PHP、QT、Linux、iOS、C++、Java、python、web、C#、EDA、proteus、RTOS等项目的源码。 【项目质量】：所有源码都经过严格测试，可以直接运行。功能在确认正常工作后才上传。 【适用人群】：适用于希望学习不同技术领域的小白或进阶学习者。可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。 【附加价值】：项目具有较高的学习借鉴价值，也可直接拿来修改复刻。对于有一定基础或热衷于研究的人来说，可以在这些基础代码上进行修改和扩展，实现其他功能。 【沟通交流】：有任何使用上的问题，欢迎随时与博主沟通，博主会及时解答。鼓励下载和使用，并欢迎大家互相学习，共同进步。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值