$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 城市地铁规划的更多相关文章

  1. 【LibreOJ】#6395. 「THUPC2018」城市地铁规划 / City 背包DP+Prufer序

    [题目]#6395. 「THUPC2018」城市地铁规划 / City [题意]给定n个点要求构造一棵树,每个点的价值是一个关于点度的k次多项式,系数均为给定的\(a_0,...a_k\),求最大价值 ...

  2. 大数据freestyle: 共享单车轨迹数据助力城市合理规划自行车道

    编者按:近年来,异军突起的共享单车极大地解决了人们共同面临的“最后一公里”难题,然而,共享单车发展迅猛,自行车道建设却始终没有能够跟上脚步.幸运的是摩拜单车大量的轨迹数据为我们提供了一种新的思路:利用 ...

  3. ZROI 暑期高端峰会 A班 Day6 DP

    [THUPC2018]城市地铁规划 (日常讲题之前 YY--) 一眼出 \(O(n^3+nk)\) 做法. \(dp[i][j]\) 表示前 \(i\) 个点,前 \(i\) 个点度数和为 \(j\) ...

  4. zhengrui集训笔记2

    Day_6 计算几何 点积\Large 点积点积 叉积\Large 叉积叉积 极角\Large 极角极角 < π\piπ :叉积判断 else :atan2 旋转\Large 旋转旋转 左乘第一 ...

  5. 中科燕园GIS外包---地铁GIS项目

    (1)地铁保护及project地质管理     • 地铁保护     地铁交通既有运量大,速度快的特点,又有差别于其它交通方式的在地下执行的空间特殊性,因此地铁的保护显得尤为重要. 首先必须编制完整的 ...

  6. SVG的动态之美-搜狗地铁图重构散记

    搜狗地图发布了新版的移动端地铁图,改版初衷是为了用户交互体验的提升以及性能的改善.原版地铁图被用户吐槽最多的是pinch缩放不流畅.无过渡动画.拖拽边界不合理等等,大体上都是交互体验上的问题.实际上原 ...

  7. 地铁合作(hjy)

       用时大概时间 18 个小时  我的队友是 :黄珺瑜我们一开始的想法是自己建立数据库,用来查询路线,再调用一下百度api作为地图(仅作地图没有任何操作的那种):后来我们的步骤也的确是这样,我们先确 ...

  8. 机器人与机器人仿真技术(zz)

    http://www.viblue.com/archives/5587.htm 一.机器人简介: 机器人(Robot)是自动执行工作的机器装置.它既可以接受人类指挥,又可以运行预先编排的程序,也可以根 ...

  9. 2018年度最优秀mac软件及游戏推荐,个个万里挑一

    今天和大家带来2018年度最优秀Mac软件和游戏合集,个个万里挑一,2018年,风云社区(scoee.com)分享了上数千款优秀的Mac软件和游戏,结合用户反馈,精选出各个类别的优秀的Mac软件,推荐 ...

随机推荐

  1. sql server分区

    1. 创建分区 分区步骤:1.创建分区函数  2.创建分区架构 3.创建分区索引(聚集)   --1. 创建分区函数 DECLARE @dt datetime SET @dt = '20030901' ...

  2. loadrunder脚本篇——执行操作系统命令

    思路: 用loadrunner system()函数 函数原型: int system( const char *string ); 示例一:在指定目录下创建指定文件 Action() { char ...

  3. 【转】dmesg 时间转换

    dmesg 时间转换 dmesg 输出的格式不易查看,可以通过命令进行转换. 记录如下: 时间查看: date -d "1970-01-01 UTC `echo "$(date + ...

  4. deeplink

    http://www.cnblogs.com/shadajin/p/5724117.html Deeplink,简单讲,就是你在手机上点击一个链接之后,可以直接链接到app内部的某个页面,而不是app ...

  5. EGLImage与纹理

    http://blog.csdn.net/sunnytina/article/details/51895406 Android使用Direct Textures提高glReadPixels.glTex ...

  6. 0530JavaScript基础2

    常用内置对象 所谓内置对象就是ECMAScript提供出来的一些对象,我们知道对象都是有相应的属性和方法 数组Array(部分相当于列表) 1.数组的创建方式 var colors = ['red', ...

  7. Vue.js学习笔记 第五篇 事件处理

    监听事件 <!DOCTYPE html> <html> <head> <meta charset="utf-8"> <titl ...

  8. HRBUST 1717 字典树模板

    之前写字典树虽然很熟也能变化 但是一直是到了场上再乱写 写的很长 于是准备写个短点的板子 于是选了个水题 然而写出了1W个bug insert和query反而写的没什么问题.. L c[100050] ...

  9. mssql 函数STUFF 的用法

    DECLARE @limitDay INT;SET @limitDay = 92;IF DATEDIFF(DAY, '2017-12-13 00:00:00', '2017-12-13 18:00:0 ...

  10. JavaWeb -- 文件上传下载示例

    1. 上传简单示例 Jsp <%@ page language="java" import="java.util.*" pageEncoding=&quo ...