sepNINE的专栏

//poj 1485//fenix_bao#include <iostream>using namespace std;const int MAXN=256,MAXK=32;int pos[MAXN],dp[MAXN][MAXK],dis[MAXN][MAXN],one_cover[MAXN][MAXN];int sum_dis[MAXN][MAXN],sum_dis_r...

//poj 1015//sep9#includeusing namespace std;const int delta=400;int a[202],b[202];int dp[202][22][2*delta+10];int path[202][22][2*delta+10];int n,m;void print(int x,int y,int z){ if(x==0)

//poj 2677//sep9#include #include using namespace std;const int maxN=128;struct P{ double x,y;}p[maxN];int n;double dp[maxN][maxN],dis[maxN][maxN];double solve(){ dp[1][2]=dis[1][2]; fo

//poj 2625//sep9#include using namespace std;typedef __int64 INT;INT n,up,down;INT gcd(INT x,INT y){ return x%y==0?y:gcd(y,x%y);}void add(INT x,INT y){ if(down==0){ up=x; down=y; re

//poj 2564//sep9#include using namespace std;const int MAXN=25000;int n,idx;char tmp[20],s[MAXN+10][20];int dp[MAXN+10];struct TRIE{ int c[27]; bool vis; int ids;}a[8*MAXN+10];void inse

题意：给二叉搜索树的内部节点p1,p2...pn和外部节点q0,q1,..qn,求一棵比较数期望最小的二叉搜索树，输出期望。分析：n^3的动态规划很容易想到，方程满足四边形优化的条件，可优化的n^2.代码（n^3）://poj 1784//sep9#include using namespace std;const int maxN=256;int dp[maxN][

//poj 2250//sep9#include using namespace std;struct WORD{ char s[32];}a[128],b[128];int dp[128][128],match[128][128];void print(int x,int y){ if(x==0||y==0) return ; if(match[x][y]==0){

//poj 2288//sep9#include using namespace std;typedef __int64 INT;INT v[16];INT dp[(1<<13)+10][13][13],path[(1<<13)+10][13][13];bool map[16][16];int n,m;INT rec(int s,int now,int pre){ if(d

//poj 2353//sep9#includeusing namespace std;int dp[128][512],cost[128][512],path[128][512];int m,n;void print(int row,int col){ if(path[row][col]==col) print(row-1,col); else if(path[row][c

//poj 3181//sep9#include using namespace std;struct H { int a[10],len; void clear(){ memset(a,0,sizeof(a));len=1; } void print(){ printf("%d",a[len]

题意：良心中文题不赘述了。分析：状态压缩dp,每行的状态和前两行有关。代码：//poj 1185//sep9#include using namespace std;int n,m;char g[128][16];int row_state[128][128];int dp[128][64][64];bool check(int row,int st){ f

?/poj 2404//sep9#include using namespace std;const int MAX=INT_MAX/2;const int MAXN=16;int n,m;int deg[MAXN],d[MAXN][MAXN],dp[1<<MAXN];int rec(int s){ if(s==0) return dp[s]=0; if(dp[s]>

//poj 3272//sep9#include using namespace std;const int MAXN=5012;const int MAXM=50012;int n,m,e,rev_e;int head[MAXN],rev_head[MAXN],flag[MAXN],dp[MAXN],rev_dp[MAXN];struct EDGE{ int v,next;}

题意：poj的良心中文题，不赘述了。分析：难点是想到能把题中的最优单连通子集问题转化为求树的最大权子树问题，之后我用的是孩子兄弟表示，当然在建树的时候直接dp也行。树的最大权子树感觉是最简单的树型dp了。代码：//poj 1192//sep9#include #include #include using namespace std;const int maxN=1

题意：给在商店里要买的商品的种类，数量，价格，和一些组合购买选项（如同时买2个a，1个b，只要x元），求最小花费。分析：因为要买的商品种类数b是变量（b代码：//poj 1170//sep9#include using namespace std;int b,c,k,p,s;int dp[8][8][8][8][8];int tot_idx[8];int idx[

//poj 1404//sep9#includeusing namespace std;int k,l;int frequency[128];int cost[128][128],price[128][128],idx[128][128];char keys[128],letters[128];void ini(){ memset(price,0x40,sizeof(pric

题意：给添加和删除每个字符的代价，求将一个字符串变成回文串的最小代价。分析：区间dp，采用了备忘录方法。代码：//poj 3280//sep9#includeusing namespace std;const int maxN=2048;int n,m;int dp[maxN][maxN];char s[maxN];int add_cost[256],del_c

//poj 3790//sep9#includeusing namespace std;int dp[1024];int rec(int n){ if(n==0) return 1; if(dp[n]!=-1) return dp[n]; int sum=0; if(n%2){ for(int i=1;i<=n;i+=2) sum+=rec((n-i)/2);

//poj 3797//sep9#include using namespace std;int dp[2][128];int n,m;int solve(){ if(m>n) swap(m,n); memset(dp,0,sizeof(dp)); int *crt=dp[0],*nxt=dp[1]; crt[0]=1; for(int i=n-1;i>=0;--i)

//poj 1848//sep9include #include using namespace std;const int MAXN=128;const int INF=100000;vector g[2*MAXN];int dp[MAXN][3];void dfs(int u,int p){ vector son; int sum=0; for(int i=0;i<

题意：       设dp[i]=min dp[j]+(S+sumT[i]-sumT[j])*F[i]  1分析：       dp[i]=min (S+sumT[i])*sumF[i]+dp[j]-sumT[j]*F[i],由于F[i]随i减小单增，sumT[j]随i减小单增，按照i从N到1倒推dp[i]满足单调队列优化的条件dp[i]=S(i)+min{fj(i)|i代码：

//poj 1322 AC//sep9#include using namespace std;int c,n,m;double dp[2][128];int main(){ while(scanf("%d",&c)==1&&c){ scanf("%d%d",&n,&m); if(m>c||m>n||(m+n)%2){ puts("0.000"); conti

//poj 1505//sep9#include using namespace std;const int MAXN=512;int sum[MAXN],path[MAXN],dp[MAXN][MAXN];int main(){ int cases; scanf("%d",&cases); while(cases--){ int m,k; scanf("%d%d",

//poj 1189//sep9#include using namespace std;typedef long long ll;int n,m;char g[64][64];ll dp[64][64];ll gcd(ll a,ll b){ return a%b==0?b:gcd(b,a%b);}void print(ll a,ll b){ ll d=gcd(a,b

//poj 2151//sep9#include using namespace std;int n,m,t;double p[1024][35],dp[1024][35][35],s[1024][35];int main(){ while(scanf("%d%d%d",&m,&t,&n)==3){ if(!(n+m+t)) break; for(int i=1;i<=t;

题意：给n个用左端点、右端点。高度标志的平台，问一个下落速度，移动速度均为1的点下落到地面的最小时间，其中每次下落距离不能超过输入值MAX.分析：动态规划ldp[i]表示从第i块木板左边下落到地面的最小时间，rdp[i]表示从第i块木板右边下落到地面的最小时间。代码：//poj 1661//sep9#include #include #include using na

//poj 1036//sep9//这种dp效率低 会爆内存#include using namespace std;const int MAXT=30012;const int MAXS=112;const int MAXN=112;int n,k,t;int a[MAXT][MAXS];int dp[MAXT][MAXS];int x[MAXN],y[MAXN],z[MA

//poj 1038#include using namespace std;int n,m,k,x,y,ans,o;int map[180][16],f[2][60000],three[16],p[16],q[16];void init(){ three[0]=1; for(int i=1;i<=11;++i) three[i]=three[i-1]*3; }int

//poj 1037//sep9#include using namespace std;typedef long long ll;ll dp[32][32][2];bool vis[32];int main(){ memset(dp,0,sizeof(dp)); dp[1][1][0]=dp[1][1][1]=1; for(int len=2;len<=20;++len)

//poj 3356//sep9#include using namespace std;const int MAXN=1024;char a[MAXN],b[MAXN];int dp[MAXN][MAXN];int main(){ int lenA,lenB; while(scanf("%d%s%d%s",&lenA,&a[1],&lenB,&b[1])!=EOF){

//poj 1082//sep9#include using namespace std;int dp[200][14][33];int month_day(int y,int m){ if(y%4==0&&y%100!=0||y%400==0) if(m==2) return 29; if(m==1||m==3||m==5||m==7||m==8||m==10||m=

//poj 1252//sep9#include using namespace std;const int MAXV=1200;int val[8];int dp[MAXV+10];int main(){ int cases; scanf("%d",&cases); while(cases--){ for(int i=0;i<6;++i) scanf("%d",

//poj 1155//sep9#include #include using namespace std;const int MAXN=3012;int n,m,ans;struct Node{ int v,w;};vector g[MAXN];int num[MAXN],dp[MAXN][MAXN];void dfs(int u){ for(int i=g[u].

题意：给字符串a,b,求长度最短的包含a且包含b的字符串。分析：是个LCS的对偶问题，dp的方程也很像。代码：//poj 2264//sep9#include using namespace std;const int MAXN=128;int dp[MAXN][MAXN],path[MAXN][MAXN];char a[MAXN],b[MAXN];void pr

题意：等价于在一棵树中，找m个叶子节点，使他们连线（重复部分只计一次）的距离和最大。分析：m==2时求树的直径倒是简单，m>2时就要用树形dp了，有n*m^2的树形dp.代码：//poj 1849//sep9#include #include using namespace std;const int MAXN=100024;vector g[MAXN],val[M

//poj 1283//sep9#include using namespace std;typedef __int64 INT;const int MAXN=200;INT dp[MAXN+10][MAXN+10];int main(){ int N,K; memset(dp,0,sizeof(dp)); for(int i=0;i<=MAXN;++i) dp[i]

代码：//poj 2096//sep9#include #include using namespace std;double E[1005][1005];int main(){ int n,s; while(scanf("%d%d",&n,&s)==2){ for(int i=0;i<=n+1;++i) for(int j=0;j<=s+1;++j)

题意：给一个迷宫，求左下角到右下角的路径数。分析：插头dp的模板题，建议先看cdq的论文再看代码，这份代码在模板基础上略微有改动。论文地址http://wenku.baidu.com/view/ed2b3e23482fb4daa58d4b74.html代码：#include using namespace std;const int maxD=16;const int H

题意：有2*n只队伍参加n论的淘汰赛，给出任意两只队伍交战的各自的胜率，求最后剩下概率最高的队伍。分析：概率dp,dp[i][j]表示第i轮结束后队伍j还在场上的概率，dp[i][j]=dp[i-1][j]*sum(dp[i-1][k]p[j][k])(1>i==(k-1)>>i&&(j-1)>>(i-1)!=(k-1)>>(i-1)，把每个队的id用二进制写出来画画图就知道为什么这样

dp水题，直接贴代码。//poj 3670//sep9#includeusing namespace std;const int maxN=30024; int n;int a[maxN],b[maxN];int dp[maxN][4];int solve(int a[]){ int i,j; for(i=0;i<=n;++i) for(j=1;j<=3;++j)