HDU-4219-Randomization?
题目描述
给定一棵\(n\)个节点的树,每条边的权值为\([0,L]\)之间的随机整数,求这棵树两点之间最长距离不超过\(S\)的概率。
Input
第一行三个整数\(n,L,S\)
接下来n-1行,每行两个整数,表示树上的一条边。
- \(1 <= T <= 512\)
- \(1 <= N <= 64\)
- \(1 <= L <= 8, 1 <= S <= 512\)
加强版:
- \(T=1\)
- \(1 <= N <= 100\)
- \(1 <= L <= 10, 1 <= S <= 2000\)
Output
输出树上任意两点距离都不超过\(S\)的概率,误差在\(10^{-6}\)以内都在算正确
Sample Input
3
2 3 2
1 2
4 3 4
1 2
2 3
3 4
7 4 10
1 2
2 3
4 5
2 6
4 7
4 6
Sample Output
Case 1: 0.750000
Case 2: 0.500000
Case 3: 0.624832
让我们先来考虑一下未加强的数据该怎么去写?
题目中树上两点最大距离就是树的直径。
因此,对于一个节点,我们需要维护信息是该节点到其子树的最大值的距离的概率。
同时,我们还要把儿子向父亲转移,即维护儿子到另一个儿子的距离的概率。
那么,我们就可以直接树形\(DP\)了。
\(dp[i][j]\)表示节点\(i\)到其子树的最大值为\(j\)的概率。
\(Dis[i]\)表示当前节点的子树到父亲的距离为\(i\)的概率。
\(temp[i]\)进行储存\(dp\)值,因为一个\(dp\)值会被多次利用。
总时间复杂度:\(O(n*S*(S+L))\)。
对于这个算法是无法过加强版的。
我们可以发现,这个\(dp\)的时间复杂度主要累积到\(temp\)值的累积上。
并且\(temp\)值的累积一定是连续的一段区间。
于是我们就可以利用前缀和优化,最后再把重复的给减掉。
时间复杂度:\(O(n*S*L)\)。
代码如下
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define u64 unsigned long long
#define reg register
#define Raed Read
#define debug(x) cerr<<#x<<" = "<<x<<endl;
#define rep(a,b,c) for(reg int a=(b),a##_end_=(c); a<=a##_end_; ++a)
#define ret(a,b,c) for(reg int a=(b),a##_end_=(c); a<a##_end_; ++a)
#define drep(a,b,c) for(reg int a=(b),a##_end_=(c); a>=a##_end_; --a)
#define erep(i,G,x) for(reg int i=(G).Head[x]; i; i=(G).Nxt[i])
inline int Read() {
int res = 0, f = 1;
char c;
while (c = getchar(), c < 48 || c > 57)if (c == '-')f = 0;
do res = (res << 3) + (res << 1) + (c ^ 48);
while (c = getchar(), c >= 48 && c <= 57);
return f ? res : -res;
}
template<class T>inline bool Min(T &a, T const&b) {
return a > b ? a = b, 1 : 0;
}
template<class T>inline bool Max(T &a, T const&b) {
return a < b ? a = b, 1 : 0;
}
const int N = 2e3+5, M = 1005, mod = 1e9 + 7;
bool MOP1;
int n,l,s,Case;
struct Link_list {
int Tot,Head[N],Nxt[M<<1],to[M<<1],cost[M<<1];
inline void Init(void) {
Tot=0;
memset(Head,0,sizeof Head);
}
inline void AddEdgepair(int a,int b) {
to[++Tot]=b,Nxt[Tot]=Head[a],Head[a]=Tot;
to[++Tot]=a,Nxt[Tot]=Head[b],Head[b]=Tot;
}
} G;
struct T290 {
double Dis[N],dp[N][N],temp[N];
void dfs(int x,int pre) {
rep(i,1,s)dp[x][i]=0;
dp[x][0]=1;
rep(j,0,s)Dis[j]=temp[j]=0;
erep(i,G,x) {
int y=G.to[i];
if(y==pre)continue;
dfs(y,x);
rep(j,0,s)Dis[j]=temp[j]=0;
rep(j,0,l)rep(k,0,s)if(j+k<=s)Dis[j+k]+=dp[y][k]*1.0/(l+1);
else break;
rep(j,0,s)rep(k,0,s-j)temp[max((int)k,(int)j)]+=dp[x][k]*Dis[j];
rep(j,0,s)dp[x][j]=temp[j];
}
}
inline void solve(void) {
ret(i,1,n) {
int a=Raed(),b=Read();
G.AddEdgepair(a,b);
}
dfs(1,0);
double Ans=0;
rep(i,0,s)Ans+=dp[1][i];
printf("Case %lld: %.6lf\n",++Case,Ans);
}
} Pbl;
struct T2100 {
double Dis[N],dp[N][N],temp[N],Sum1[N],Sum2[N];
void dfs(int x,int pre) {
rep(i,1,s)dp[x][i]=0;
dp[x][0]=1;
erep(i,G,x) {
int y=G.to[i];
if(y==pre)continue;
dfs(y,x);
rep(j,0,s)Dis[j]=temp[j]=0;
rep(j,0,l)rep(k,0,s)if(j+k<=s)Dis[j+k]+=dp[y][k]*1.0/(l+1);
else break;
Sum1[0]=dp[x][0];
Sum2[0]=Dis[0];
rep(i,1,s) {
Sum1[i]=Sum1[i-1]+dp[x][i];
Sum2[i]=Sum2[i-1]+Dis[i];
}
rep(j,0,s) {
int k=min(j,s-j);
temp[j]+=Sum1[k]*Dis[j];
}
rep(j,0,s) {
int k=min(j,s-j);
temp[j]+=Sum2[k]*dp[x][j];
}
rep(j,0,s)if(j*2<=s)temp[j]-=dp[x][j]*Dis[j];
rep(j,0,s)dp[x][j]=temp[j];
}
}
inline void solve(void) {
ret(i,1,n) {
int a=Raed(),b=Read();
G.AddEdgepair(a,b);
}
dfs(1,0);
double Ans=0;
rep(i,0,s)Ans+=dp[1][i];
printf("Case %lld: %.6lf\n",++Case,Ans);
}
} P100;
bool MOP2;
inline void _main(void) {
int T=Read();
while(T--) {
G.Init();
n=Raed(),l=Read(),s=Read();
P100.solve();
}
}
signed main() {
#define offline1
#ifdef offline
freopen("random.in", "r", stdin);
freopen("random.out", "w", stdout);
_main();
fclose(stdin);
fclose(stdout);
#else
_main();
#endif
return 0;
}
HDU-4219-Randomization?的更多相关文章
- HDOJ 2111. Saving HDU 贪心 结构体排序
Saving HDU Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total ...
- 【HDU 3037】Saving Beans Lucas定理模板
http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...
- hdu 4859 海岸线 Bestcoder Round 1
http://acm.hdu.edu.cn/showproblem.php?pid=4859 题目大意: 在一个矩形周围都是海,这个矩形中有陆地,深海和浅海.浅海是可以填成陆地的. 求最多有多少条方格 ...
- HDU 4569 Special equations(取模)
Special equations Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u S ...
- HDU 4006The kth great number(K大数 +小顶堆)
The kth great number Time Limit:1000MS Memory Limit:65768KB 64bit IO Format:%I64d & %I64 ...
- HDU 1796How many integers can you find(容斥原理)
How many integers can you find Time Limit:5000MS Memory Limit:32768KB 64bit IO Format:%I64d ...
- hdu 4481 Time travel(高斯求期望)(转)
(转)http://blog.csdn.net/u013081425/article/details/39240021 http://acm.hdu.edu.cn/showproblem.php?pi ...
- HDU 3791二叉搜索树解题(解题报告)
1.题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=3791 2.参考解题 http://blog.csdn.net/u013447865/articl ...
- hdu 4329
problem:http://acm.hdu.edu.cn/showproblem.php?pid=4329 题意:模拟 a. p(r)= R'/i rel(r)=(1||0) R ...
- HDU 2586
http://acm.hdu.edu.cn/showproblem.php?pid=2586 题意:求最近祖先节点的权值和 思路:LCA Tarjan算法 #include <stdio.h&g ...
随机推荐
- js文件夹上传
文件夹上传:从前端到后端 文件上传是 Web 开发肯定会碰到的问题,而文件夹上传则更加难缠.网上关于文件夹上传的资料多集中在前端,缺少对于后端的关注,然后讲某个后端框架文件上传的文章又不会涉及文件夹. ...
- 洛谷 P2473 [SCOI2008]奖励关 ( 期望DP )
题目链接 题意 : 中文题.点链接 分析 : 第一道有关概率期望的DP 有个大部分情况下通用的结论 概率正推.期望反推 原因不明.其实是没有查到较好的解释 这题由于有一些取物品的先决条件在这里 而且观 ...
- TTTTTTTTTTTT CF 653D 送邮递员
链接:给一张n个点m条带权边的有向图,有x个人从起点出发到终点,每个人带的都带相同重量的货物, 规定一条边最多能经过其上权的重量的货物,问最多能带多重的货物? 2 ≤ n ≤ 50, 1 ≤ m ≤ ...
- Pod初始化容器之Init Container
Init 容器的介绍 Pod能够具有多个容器,应用运行在容器里面,但是它也可能有一个或多个先于应用容器启动的 Init容器Init 容器与普通的容器非常像,除了如下两点: c Init 容器总是运行 ...
- 12.Python数值类型(整形、浮点型和复数)及其用法
实际开发中,我们经常需要使用数字记录游戏中用户的得分.游戏中角色的生命值.伤害值等信息,Python 语言提供了数值类型用于保存这些数值. 需要注意的是,Python 中这些数值类型都是不可改变的,也 ...
- SVN提交大量无效文件补救方法
有的时候,使用SVN时候会发现,由于系统编译器的问题,会自动生成大量.class文件, 或者一些多余的配置文件,这里主要就是整理一下,当如果手误,将这些多余文件都提交到了svn上面的补救方法. 可以在 ...
- 使用宝塔面板 配置nginx 访问ftp服务器下面的图片
如果 你在服务器上 运行war项目 可以在tomcat 配置访问的: tomcat 也贴出来吧! 一.tomca配置访问,需要更改配置文件server.xml ,如果找不到,自己好好找一下 一般在 ...
- 记一次elastic-job使用
当当的elastic-job定时任务 业务场景是定时从微信取accesstoken和jsticket,因为都只有7200秒的有效时间,所以设置了定时任务,定时将得到的数据存到redis缓存中 问题1: ...
- 20175215 2018-2019-2 第二周java课程学习总结
一.学生免费申请使用IDEA 下载好IDEA后,设置到最后有一个界面, 我们需要到IDEA官网进行IDEA免费试用权的申请,如果有学校的邮箱,使用学校的邮箱注册并证明是自己的就可以直接通过申请.如下图 ...
- ASP.NET postback with JavaScript (UseSubmitBehavior)
ASP.NET postback with JavaScript Here is a complete solution Entire form tag of the asp.net page < ...