传送门

题目大意

给定一张$n$个点的有向图$(n\leq 1000)$,每次操作定义为等概率随机选定一个未被标记的点,从它出发将所有它和它能到达的点进行标记,询问期望进行多少次操作。

题解

考虑每个点被选定的概率。

它被直接选定当且仅当它在能到达它的点被选定之前选定。

对于一个点$x$设能到达它的点(包括$x$本身)有$G_x$个,那么$x$被直接选定的概率为$\frac{1}{G_x}$,可以理解为它对答案的贡献期望是$\frac{1}{G_x}$。

处理$G_x$即可。

可以先使用$tarjan$将原图缩成有向无环图,进行拓扑排序时维护能到达一个点的点集,利用$bitset$优化传递。

设一个强联通分量$x$的点数为$Size_x$,能到达$x$的点数之和(包括$x$)为$G_x$。

则答案$=\sum \frac{Size_x}{G_x}$。

复杂度大概是$O(T\frac{n^3}{32})$

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<bitset>
#include<cmath>
#define LL long long
#define M 1020
using namespace std;
int read(){int nm=0,fh=1; char cw=getchar();for(;!isdigit(cw);cw=getchar()) if(cw=='-') fh=-fh;for(;isdigit(cw);cw=getchar()) nm=nm*10+(cw-'0');return nm*fh;
}
int n,m,c[M][M],sz[M],vis[M],dfn[M],low[M],cnt,S[M],top,tk[M],hd,tl,q[M];
int fs[M],to[M*M],tmp,nt[M*M],be[M],tot,ct[M],ind[M]; bool gt[M][M]; bitset<M>G[M];
void link(int x,int y){nt[tmp]=fs[x],fs[x]=tmp,to[tmp++]=y,ind[y]++;}
void dfs(int x){cnt++,dfn[x]=low[x]=cnt,S[++top]=x,S[top+1]=0;for(int i=1;i<=ct[x];i++){if(!dfn[c[x][i]]) dfs(c[x][i]);if(!be[c[x][i]]) low[x]=min(low[x],low[c[x][i]]);}if(low[x]<dfn[x]) return; sz[++tot]=0,G[tot].reset(),ind[tot]=0;while(S[top+1]!=x) be[S[top]]=tot,sz[tot]++,G[tot].set(S[top]),top--;
}
int main(){for(int tt=1,T=read();tt<=T;tt++){tot=cnt=hd=tl=tmp=0,n=read(),memset(vis,0,sizeof(vis)),memset(dfn,0,sizeof(dfn));memset(fs,-1,sizeof(fs)),memset(gt,false,sizeof(gt)),memset(be,0,sizeof(be));for(int i=1;i<=n;i++) for(ct[i]=0,m=read();m;--m) c[i][++ct[i]]=read();for(int i=1;i<=n;i++) if(!dfn[i]) dfs(i);for(int i=1;i<=n;i++){for(int j=1;j<=ct[i];j++)if(be[i]!=be[c[i][j]]) gt[be[i]][be[c[i][j]]]=true;}for(int i=1,res=0;i<=tot;res=0,i++){for(int j=1;j<=tot;j++)  if(i!=j&&gt[j][i]) res=1,link(j,i);if(!res) q[tl++]=i;}while(hd<tl){int x=q[hd++];for(int i=fs[x];i!=-1;i=nt[i]){G[to[i]]|=G[x]; if(!(--ind[to[i]])) q[tl++]=to[i];}}double ans=0.0;for(int i=1;i<=tot;i++) ans+=(sz[i]*1.0)/(G[i].count()*1.0);printf("Case #%d: %.5f\n",tt,ans);}return 0;
}

  

转载于:https://www.cnblogs.com/OYJason/p/9745058.html

Hdu 5036 Explosion相关推荐

  1. hdu 5036 Explosion bitset优化floyd

    http://acm.hdu.edu.cn/showproblem.php?pid=5036 题意就是给定一副有向图,现在需要走遍这n个顶点,一开始出发的顶点是这n个之中的随便一个. 如果走了1,那么 ...

  2. 2014 ACM/ICPC 北京赛区网络赛解题报告汇总

    首页 算法竞赛» 信息聚合 ONLINE JUDGE 书刊杂志 BLOG» 新闻故事» 招聘信息» 投稿须知 2014 ACM/ICPC 北京赛区网络赛解题报告汇总 九月 21st, 2014 | P ...

  3. HDU 4389 - X mod f(x)

    题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=4389 2012多校,第9场,1010 . 问题是,询问区间内 存在多少个 哈沙德数(Harshad ...

  4. hdu 4389 囧,打表

    http://acm.hdu.edu.cn/showproblem.php?pid=4389 题意 :一个数能被他各个位数之和整除则符合要求,给L,R,问区间里有多少个数符合要求. 囧,居然打表就能过 ...

  5. HDU——1106排序(istringstream的使用、STLvector练习)

    排序 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submiss ...

  6. hdu 5438 Ponds 拓扑排序

    Ponds Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/contests/contest_showproblem ...

  7. HDU 1248 寒冰王座(全然背包:入门题)

    HDU 1248 寒冰王座(全然背包:入门题) http://acm.hdu.edu.cn/showproblem.php?pid=1248 题意: 不死族的巫妖王发工资拉,死亡骑士拿到一张N元的钞票 ...

  8. hdu 1312 Red and Black 解题报告

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1312 第二条深搜,题目并不难,但是做了我好久好久,由于一个细节,让我赌上了一个晚上的时间. 题目大意: ...

  9. HDU 1429 胜利大逃亡(续) (BFS+位压缩)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1429 胜利大逃亡(续) Time Limit: 4000/2000 MS (Java/Others)  ...

最新文章

  1. LINQ圣经——《LINQ实战》
  2. BGP 最佳路径选择之 -- Origin
  3. MySQL分区的限制(最多有多少个分区)
  4. 使用PostSharp开始AOP
  5. 存储Image到XML、显示Image
  6. CTFshow 反序列化 web255
  7. 编程之美-队列中取最大值操作问题
  8. To B设计系统 - 在平平淡淡中开花结果
  9. wireshark抓包过滤指定的字节数据
  10. .NET Core HttpClient请求异常分析
  11. 使用setCustomView建立标题居中的自定义ActionBar
  12. javatodo框架中怎么配置路由
  13. JS获取IP、MAC和主机名的五种方法
  14. 中压变频器行业现状调研及趋势分析报告
  15. 计算机操作系统(第四版)第二章 习题答案
  16. DirectShow之视频渲染
  17. GitHut上最受关注的Objective-C项目
  18. 百位LOL英雄联盟角色合集
  19. 用户反馈驱动抖音产品体验优化实践
  20. Outline浏览加速-实操分享今日随笔

热门文章

  1. 如何查看服务器gpu性能,ubuntu服务器查看GPU和CPU实时使用情况
  2. ue查找文件中每行第二个单词_UI设计和UE/UX设计有什么区别?它们的晋升路径是什么?...
  3. 信息安全工程师笔记-恶意代码防范技术原理
  4. Spring Boot笔记-get请求发送json数据(方便前端vue解析)
  5. HTTP中CORS跨域请求的实现(C++|Qt框架实现)
  6. Qt + Python + OpenCV图标替换工具 之 项目介绍(一)
  7. C++|STL学习笔记-对STL中关联容器map的进一步认识
  8. html美甲预约网页制作,教程.html · NFS1077/FNM美甲店 - Gitee.com
  9. php调用视频功能,phpcms如何调用视频?
  10. javascript 动态修改css样式方法汇总(四种方法)