DP思路
在这里记录一些在大神们的博客,以及自己做过的一些DP的神奇思路吧
1.2015/04 NEUQ 月赛 转自:http://zyfzyf.is-programmer.com/posts/89993.html
E.又被DP卡住了。 感觉是必须记录和的具体差值的。因为只有最值无法保证子问题最优的性质。 当然如果既记录具体差值,又记录交换了多少次的话不仅MLE而且TLE。 然后我就弃疗了,觉得一定是防AK神题! 然后学姐发了题解,顿时感觉好巧妙啊。 我们只记录上面的和。用f[i][j]表示前i个数组成和为j最少需要交换多少次,这个满足子问题最优的性质,而且可以很简单的DP。 太神了,Orz! 这次是败在了状态的设计上,以后要多下功夫。(好像脑袋里闪过这个思路?但立马被否决了?) 因为E比较好,所以搬运一下题面:
陈船长和xzx每人都有N个玩具,序号从1到N。不同的玩具可能并不一样好玩。
其中陈船长的第i个玩具的好玩度为C[i],xzx的第i个玩具的好玩度为X[i], 好玩度可以是负数,可以理 解为不好玩的程度。
你可以交换陈船长的第i个玩具和xzx的第i个玩具,为了让他们两个玩的尽可能一样高兴,他们希望经过 不大于K次交换后陈船长所有玩具的好玩度总和与xzx的总和之差的绝对值最小,即 |Sum(C) - Sum(X)| 最小。现在,请你告诉他俩最小值是多少。
多组测试数据,第一行一个整数T(T < ), 表示测试数据组数。
每组数据第一行两个整数N( < N < ), K( < K <= N)。
接下来一行N个整数,第i个数代表C[i](- <= C[i] <= )。
再接下来一行N个整数, 第i个数代表X[i](- <= X[i] <= )。
数组下标从0开始。 比赛地址
2. 【模拟试题1】【20150514】
三种物品的背包:
1. $v(x)=A*x^2-B*x$ 价值随所分配的体积的变化而变化……
2. 多重背包
3. 完全背包
其实是个傻逼题,因为数据规模小,暴力就能过,然而由于没见过第一种物品的价值函数,加上题目描述不清楚,所以自己傻逼了(不知道怎么乱搞了一下,得了70分)
其实第一种物品就暴力枚举每一个物品我们分配给它多少空间就可以,“每个体积的甲类物品只有一个”TM就是在误导人!
不过倒是学习了一种新的写背包的姿势~
//2015_05_14 pack
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
int r=,v=; char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-;
for(; isdigit(ch);ch=getchar()) v=v*-''+ch;
return r*v;
}
const int N=,M=;
/*******************template********************/
int n,m;
LL f[M]; void solve1(int a,int b){
D(j,m,) F(i,,j)
f[j]=max(f[j],f[j-i]+a*i*i-b*i);
}
void zeroone(int a,int b){
D(j,m,b)
f[j]=max(f[j],f[j-b]+a);
}
void solve2(int a,int b,int c){
int k=;
while(c>=(<<k)){
zeroone(a*(<<k),b*(<<k));
c-=(<<k);
k++;
}
if (c) zeroone(a*c,b*c);
}
void solve3(int a,int b){
F(j,b,m)
f[j]=max(f[j],f[j-b]+a);
}
int main(){
#ifndef ONLINE_JUDGE
freopen("pack.in","r",stdin);
freopen("pack.out","w",stdout);
#endif
n=getint(); m=getint();
F(i,,n){
int x=getint(), A=getint(), B=getint(),C;
if (x==){
solve1(A,B);
}else if (x==){
C=getint();
solve2(A,B,C);
}else{
solve3(A,B);
}
}
printf("%d\n",f[m]);
return ;
}
3. 【BestCoder】【Round#41】 & 【BZOJ】【3612】【HEOI 2014】平衡
整数拆分的DP姿势= =
f[i][j]表示将 i 拆成 j 个不同的数(均不超过n)的方案数
f[i][j]=f[i-j][j]+f[i-j][j-1]
if (i>n+1) f[i][j]-=f[i-(n+1)][j-1](相当于去掉了最后那一列/一个数(n+1))
4.【BZOJ】【1089】【SCOI2003】严格n元树
f[i]=f[i-1]^n+1 f[i]表示深度小于等于 i 的严格n元树的数目
ans=f[d]-f[d-1] //这个表达式感觉好神……
5.【BZOJ】【2750】【HAOI2012】Road&& CTSC D1T1 好朋友
拓扑序DP,然而2750这题的DP是较简单的一个,统计每条边在多少条最短路径上……其实可以统计 有多少条最短路径经过了x,以及y出发到达任意一个结束点有多少种走法(沿最短路)
我们可以用Dijkstra求出以 i 为起点的最短路径图,它是一个DAG,然后我们用dij扩展的顺序(一个拓扑序)来搞DP!
令a[x]表示从 i 沿最短路走到 x 的方案数,b[x]表示从 x 往出走,沿最短路走到任意一个结束结点的总方案(就是经过x的最短路条数?只不过只看后半段)
而CTSC的我还没搞出来……
UPD:
答案表达式是这样的:$$ ans[v]=\sum_{s \not = v \not =t} \frac{a[s]*a[t]*\sigma_{s,t}(v)}{\sigma_{s,t}} $$
f1[x]表示从S到x的最短路宽度,然后我们需要处理的是最短路的后一半。。
如果没有分母那个东西,f2[x]应该是 $\sum (f2[to[i]]+a[to[i]])*width[i] $
然而分母和分子是不可以分开搞的……
但是分母其实我们已经算出来了!!就是f1[x]!(我是sb一直没想到……)
所以维护f2[x]还是很容易……$f2[x]=\sum (f2[to[i]]+\frac{a[to[i]]}{f1[to[i]]})*width[i]$
sad....我需要吃药。。。脑残片QAQ
/**************************************************************
Problem: 4055
User: Tunix
Language: C++
Result: Accepted
Time:6976 ms
Memory:1520 kb
****************************************************************/ #include<queue>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
int r=,v=; char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-;
for(; isdigit(ch);ch=getchar()) v=v*-''+ch;
return r*v;
}
//#define debug
const int N=,M=,INF=;
/*******************template********************/ typedef long double db;
int head[N],next[M<<],to[M<<],len[M<<],cnt;
db b[M<<];
void ins(int x,int y,int z,db wd){
to[++cnt]=y; next[cnt]=head[x]; head[x]=cnt; len[cnt]=z; b[cnt]=wd;
}
void add(int x,int y,int z,db wd){
ins(x,y,z,wd); ins(y,x,z,wd);
}
int n,m,a[N];
typedef pair<int,int> pii;
#define mp make_pair
priority_queue<pii,vector<pii>,greater<pii> >Q;
int c[N],d[N],vis[N];
db ans[N],f1[N],f2[N]; void Dij(int S){
F(i,,n) d[i]=INF,vis[i]=,c[i]=;
d[S]=; Q.push(mp(,S));
int t=;
while(!Q.empty()){
int x=Q.top().second; Q.pop();
if (vis[x]) continue;
vis[x]=++t;
for(int i=head[x];i;i=next[i])
if (d[to[i]]>d[x]+len[i]){
d[to[i]]=d[x]+len[i];
Q.push(mp(d[to[i]],to[i]));
}
}
F(i,,n) f1[i]=f2[i]=; f1[S]=;
F(i,,n) c[vis[i]]=i;
// F(i,1,t) printf("%d ",c[i]); puts("");
F(i,,t){
int x=c[i];
for(int j=head[x];j;j=next[j])
if (d[to[j]]==d[x]+len[j]) f1[to[j]]+=f1[x]*b[j];
}
D(i,t,){
int x=c[i];
for(int j=head[x];j;j=next[j])
if (d[to[j]]==d[x]+len[j])
f2[x]+=b[j]*(f2[to[j]]+(db)a[to[j]]/(db)f1[to[j]]);
}
F(i,,n) if (i!=S) ans[i]+=(db)a[S]*f1[i]*f2[i];
}
int main(){
#ifndef ONLINE_JUDGE
freopen("misc.in","r",stdin);
freopen("misc.out","w",stdout);
#endif
n=getint(); m=getint();
F(i,,n) a[i]=getint();
F(i,,m){
int x=getint(),y=getint(),z=getint();
double t1; scanf("%lf",&t1);
add(x,y,z,t1);
}
F(i,,n)
Dij(i);
F(i,,n){
double x=ans[i];
printf("%.9f\n",x);
}
return ;
}
(CTSC D1T1)
6.树形DP【模拟试题3】【20150527】
(1)树形状压DP,其实在看题解之前我似乎并没有搞懂这题在干什么……
对于节点 i ,我们考虑f[i][j]表示 i 这棵子树中,分部包含情况为 j 的最大收益,因为一个分部管的是从x到根的所有点,所以对于一个点来说,管它的就是子树中的所有分部,所以就可以dp啦~依次枚举每个儿子的子树中有哪些分部,用一个辅助数组,我们可以搞:
c[j]=f[x][j];
c[j]=max(c[j],f[y][k]+f[x][j^k]);
f[x][j]=c[j];
依次搞下来即可。
初始化就是分部全部在 x 节点。计算完后要加上符合的收益(T种中的某一些……)
(2)无向图最大权路径覆盖= =?然而其实并不是费用流之类……
f[x][0]表示x是孤立点时,x子树内最大收益;
f[x][1]表示x子树内有一条链与它相连;
f[x][2]表示x子树内有两条链与它相连
前两种情况是可以继续向上连的,而第三种情况就不可以了……
转移……Orz regina8023 其实是个贪心。
先假设所有点都没跟x连,得到f[x][0],然后再考虑一下将某个儿子连到父亲,会使答案改变多少,贪心地选出1-2个,即可得到f[x][1]和f[x][2]。
感觉这个先假设全选、【求差】然后在差值中选最优的思路好神啊。。。(其实是我太弱了没见过
7.【BZOJ】【2878】【NOI2012】迷失游乐园
期望&树形&基环树DP的好题QAQ,题解太长我就不复制了,感觉这个先求出来down再求up的方法好神奇……
8.【TYVJ 五月图论专项有奖比赛】
第三题:求树的双重心?。。。树形dp很神奇,Orz zyf。
枚举断掉的点,然后找最大的儿子往过走,更新答案。。。QwQ说不清了
9.【BZOJ】【4145】【AMPPZ2014】The Prices
状压->01背包
DP思路的更多相关文章
- Codeforces 1204D2. Kirk and a Binary String (hard version) (dp思路)
题目链接:http://codeforces.com/contest/1204/problem/D2 题目是给定一个01字符串,让你尽可能多地改变1变为0,但是要保证新的字符串,对任意的L,R使得Sl ...
- 洛谷P1244 青蛙过河 DP/思路
又是一道奇奇怪怪的DP(其实是思路题). 原文戳>>https://www.luogu.org/problem/show?pid=1244<< 这题的意思给的挺模糊,需要一定的 ...
- 很好的DP思路,字符串比较次数
题目: https://leetcode.com/problems/distinct-subsequences/?tab=Description 一般没有明显思路的情况下,都要想想DP,用下Divid ...
- BZOJ.3227.[SDOI2008]红黑树tree(树形DP 思路)
BZOJ orz MilkyWay天天做sxt! 首先可以树形DP:\(f[i][j][0/1]\)表示\(i\)个点的子树中,黑高度为\(j\),根节点为红/黑节点的最小红节点数(最大同理). 转移 ...
- 洛谷P2362 围栏木桩----dp思路
在翻dp水题的时候找到的有趣的题0v0 原文>>https://www.luogu.org/problem/show?pid=2362<< 题目描述 某农场有一个由按编号排列的 ...
- AGC033 D~F——[ 值放到角标的DP ][ 思路+DP ][ 思路 ]
地址:https://atcoder.jp/contests/agc033/ D Complexity dp[ i ][ j ][ k ][ l ] 表示左上角是 ( i , j ) .右下角是 ( ...
- [状压DP思路妙题]图
源自 luhong 大爷的 FJ 省冬令营模拟赛题 Statement 给定一个 \(n\) 个点 \(m\) 条边的图,没有重边与自环 每条边的两端点编号之差不超过 \(12\) 求选出一个非空点集 ...
- ZOJ 3647 Gao the Grid dp,思路,格中取同一行的三点,经典 难度:3
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4837 三角形的总数=格子中任取3个点的组合数-同一横行任取3个点数目-同一纵行 ...
- BZOJ.1566.[NOI2009]管道取珠(DP 思路)
BZOJ 洛谷 考虑\(a_i^2\)有什么意义:两个人分别操作原序列,使得得到的输出序列都为\(i\)的方案数.\(\sum a_i^2\)就是两人得到的输出序列相同的方案数. \(f[i][j][ ...
随机推荐
- 洛谷P4742 [Wind Festival]Running In The Sky [Tarjan缩点,DAGDP]
题目传送门 Running In The Sky 格式难调,题面就不放了. 分析: 一句话题意:给定一张带点权的有向图,求最长点权路径及该路径上的最大点权. 很明显的$DAGDP$,因此需要缩点,将该 ...
- q.size()
在队列中,q.empty()比q.size()要慢,如果可能就用q.size(); 比如spfa的 while(q.size()>0) better than while(!q.empty())
- 010.Zabbix的zatree插件安装
一 zatree简介 zatree 是来自国内58公司开发的监控软件zabbix的一个插件,主要功能是提供host group的树形展示和在item里指定关键字查询及数据排序. 二 安装前准备 2.1 ...
- Javascript的字符串(String)操作学习
1.bold() 方法用于把字符串显示为粗体.语法: stringObject.bold() 如下,对str进行bold操作之后,实际上时对这个字符串加了<b>标签,在文档中将以粗体进行展 ...
- 循序渐进学.Net Core Web Api开发系列【14】:异常处理
系列目录 循序渐进学.Net Core Web Api开发系列目录 本系列涉及到的源码下载地址:https://github.com/seabluescn/Blog_WebApi 一.概述 本篇介绍异 ...
- 设计模式-装饰者模式(Decorator Pattern)
本文由@呆代待殆原创,转载请注明出处. 此设计模式遵循的设计原则之一:类应该支持扩展,而拒绝修改(Open-Closed Principle) 装饰者模式简述 装饰者模式通过组合的方式扩展对象的特性, ...
- centos7 静默安装oracle
系统centos7.4 mini 关闭selinux.firewalld 配置主机名: hostnamectl set-hostname --static oracle 之前说oracle不认cen ...
- [BZOJ4129]Haruna’s Breakfast(树上带修改莫队)
BZOJ3585,BZOJ2120,BZOJ3757三合一. 对于树上路径问题,树链剖分难以处理的时候,就用树上带修改莫队. 这里的MEX问题,使用BZOJ3585的分块方法,平衡了时间复杂度. 剩下 ...
- 喵哈哈村的魔法考试 Round #1 (Div.2) 题解
喵哈哈村的魔法考试 Round #1 (Div.2) 题解 特别感谢出题人,qscqesze. 也特别感谢测题人Xiper和CS_LYJ1997. 没有他们的付出,就不会有这场比赛. A 喵哈哈村的魔 ...
- Codeforces Round #372 (Div. 1) A. Plus and Square Root 数学题
A. Plus and Square Root 题目连接: http://codeforces.com/contest/715/problem/A Description ZS the Coder i ...