THUPC2018 城市地铁规划
$n$ 个点,你可以随意连成一棵树,一个点的贡献为 $F(度数) \space mod \space 59393$ ,$F$ 为给定多项式函数,不超过 $10$ 次
求这 $n$ 个点的最大贡献,和最后连出来的那棵树
$n \leq 3000$
sol:
看到这种跟树度数有关的题大概是要上 prufer 序列?
对 prufer 序列进行 dp,每个点大概相当于一个物品,由于 prufer 序列可以任意放,大概还是个完全背包
于是可以写出一个朴素的转移式:
$f_{(i,j)}$ 表示 prufer 序列的前 $i$ 项,目前出现的最后一个数出现了 $j$ 次的最大贡献
每次可以转移到 $f_{(i+1,j+1)}$ (填一个一样的)或者 $f_{(i+1,1)}$ (填一个新的)
或者优秀一点(平时刷题啥都敢写系列),直接对于 prufer 序列上每一个位置,分配它连了多少个叶子,这样空间是一维的,也少了很多分类讨论
后面把 prufer 序列变成一棵树就...用个 set 模拟一下就完事了
#include<bits/stdc++.h>
#define LL long long
#define rep(i, s, t) for(register LL i = (s), i##end = (t); i <= i ## end; ++i)
#define dwn(i, s, t) for(register LL i = (s), i##end = (t); i >= i ## end; --i)
using namespace std;
inline int read() {
int x = , fv = ; char ch = getchar();
for(;!isdigit(ch);ch=getchar())if(ch == '-') fv=-fv;
for(;isdigit(ch);ch=getchar())x = * x + ch - '';
return x * fv;
}
const int maxn = , mod = ;
inline int inc(int x, int y) { x += y; if(x >= mod) x -= mod; return x; }
inline int dec(int x, int y) { x -= y; if(x < ) x += mod; return x; }
inline int mul(int x, int y) { return 1LL * x * y % mod; }
inline int power(int a, int b) {
int res = ;
for (; b; b >>= , a = mul(a, a))
if (b & ) res = mul(res, a);
return res;
}
int n, k, a[], fv[maxn], p[maxn];
LL f[maxn];
inline int getval(int x) {
int res = ;
for (int i = , j = ; i <= k; j = mul(j, x), i++) res = inc(res, mul(j, a[i]));
return res;
}
multiset<int> S, leaves;
int main() {
n = read(); k = read();
rep(i, , k) a[i] = read();
rep(i, , n) fv[i] = getval(i);
if (n == ) {
cout << << " " << a[] << endl;
return ;
}
if (n == ) {
cout << << " " << inc(fv[], fv[]) << endl << << " " << << endl;
return ;
}
rep(i, , n-) rep(j, , i)
if (f[i - j] + 1LL * fv[] * (j - ) + fv[j + ] > f[i])
f[i] = f[i - j] + 1LL * (j - ) * fv[] + fv[j + ], p[i] = i - j; printf("%d %lld\n", n-, f[n-] + fv[] + fv[]);
for (int i = n - , j = i, cnt = ; i; ++cnt, j = i = p[i])
while (j > p[i]) S.insert(cnt), --j;
rep(i, , n) if (!S.count(i)) leaves.insert(i);
rep(i, , n-) {
int u = *leaves.begin(); leaves.erase(u);
int v = *S.begin(); S.erase(S.find(v));
printf("%d %d\n", u, v);
if (!S.count(v)) leaves.insert(v);
}
printf("%d %d\n", *leaves.begin(), *leaves.rbegin());
}
THUPC2018 城市地铁规划的更多相关文章
- 【LibreOJ】#6395. 「THUPC2018」城市地铁规划 / City 背包DP+Prufer序
[题目]#6395. 「THUPC2018」城市地铁规划 / City [题意]给定n个点要求构造一棵树,每个点的价值是一个关于点度的k次多项式,系数均为给定的\(a_0,...a_k\),求最大价值 ...
- 大数据freestyle: 共享单车轨迹数据助力城市合理规划自行车道
编者按:近年来,异军突起的共享单车极大地解决了人们共同面临的“最后一公里”难题,然而,共享单车发展迅猛,自行车道建设却始终没有能够跟上脚步.幸运的是摩拜单车大量的轨迹数据为我们提供了一种新的思路:利用 ...
- ZROI 暑期高端峰会 A班 Day6 DP
[THUPC2018]城市地铁规划 (日常讲题之前 YY--) 一眼出 \(O(n^3+nk)\) 做法. \(dp[i][j]\) 表示前 \(i\) 个点,前 \(i\) 个点度数和为 \(j\) ...
- zhengrui集训笔记2
Day_6 计算几何 点积\Large 点积点积 叉积\Large 叉积叉积 极角\Large 极角极角 < π\piπ :叉积判断 else :atan2 旋转\Large 旋转旋转 左乘第一 ...
- 中科燕园GIS外包---地铁GIS项目
(1)地铁保护及project地质管理 • 地铁保护 地铁交通既有运量大,速度快的特点,又有差别于其它交通方式的在地下执行的空间特殊性,因此地铁的保护显得尤为重要. 首先必须编制完整的 ...
- SVG的动态之美-搜狗地铁图重构散记
搜狗地图发布了新版的移动端地铁图,改版初衷是为了用户交互体验的提升以及性能的改善.原版地铁图被用户吐槽最多的是pinch缩放不流畅.无过渡动画.拖拽边界不合理等等,大体上都是交互体验上的问题.实际上原 ...
- 地铁合作(hjy)
用时大概时间 18 个小时 我的队友是 :黄珺瑜我们一开始的想法是自己建立数据库,用来查询路线,再调用一下百度api作为地图(仅作地图没有任何操作的那种):后来我们的步骤也的确是这样,我们先确 ...
- 机器人与机器人仿真技术(zz)
http://www.viblue.com/archives/5587.htm 一.机器人简介: 机器人(Robot)是自动执行工作的机器装置.它既可以接受人类指挥,又可以运行预先编排的程序,也可以根 ...
- 2018年度最优秀mac软件及游戏推荐,个个万里挑一
今天和大家带来2018年度最优秀Mac软件和游戏合集,个个万里挑一,2018年,风云社区(scoee.com)分享了上数千款优秀的Mac软件和游戏,结合用户反馈,精选出各个类别的优秀的Mac软件,推荐 ...
随机推荐
- oracle导入数据库报错:IMP-00019: 由于 ORACLE 错误 12899 而拒绝行 IMP-00003: 遇到 ORACLE 错误 12899
主要是字符集 成 导入的Oracle服务器的字符集 对应不上. 以下方案为是修改服务器的字符集. 这样会影响之前的Oracle其他数据库的数据显示(正式服务器慎用) 个人认为应该修改导出文件的字 ...
- javaweb学习纲要
Java Web 概述 大纲: 1.C/S体系结构 2.B/S体系机构 3.两种体系结构比较 4.主流的Web程序应用平台 5.java web学习路线图 1.C/S体系结构 C/S是Client/ ...
- 每天一个Linux命令(52)telnet命令
执行telnet指令开启终端机阶段作业,并登入远端主机. (1)用法: 用法: telnet [参数] [主机] (2)功能: 功能: telnet命令通常 ...
- 设计模式—迭代器Iterator模式
什么是迭代器模式? 让用户通过特定的接口访问容器的数据,不需要了解容器内部的数据结构. 首先我们先模仿集合中ArrayList和LinkedList的实现.一个是基于数组的实现.一个是基于链表的实现, ...
- pppoe白皮书
转:https://blog.csdn.net/achejq/article/details/19478811 PPPoE技术白皮书 关键词:PPP,Ethernet,PPPoE 摘 要:PPP ...
- @MarkFan 口语练习录音 20140401
Hi,everybody 对于未来,我只梦想最好的情况, 并定下最踏实的计划,而绝不花时间在无谓的担心上, 因为我知道,只要把我对自己的承诺付诸实践, 我的未来将不会只是一个梦…… 这是引用考拉小巫的 ...
- nginx location proxy pass
nginx: 192.168.1.23作为nginx反向代理机器 目标机器192.168.1.5上部署一个8090端口的nginx [root@localhost conf.d]# cat test. ...
- java基础之常量与变量
概要:通过这段时间的工作,发现自己的基础还是很薄弱的,so,you know 常量 一种特殊的变量,程序运行过程中不能改变的值 语法格式:final 数据类型 常量名称 = 常量值 例子:fina i ...
- groupby和agg的使用
先来看一段代码: 分析下groupby和agg的联合使用: reset_index()表示重新设置索引 agg传进来的统计特征: 按照A这一列作聚合,C这一列作统计 注意:df = df.groupb ...
- C语言下文件目录查看
C语言下文件目录遍历通常会用到下面这些函数 _access() /* 判断文件或文件夹路径是否合法 */ _chdir() /* 切换当前工作目录 */ _findfirst() / ...