Solution -「HDU #6566」The Hanged Man
\(\mathcal{Description}\)
Link.
给定一棵含 \(n\) 个点的树,每个结点有两个权值 \(a\) 和 \(b\)。对于 \(k\in[1,m]\),分别求
\]
其中 \(S\) 是树上的一个独立点集。
测试数据组数 \(\le20\),\(n\le50\),\(m\le5\times10^3\)。
\(\mathcal{Solution}\)
一个 naive 的 DP 方法是,令 \(f(u,i,0/1)\) 表示 \(u\) 子树内,入 \(S\) 的点 \(a\) 值之和为 \(i\),且 \(u\) 入/未入集时最大的 \(b\) 值之和及其方案数。这种做法慢在需要 \(\mathcal O(m^2)\) 合并两棵子树的 DP 信息。我们尝试转化成一种仅需要单点更新的 DP。
考虑点分树,因为结点 \(u\) 在原树上的邻接点必然是点分树上 \(u\) 的祖先或子树内的孩子,所以可以在点分树上以 DFN 的顺序 DP,令状态 \(f(i,s,j)\) 表示考虑了 DFN 在 \([1,i]\) 之间的结点,设 \(u\) 是 DFN 为 \(i\) 的结点,那么 \(s\) 为 \(u\) 的祖先中被选结点的深度集合,\(j\) 为 \(a\) 值之和,状态表示最大值及方案数。
以上,每次暴力插入一个结点更新 DP 数组,就能做到 \(\mathcal O(Tnm2^{\log n})=\mathcal O(Tn^2m)\) 了。
(我写 T 了,也许是代码问题,抱歉 qwq。
\(\mathcal{Code}\)
会 T 掉的可怜代码。
/* Clearink */
#include <cstdio>
#define rep( i, l, r ) for ( int i = l, repEnd##i = r; i <= repEnd##i; ++i )
#define per( i, r, l ) for ( int i = r, repEnd##i = l; i >= repEnd##i; --i )
typedef long long LL;
inline void chkmax( int& a, const int b ) { a < b && ( a = b, 0 ); }
const int MAXN = 100, MAXM = 5e3;
int n, m, a[MAXN + 5], b[MAXN + 5];
int siz[MAXN + 5], wgt[MAXN + 5], ban[MAXN + 5], dep[MAXN + 5];
bool vis[MAXN + 5];
struct Graph {
int ecnt, head[MAXN + 5], to[MAXN * 2 + 5], nxt[MAXN * 2 + 5];
inline void operator () ( const int u, const int v ) {
to[++ecnt] = v, nxt[ecnt] = head[u], head[u] = ecnt;
}
inline void clear( const int n ) {
ecnt = 0;
rep ( i, 1, n ) head[i] = 0;
}
} srcT, divT;
#define adj( T, u, v ) \
for ( int i = T.head[u], v; v = T.to[i], i; i = T.nxt[i] )
struct Atom {
int val; LL cnt;
Atom(): val( 0 ), cnt( 0 ) {}
Atom( const int v, const LL c ): val( v ), cnt( c ) {}
inline Atom operator + ( const int x ) {
return cnt ? Atom( val + x, cnt ) : Atom();
}
inline Atom& operator += ( const Atom& u ) {
if ( val == u.val ) cnt += u.cnt;
if ( val < u.val ) *this = u;
return *this;
}
} f[MAXN * 4 + 5][MAXM + 5];
inline void findG( const int u, const int fa, const int all, int& rt ) {
siz[u] = 1, wgt[u] = 0;
adj ( srcT, u, v ) if ( !vis[v] && v != fa ) {
findG( v, u, all, rt ), siz[u] += siz[v];
chkmax( wgt[u], siz[v] );
}
chkmax( wgt[u], all - siz[u] );
if ( !rt || wgt[u] < wgt[rt] ) rt = u;
}
inline void build( const int u ) {
vis[u] = true;
adj ( srcT, u, v ) if ( !vis[v] ) {
int rt = 0; findG( v, 0, siz[v], rt );
divT( u, rt ), dep[rt] = dep[u] + 1;
build( rt );
}
}
inline void solve( const int u, int& lasd ) {
int d = dep[u];
rep ( i, 1 << d, ( 1 << lasd << 1 ) - 1 ) {
Atom *cf = f[i & ( ( 1 << d ) - 1 )], *lf = f[i];
rep ( j, 0, m ) cf[j] += lf[j], lf[j] = Atom();
}
rep ( i, 0, ( 1 << d ) - 1 ) {
if ( i & ban[u] ) continue;
Atom *cf = f[i | 1 << d], *lf = f[i];
rep ( j, 0, m - a[u] ) cf[j + a[u]] += lf[j] + b[u];
}
lasd = d;
adj ( divT, u, v ) solve( v, lasd );
}
int main() {
int T; scanf( "%d", &T );
rep ( cas, 1, T ) {
scanf( "%d %d", &n, &m );
srcT.clear( n ), divT.clear( n );
rep ( i, 1, n ) vis[i] = false;
// f[][] is cleared when outputing;
// ban[] is cleared when calculating.
rep ( i, 1, n ) scanf( "%d %d", &a[i], &b[i] );
for ( int i = 2, u, v; i <= n; ++i ) {
scanf( "%d %d", &u, &v );
srcT( u, v ), srcT( v, u );
}
int rt = 0; findG( 1, 0, n, rt );
build( rt );
rep ( u, 1, n ) {
ban[u] = 0;
adj ( srcT, u, v ) ban[u] |= ( dep[v] < dep[u] ) << dep[v];
}
f[0][0].cnt = 1;
int lasd = 0; solve( rt, lasd );
printf( "Case %d:\n", cas );
rep ( i, 1, m ) {
Atom ans;
rep ( j, 0, ( 1 << lasd << 1 ) - 1 ) {
ans += f[j][i], f[j][i] = Atom();
}
printf( "%lld%c", ans.cnt, i ^ m ? ' ' : '\n' );
}
}
return 0;
}
Solution -「HDU #6566」The Hanged Man的更多相关文章
- Solution -「HDU 6875」Yajilin
\(\mathcal{Description}\) Link.(HDU 裂开了先放个私链 awa.) 在一个 \(n\times n\) 的方格图中,格子 \((i,j)\) 有权值 \(w_ ...
- Solution -「HDU 5498」Tree
\(\mathcal{Description}\) link. 给定一个 \(n\) 个结点 \(m\) 条边的无向图,\(q\) 次操作每次随机选出一条边.问 \(q\) 条边去重后构成生成 ...
- Solution -「HDU 6643」Ridiculous Netizens
\(\mathcal{Description}\) Link. 给定一棵含有 \(n\) 个结点的树,点 \(u\) 有点权 \(w_u\),求树上非空连通块的数量,使得连通块内点权积 \(\ ...
- Solution -「HDU 1788」CRT again
\(\mathcal{Description}\) Link. 解同余方程组: \[x\equiv m_i-a\pmod{m_i} \] 其中 \(i=1,2,\dots,n\). \ ...
- Solution -「ARC 104E」Random LIS
\(\mathcal{Description}\) Link. 给定整数序列 \(\{a_n\}\),对于整数序列 \(\{b_n\}\),\(b_i\) 在 \([1,a_i]\) 中等概率 ...
- Solution -「HDU」Professor Ben
Description 有 \(Q\) 个询问.每次给定一个正整数 \(n\),求它的所有因数的质因数个数的和. Solution 就讲中间的一个 Trick. 我们定义正整数 \(x\) 有 \(f ...
- Solution -「CTS 2019」「洛谷 P5404」氪金手游
\(\mathcal{Description}\) Link. 有 \(n\) 张卡牌,第 \(i\) 张的权值 \(w_i\in\{1,2,3\}\),且取值为 \(k\) 的概率正比于 \ ...
- Solution -「BZOJ 3812」主旋律
\(\mathcal{Description}\) Link. 给定含 \(n\) 个点 \(m\) 条边的简单有向图 \(G=(V,E)\),求 \(H=(V,E'\subseteq E)\ ...
- Solution -「CF 1342E」Placing Rooks
\(\mathcal{Description}\) Link. 在一个 \(n\times n\) 的国际象棋棋盘上摆 \(n\) 个车,求满足: 所有格子都可以被攻击到. 恰好存在 \(k\ ...
随机推荐
- BIO、NIO、AIO --- 个人理解
1.前言 什么是 BIO.NIO.AIO ,不难看出,都是共同的字符IO , IO的意思是input output ,即输入输出 , 那么 B . N .A 分别指不同的io模型 ,而io又分为 ...
- 理解闭包--js面向对象编程
什么是闭包? 先看一段代码: function a(){ var n = 0; function inc() { n++; console.log(n); } inc(); inc(); } a(); ...
- python实现掘金定时签到抽奖
python实现掘金定时签到抽奖 一. 概述 这里记录一下使用 python 实现掘金定时签到抽奖.首先需要登录掘金,进入签到页面,按 F12 打开浏览器的调试面板,选择 Network,选择 XHR ...
- BUU XSS COURSE 1 & [CISCN2019 华东北赛区]Web2
BUU XSS COURSE 1 & [CISCN2019 华东北赛区]Web2 XSS的题目没怎么做过,比赛基本上也没有(=_=),总结下今天做的两题 都是看赵总视频现学现做的,这里附上学习 ...
- TextBox,RichTextBox设置行高
/// <summary> /// 设置行距 /// </summary> /// <param name="ctl">控件</param ...
- React教程
教程 一.demo <!DOCTYPE html> <html> <head> <meta charset="UTF-8"/> &l ...
- lvs的三种模式
一.NAT模式(VS-NAT) 原理:就是把客户端发来的数据包的IP头的目的地址,在负载均衡器上换成其中一台RS的IP地址,并发至此RS来处理,RS处理完成后把数据交给经过负载均衡器,负载均衡器再把数 ...
- 【webpack4.0】---webpack的基本使用(一)
一.初识webpack 1.什么是webpack? WebPack可以看做是模块打包机:它做的事情是,分析你的项目结构,找到JavaScript模块以及其它的一些浏览器不能直接运行的拓展语言(Scss ...
- 笔记:Bridging the Gap Between Relevance Matching and Semantic Matching for Short Text Similarity Modeling
笔记:Bridging the Gap Between Relevance Matching and Semantic Matching for Short Text Similarity Model ...
- React之redux学习日志(redux/react-redux/redux-saga)
redux官方中文文档:https://www.redux.org.cn/docs/introduction/CoreConcepts.html react-redux Dome:https://co ...