【洛谷】P4202 [NOI2008]奥运物流
【洛谷】P4202 [NOI2008]奥运物流
感觉有点降智
首先设环长为\(len\),很容易推导出
\]
我好像觉得似乎\(len\)越大底下越小了,太降智了,len越大底下越大
所以环长越小越好,所以我们如果动了一条环边,就是把它尽量往1连
枚举环上,强制要求环上上的点不选,除开1到1的后继的那条边,我们就获得了一棵树
我们希望\(dp\)出改了m个点最大是多少,显然改就一定直接连到1
设\(f[u][j][d]\)为\(u\)这个点,子树里改了\(j\)个点,我们硬点\(u\)到\(1\)的距离是\(d\)
转移的话分\(u\)有没有被连到1
如果被连到1了,那么就从\(max(f[v][k][1],f[v][k][2])\)转移过来,要加上\(C_{u} \times K\)
否则就从\(max(f[v][k][d + 1],f[v][k][1])\)转移过来,要加上\(C_{u} \times K^{d}\)
第二维是个背包,直接做就行了
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define ba 47
#define MAXN 5005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
int N,M;
int fa[65];
db K[65],C[65],f[65][65][65],g[65];
bool nouse[65];
void dfs(int u,int dep) {
for(int v = 2 ; v <= N ; ++v) {
if(fa[v] == u) dfs(v,dep + 1);
}
for(int d = dep ; d >= min(dep,2) ; --d) {
memset(g,0,sizeof(g));
for(int v = 1 ; v <= N ; ++v) {
if(fa[v] == u) {
for(int j = M ; j >= 0 ; --j) {
for(int k = j ; k >= 0 ; --k) {
g[j] = max(g[j],g[k] + max(f[v][j - k][d + 1],f[v][j - k][1]));
}
}
}
}
for(int j = 0 ; j <= M ; ++j) f[u][j][d] = g[j] + C[u] * K[d];
}
if(!nouse[u] && fa[u] != 1) {
memset(g,0,sizeof(g));
for(int v = 1 ; v <= N ; ++v) {
if(fa[v] == u) {
for(int j = M ; j >= 0 ; --j) {
for(int k = j; k >= 0 ; --k) {
g[j] = max(g[j],g[k] + max(f[v][j - k][2],f[v][j - k][1]));
}
}
}
}
for(int j = 0 ; j < M ; ++j) f[u][j + 1][1] = g[j] + C[u] * K[1];
}
}
void Init() {
read(N);read(M);scanf("%lf",&K[1]);
for(int i = 1 ; i <= N ; ++i) {
read(fa[i]);
}
K[0] = 1;
for(int i = 2 ; i <= N ; ++i) {
K[i] = K[i - 1] * K[1];
}
for(int i = 1 ; i <= N ; ++i) {
scanf("%lf",&C[i]);
}
}
void Solve() {
int p = 1;
db ans = 0;
for(int len = 2 ; len ; ++len) {
p = fa[p];
if(p == 1) break;
memset(nouse,0,sizeof(nouse));
int t = 1,l = 0;
while(l != len) {
nouse[t] = 1;
t = fa[t];
++l;
}
int rec = 0;
if(fa[p] != 1) {--M;rec = fa[p];fa[p] = 1;}
memset(f,0,sizeof(f));
dfs(1,0);
db res = f[1][M][0] / (1 - K[len]);
ans = max(ans,res);
if(rec) {++M;fa[p] = rec;}
}
printf("%.2lf\n",ans);
}
int main(){
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
Solve();
}
【洛谷】P4202 [NOI2008]奥运物流的更多相关文章
- bzoj 1065: [NOI2008] 奥运物流
1065: [NOI2008] 奥运物流 Description 2008北京奥运会即将开幕,举国上下都在为这一盛事做好准备.为了高效率.成功地举办奥运会,对物流系统 进行规划是必不可少的.物流系统由 ...
- [NOI2008]奥运物流
题目 洛谷 BZOJ 做法 单环有向图毒瘤题 不考虑环和改变后继:\(\sum\limits{i=1}^n C_i\cdot K^{dep(i)}\) 考虑环无穷等比求极m:\(R(1)=\sum\l ...
- 1065: [NOI2008]奥运物流 - BZOJ
Sample Input4 1 0.52 3 1 310.0 10.0 10.0 10.0Sample Output30.00 推荐题解:http://blog.csdn.net/whjpji/art ...
- 洛谷P3980 [NOI2008]志愿者招募
题解 最小费用最大流 每一天是一条边\((inf-a[i], 0)\) 然后对于一类志愿者, 区间两端连一条\((inf, c[i])\) \(S\)向第一个点连\((inf, 0)\) 最后一个点向 ...
- 【洛谷 P1772】 [ZJOI2006]物流运输(Spfa,dp)
题目链接 \(g[i][j]\)表示不走在\(i\text{~}j\)时间段中会关闭的港口(哪怕只关\(1\)天)从\(1\)到\(m\)的最短路. \(f[i]\)表示前\(i\)天的最小花费.于是 ...
- 洛谷 P1477 [NOI2008]假面舞会
题目链接 题目描述 一年一度的假面舞会又开始了,栋栋也兴致勃勃的参加了今年的舞会. 今年的面具都是主办方特别定制的.每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具.每个面具都有一个编号,主办方 ...
- 【BZOJ1065】【NOI2008】奥运物流(动态规划)
[BZOJ1065][NOI2008]奥运物流(动态规划) 题面 BZOJ 洛谷 题解 先不考虑环的情况,于是变成了一棵树. 这样子我们答案的贡献是\(\sum_{i=1}^nC_i\times k^ ...
- 【洛谷】P3980 [NOI2008]志愿者招募
[洛谷]P3980 [NOI2008]志愿者招募 我居然现在才会用费用流解线性规划-- 当然这里解决的一类问题比较特殊 以式子作为点,变量作为边,然后要求就是变量在不同的式子里出现了两次,系数一次为+ ...
- 洛谷 P1656 炸铁路
P1656 炸铁路 题目提供者kkksc03 标签图论搜索/枚举洛谷原创 难度普及/提高- 题目描述 因为某国被某红色政权残酷的高压暴力统治.美国派出将军uim,对该国进行战略性措施,以解救涂炭的生灵 ...
随机推荐
- [Luogu] 让我们异或吧
https://www.luogu.org/problemnew/show/P2420 异或满足 A ^ B = B ^ A A ^ A = 0 0 ^ A = A #include <cstd ...
- 「LOJ 121」「离线可过」动态图连通性「按时间分治 」「并查集」
题意 你要维护一张\(n\)个点的无向简单图.你被要求执行\(m\)条操作,加入删除一条边及查询两个点是否连通. 0:加入一条边.保证它不存在. 1:删除一条边.保证它存在. 2:查询两个点是否联通. ...
- Day1 读题解题提升
The 2014 ACM-ICPC Asia Mudanjiang Regional Contest 昨晚做了训练赛,然后读题又自闭了QAQ. Average Score ZOJ - 3819 题意: ...
- 内存管理4-Autoreleasepool
自动释放池是OC里面的一种内存回收机制,一般可以将一些临时变量添加到自动释放池中,统一回收释放,当自动释放池销毁时,池里面的所有对象都会调用一次release,也就是计数器会减1,但是自动释放池被销毁 ...
- JavaWeb_(Spring框架)整合Mybatis加入事务操作数据库
整合Mybatis a)导包: i.Spring:基本包.aop.aspects.jdbc.tx.test: ii.Mybatis:mybatis-3.4.6 iii.整合包:mybatis-spri ...
- Kafka - SASL认证
kafka SASL认证配置 1.找到kafka安装根目录,在config文件夹下创建kafka_server_jaas.conf,写入 KafkaServer { org.apache.kafka. ...
- Flutter点击两次返回键退出APP
在APP中一些页面为了防止用户操作失误点击到返回键导致退出APP,可以设置其一定时间内点击两次返回键才允许退出APP,完成这个功能可以通过WillPopScope和SystemNavigator.po ...
- 黑马vue---61、为什么vue组件的data要是一个函数
黑马vue---61.为什么vue组件的data要是一个函数 一.总结 一句话总结: 因为js中以函数为变量作用域,所以这样可以保证每个组件的数据不互相影响 二.why components data ...
- 【matlab】模拟变焦拼接代码备份
1.初版,边缘未处理. % % In----near % If----far % In=imread('D:\文件及下载相关\桌面\模拟变焦拼接\Matlab_code\nearframe\frame ...
- Cloudera Manager 概念
cloudera公司发布的CDH集群,使用Cloudera Manager来管理整个集群,使用过程中主要涉及到几个关键概念:cluster.service.role.host.直接上图,直观理解几个概 ...