[HEOI2013] SAO(dp + 组合数 + 前缀和)

problem

luogu-P4099

solution

两篇前提题解:排列计数,老C的键盘

想必已经看了 CQOI2017 老C的键盘 一题题解了。

这里直接考虑优化状态转移方程。

我们发现 f(u,i),f(v,j)f(u,i),f(v,j)f(u,i),f(v,j),当枚举 jjj 后,对应的 kkk 是一段连续区间。

f(u,i)∗f(v,j)∗(k−1i−1)∗(sizu+sizv−ksizu−i)f(u,i)*f(v,j)*\binom{k-1}{i-1}*\binom{siz_u+siz_v-k}{siz_u-i}f(u,i)f(v,j)(i1k1)(sizuisizu+sizvk) 转移贡献变化的只有 kkk,且有两项都与之挂钩。

我们不妨反过来固定 kkk,能转移到这个 kkk 的肯定也是一段连续的 jjj

所以我们求出来 f(v,j)f(v,j)f(v,j) 后,直接做个前缀和,f(v,j)=∑j′≤jf(v,j′)f(v,j)=\sum_{j'\le j}f(v,j')f(v,j)=jjf(v,j)

以要求 hu<hvh_u<h_vhu<hv 为例。

其原来是 i+j≤k≤i+sizvi+j\le k\le i+siz_vi+jki+sizv,固定 i,ki,ki,k 后,能转移过来的 jjj 范围则为 1∼k−i1\sim k-i1ki

另一种则是 k−i<j≤sizvk-i<j\le siz_vki<jsizv

两种的 kkk 范围也稍稍不同。

具体可以看代码实现。

code

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define mod 1000000007
#define maxn 1005
char s[5];
int n;
int c[maxn][maxn], f[maxn][maxn], g[maxn], siz[maxn];
vector < pair < int, int > > G[maxn];void dfs( int u, int fa ) {f[u][1] = siz[u] = 1;for( int i = 0;i < G[u].size();i ++ ) {int v = G[u][i].first, w = G[u][i].second;if( v == fa ) continue;dfs( v, u );memcpy( g, f[u], sizeof( f[u] ) );memset( f[u], 0, sizeof( f[u] ) );for( int i = 1;i <= siz[u];i ++ )if( w ) {for( int k = i + 1;k <= i + siz[v];k ++ )(f[u][k] += f[v][k - i] % mod * g[i] % mod * c[k - 1][i - 1] % mod * c[siz[u] + siz[v] - k][siz[u] - i]) %= mod;}else {for( int k = i;k < i + siz[v];k ++ )(f[u][k] += (f[v][siz[v]] - f[v][k - i]) % mod * g[i] % mod * c[k - 1][i - 1] % mod * c[siz[u] + siz[v] - k][siz[u] - i]) %= mod;}siz[u] += siz[v];}for( int i = 1;i <= siz[u];i ++ ) (f[u][i] += f[u][i - 1]) %= mod;
}signed main() {for( int i = 0;i <= 1000;i ++ ) {c[i][0] = c[i][i] = 1;for( int j = 1;j < i;j ++ )c[i][j] = (c[i - 1][j - 1] + c[i - 1][j]) % mod;}int T;scanf( "%lld", &T );while( T -- ) {scanf( "%lld", &n );memset( f, 0, sizeof( f ) );for( int i = 0;i < n;i ++ ) G[i].clear();for( int i = 1, u, v;i < n;i ++ ) {scanf( "%lld %s %lld", &u, s, &v );G[u].push_back( make_pair( v, s[0] == '<' ) );G[v].push_back( make_pair( u, s[0] == '>' ) );}dfs( 0, 0 );printf( "%lld\n", ( f[0][n] + mod ) % mod );}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316350.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

cf1555A. PizzaForces

cf1555A. PizzaForces A. PizzaForces 题意&#xff1a; 有三种披萨&#xff0c;第一种有六块&#xff0c;需要花费15分钟&#xff0c;第二种有8块&#xff0c;需要花费20分钟&#xff0c;第三问有10块&#xff0c;需要花费25分钟。 现在要吃x块披萨&#xff0c;问最少时间花…

DI是实现面向切面和面向抽象的前提

DI越来越重要DI就是依赖注入&#xff0c;现在来说&#xff0c;大部分框架都是以DI为基础组件的&#xff0c;每一个框架都有自己的DI组件&#xff0c;像dotnet core&#xff0c;java spring等&#xff0c;也都为自己的框架量身打造了DI工具。面向对象的几个原则依赖倒置原则&…

[CQOI2017] 老C的任务(差分 + 树状数组 / K-D tree)

problem luogu-P3755 solution 这题第一眼矩阵内的点权值和&#xff0c;马上就是 K-D tree\text{K-D tree}K-D tree 不过脑子的敲。 这其实就是个二维数点问题&#xff0c;完全可以树状数组。 将矩阵差分成四个以原点为左下角的矩阵。 然后将基站按 xxx 轴排序&#xff0…

.net core 并发下的线程安全问题

抱歉&#xff0c;其实内容并不如题&#xff01;&#xff01;&#xff01;背景&#xff08;写测试demo所出现的异常&#xff0c;供大家学习与拍砖&#xff09;&#xff1a;.net core webapi项目&#xff0c;做了一个授权的filter&#xff08;真正的生产项目的话&#xff0c;JWT很…

cf1555B. Two Tables

cf1555B. Two Tables 题意&#xff1a; 一个大矩阵空间内放置一个矩阵a&#xff0c;现在要再往这个空间内放一个矩阵b&#xff0c;a移动距离len才能放下b&#xff0c;问len最小是多少 题解&#xff1a; 不难发现左右或上下移动是最佳的&#xff0c;斜着移动是最不好的。此时…

cf1555C Coin Rows

cf1555C Coin Rows 题意&#xff1a; 有一个两行m列的地图&#xff0c;每个格子都有对应的价值&#xff0c;有a&#xff0c;b两个人&#xff0c;都从左上角到右下角&#xff0c;且都只能向右向下走&#xff0c;a先出发&#xff0c;a每到一个格子&#xff0c;就会获得这个地方…

C#并行编程(2):.NET线程池

线程 Thread在总结线程池之前&#xff0c;先来看一下.NET线程。.NET线程与操作系统(Windows)线程有什么区别&#xff1f;.NET利用Windows的线程处理功能。在C#程序编写中&#xff0c;我们首先会新建一个线程对象System.Threading.Thread&#xff0c;并为其指定一个回调方法&…

[CQOI2015] 任务查询系统(主席树)

problem luogu-P3168 solution 主席树板题。 将一个任务拆成 lil_ili​ 秒开始&#xff0c;ri1r_i1ri​1 秒结束的两个任务。 但不建议以每一秒作为一个主席树版本&#xff0c;因为一秒中可能有若干个任务开始或结束。 不妨将所有任务按时刻排序&#xff0c;然后以每个任…

ASP.NET Core launchsettings.json文件(8)《从零开始学ASP.NET CORE MVC》:

本文出自《从零开始学ASP.NET CORE MVC》推荐文章&#xff1a;ASP.NET Core 进程外(out-of-process)托管ASP.NET Core launchsettings.json文件在本视频中&#xff0c;我们将讨论在ASP.NET Core项目中launchsettings.json文件的重要性。launchsettings.json文件您将在项目根文件…

[CQOI2017] 老C的方块(网络流染色建图)

problem luogu-P3756 solution 据说要做网络流 24\text{24}24 题中的《方格取数问题》和《骑士共存问题》。 &#xff1f;&#xff1f;&#xff1f;那个不是直接最小割吗&#xff1f;哦原来是从黑白染色来理解的。我还是太水了。 这种题之所以能用网络流做&#xff0c;是因…

cf1555D. Say No to Palindromes

cf1555D. Say No to Palindromes 题意&#xff1a; 给出一个字符串&#xff0c;长度为n&#xff0c;而且都是a,b,c三个字符构成的&#xff0c;然后有m个询问 每个询问给出l r&#xff0c;问要想这个区间内任意长度字串都不是回文子串&#xff0c;至少要改多少个字符 题解&am…

江湖召集:.NET开发者们看过来,这场长沙的开发者技术大会正是为你精心准备的大餐...

看过去&#xff0c;历史的尘埃与沧海桑田古语有云“近代中国&#xff0c;湖南独撑半边天”&#xff0c;湖南长沙&#xff0c;作为湖南省的省会&#xff0c;自古以来便是各界风云人士兴起之地。随着互联网时代的到来&#xff0c;长沙&#xff0c;这座历史悠久的文化名城&#xf…

[HEOI2012] 朋友圈(最大团 + 结论 + 二分图 + 网络流)

problem luogu-P2423 solution 本题即求无向图最大团问题。这是个 NP hard\text{NP hard}NP hard 问题&#xff0c;所以必须从图的特殊性质出发&#xff0c;否则只能暴搜。 异或运算等价于二进制下不进位的加法运算。 observation1:\text{observation1}:observation1: AAA …

cf1555 E. Boring Segments

cf1555 E. Boring Segments 题意&#xff1a; 给你n个线段&#xff0c;最大点是m&#xff0c;每一个线段有一个权值w&#xff0c;你能选择线段来覆盖1-m这个区间的&#xff0c;选择的代价为最大权值和最小权值的差。问你最小的的代价是多少。 题解&#xff1a; 尺取线段树 …

C#并行编程(3):并行循环

初识并行循环并行循环主要用来处理数据并行的&#xff0c;如&#xff0c;同时对数组或列表中的多个数据执行相同的操作。在C#编程中&#xff0c;我们使用并行类System.Threading.Tasks.Parallel提供的静态方法Parallel.For和Parallel.ForEach来实现并行循环。从方法名可以看出&…

[NOI Online 2022 提高组] 丹钓战(单调栈 + 树状数组 / 主席树)

problem luogu-P8251 solution 按照题意模拟单调栈。 求出对于 iii 而言&#xff0c;当时单调栈的栈顶元素记为 pip_ipi​。 如果到 iii 时&#xff0c;栈顶已经为 pip_ipi​ 了&#xff0c;意味着这中间的所有元素要么是被 iii 弹出&#xff0c;要么就是被 iii 前面的某些…

Acwing 252. 树

Acwing 252. 树 题意&#xff1a; 给定一个有 N 个点&#xff08;编号 0,1,…,N−1&#xff09;的树&#xff0c;每条边都有一个权值&#xff08;不超过 1000&#xff09;。 树上两个节点 x 与 y 之间的路径长度就是路径上各条边的权值之和。 求长度不超过 K 的路径有多少条…

.net 4.5部署到docker容器

.NET FX 应用程序也是可以容器化的&#xff0c;容器化的选项有两个&#xff1a;部署到windows容器部署到linux容器部署到windows容器由于.net本身就是运行在windows平台的&#xff0c;所以它与windows容器也是更加适合&#xff0c;你可以以iis镜像为基础&#xff0c;去编写你的…

[NOI Online 2022 提高组] 讨论(巧妙的切入方式)

problem luogu-P8252 solution 本题最难处理的就是两个人会做的题目集合是包含关系的限制。 将所有人按会做的题数从大到小排序。 然后枚举 iii&#xff0c;只要这个人和之前某个人存在有至少一道公共的题目&#xff0c;并且保证这个人有新的题目&#xff0c;那么这两个人…

P4149 [IOI2011]Race

P4149 [IOI2011]Race 题意&#xff1a; 给一棵树&#xff0c;每条边有权。求一条简单路径&#xff0c;权值和等于 k&#xff0c;且边的数量最小。 题解&#xff1a; 用t[i]:长度为i的路径包含的最少边数 按照子树顺序&#xff0c;依次用dep[u]t[K-d[u]]更新ans&#xff0c;…