! 卡时过了 为什么我的这么慢?姿势不对??? -->谢 ws_fqk 我的Do(num[i])应该用 DO(root) 找了半天的root居然没有用....
define的教训永远忘不了了!!!
优化了好几处才过的...者...

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#define MAXN 30010
using namespace std;
int n,m,k;
bool b[MAXN];
int team[10000000],head,tail,d[MAXN];
int tot,g[MAXN],nnext[MAXN*2],num[MAXN*2],cost[MAXN*2];
int tot0,g0[MAXN],nnext0[MAXN*4],num0[MAXN*4],cost0[MAXN*4];
int tot1,g1[MAXN],nnext1[MAXN*4],num1[MAXN*4],cost1[MAXN*4];
struct H{int x,y,z;}G[MAXN*4];
bool cmp(H a,H b){if(a.x==b.x)return a.y>b.y;return a.x<b.x;}
void SPFA()
{memset(d,61,sizeof(d));d[1]=0;team[++tail]=1;while(head!=tail){int x=team[++head];b[x]=false;for(int i=g0[x];i;i=nnext0[i]){int tmp=num0[i];if(d[x]+cost0[i]<d[tmp]){d[tmp]=d[x]+cost0[i];if(!b[tmp]){b[tmp]=true;team[++tail]=tmp;}}}}
}
void Dfs_Road(int x)
{b[x]=true;for(int i=g1[x];i;i=nnext1[i]){int tmp=num1[i],z=cost1[i];if(!b[tmp]){tot++,nnext[tot]=g[x],g[x]=tot,num[tot]=tmp,cost[tot]=z;tot++,nnext[tot]=g[tmp],g[tmp]=tot,num[tot]=x;cost[tot]=z;Dfs_Road(tmp);}}
}
//************点分治***************
int root,sum,max_son[MAXN],size[MAXN],ans0,ans1,N[30000+1][2],T[30000+1][2];void Get_Root(int x,int fa)
{size[x]=1;max_son[x]=0;for(int i=g[x];i;i=nnext[i])if(num[i]!=fa&&!b[num[i]]){Get_Root(num[i],x);max_son[x]=max(max_son[x],size[num[i]]);size[x]+=size[num[i]];}max_son[x]=max(max_son[x],sum-size[x]);if(max_son[x]<max_son[root]) root=x;
}
void Dfs_T(int x,int fa,int depth,int c)
{if(depth>k) return;if(c>T[depth][0]) T[depth][0]=c,T[depth][1]=1;else if(c==T[depth][0]) T[depth][1]++;for(int i=g[x];i;i=nnext[i])if(!b[num[i]]&&num[i]!=fa)Dfs_T(num[i],x,depth+1,c+cost[i]);
}
void Do(int x)
{for(int i=0;i<=k;i++) N[i][0]=N[i][1]=0;N[0][1]=1;b[x]=true;for(int i=g[x];i;i=nnext[i])if(!b[num[i]]){for(int i=0;i<=k;i++) T[i][0]=T[i][1]=0;Dfs_T(num[i],-1,1,cost[i]);for(int i=1;i<=k;i++){if(ans0<N[k-i][0]+T[i][0])ans0=N[k-i][0]+T[i][0],ans1=N[k-i][1]*T[i][1];else if(ans0==N[k-i][0]+T[i][0])ans1+=N[k-i][1]*T[i][1];}for(int i=1;i<=k;i++){if(N[i][0]<T[i][0])N[i][0]=T[i][0],N[i][1]=T[i][1];else if(N[i][0]==T[i][0])N[i][1]+=T[i][1];}}for(int i=g[x];i;i=nnext[i])if(!b[num[i]]){root=0;sum=size[num[i]];Get_Root(num[i],-1);Do(root);}
}int main()
{cin>>n>>m>>k; k--;for(int i=1;i<=m;i++){int x,y,z;scanf("%d %d %d",&x,&y,&z);G[i]=(H){x,y,z};G[i+m]=(H){y,x,z};tot0++;nnext0[tot0]=g0[x],g0[x]=tot0,num0[tot0]=y;cost0[tot0]=z;tot0++;nnext0[tot0]=g0[y],g0[y]=tot0,num0[tot0]=x,cost0[tot0]=z;}SPFA(); sort(G+1,G+m*2+1,cmp);for(int i=1;i<=m*2;i++){int x=G[i].x,y=G[i].y,z=G[i].z;if(d[x]+z==d[y])tot1++,nnext1[tot1]=g1[x],g1[x]=tot1,num1[tot1]=y,cost1[tot1]=z;}memset(b,false,sizeof(b));Dfs_Road(1);memset(b,false,sizeof(b));root=0,sum=n,max_son[0]=n+1;Get_Root(1,-1);Do(root);cout<<ans0<<' '<<ans1;return 0;
}

转载于:https://www.cnblogs.com/ofsxb/p/5141635.html

【BZOJ 4016】[FJOI2014]最短路径树问题相关推荐

  1. bzoj 4016: [FJOI2014]最短路径树问题

    Description 给一个包含n个点,m条边的无向连通图.从顶点1出发,往其余所有点分别走一次并返回. 往某一个点走时,选择总长度最短的路径走.若有多条长度最短的路径,则选择经过的顶点序列字典序最 ...

  2. Bzoj4016/洛谷P2993 [FJOI2014] 最短路径树问题(最短路径问题+长链剖分/点分治)

    题面 Bzoj 洛谷 题解 首先把最短路径树建出来(用\(Dijkstra\),没试过\(SPFA\)\(\leftarrow\)它死了),然后问题就变成了一个关于深度的问题,可以用长链剖分做,所以我 ...

  3. 点分治问题 ----------- P2993 [FJOI2014]最短路径树问题 [最短路径树+点分治+采坑]

    题目链接 解题思路: 首先我们知道最小路径树实际上就是Dijkstra算法在找最短路的时候转移的过程就是一个最短路径树. 那么我们就可以先跑个最短路,记录一下各个最短路的路径.然后就是很裸的点分治.分 ...

  4. Bzoj4016: [FJOI2014]最短路径树问题

    题面 传送门 Sol 先\(SPFA\)求出单源最短路,\(Bfs\)建出树,字典序可以用堆解决 然后就是点分治的一眼题 开桶记录到当前根经过边长度相同的最长路,记录它的长度 自己强行\(yy\)了一 ...

  5. HDU 5385 The path(贪心、构造、最短路径树)

    HDU 5385 题目大意:给定一个图,dis表示第i个点到1点的最短路,dis1=0,给有向图上的边赋权值(1~n)满足dis1<dis2<dis3<--<disk>d ...

  6. CSP认证201609-4 交通规划[C++题解]:最短路径树、dijkstra求单源最短路、递推思想

    题目分析 来源:acwing 分析: 这题是最短路树.保持原图中所有点到根结点的最短距离不变,然后在原图中选择一些边,使所有点连通的最短路是多长. 最短路径树,是一种使用最短路径算法生成的数据结构树. ...

  7. Luogu P1198 BZOJ 1012 最大数 (线段树)

    Luogu P1198 BZOJ 1012 最大数 (线段树) 手动博客搬家: 本文发表于20170821 14:32:05, 原地址https://blog.csdn.net/suncongbo/a ...

  8. POJ3013 Big Christmas Tree(最短路径树)

    题目大概说给一张点和边都有权的图,现在要求其一棵以1结点为根的生成树使树的边权和最小,树边权 = 对应的图边权 * 树边末端点为根的子树所有结点对于图顶点的点权和. 要求∑(边权*子树点权和),等价于 ...

  9. 最小生成树、最短路径树

    一.最小生成树与最短路径树的区别 最小生成树能够保证整个拓扑图的所有路径之和最小,但不能保证任意两点之间是最短路径. 应用如网络部线,把所有的电脑(服务器?)都连起来用的网线(光纤?)最少,即用最小的 ...

最新文章

  1. 基于OpenCV提取特定区域方法汇总
  2. 机器学习模型调参指南(附代码)
  3. 从微信AI首席顾问到金融文档智能
  4. 深度理解特斯拉自动驾驶解决方案 2:向量空间
  5. while循环语句用法_EXCEL退出DO LOOP循环语句|until、while
  6. node 获取表单数据 为空_像声明类型一样写表单——基础功能
  7. 支持的vCenter Server High Availability选项(2096800)
  8. java.sql.SQLException: null, message from server: “Host ‘xxx.xxx.xxx.xxx‘ is not allowed to
  9. android自定义矩形框,Android中ImageCropper矩形、圆形 裁剪框的实现方法
  10. MongoDBTool-最新进展报告
  11. centos7.0 安装docker
  12. H5调用安卓以及IOS前置摄像头
  13. ps 蒙版抠图更换背景色或背景图片
  14. [BZOJ]4987: Tree 树形DP
  15. Fedora 12下使用HP USB打印机(HP laserjet P1007)
  16. 培训班出来的程序员怎么就不行了?
  17. 年轻不要停止奋斗的脚步
  18. python func函数用法_python函数局部变量用法实例分析
  19. Omni Core v0.11.0 rpc-api
  20. Flutter的基础知识之Dart语法

热门文章

  1. php跨域访问java,案例:PHP Ajax 跨域最佳解决方案
  2. redis hash删除所有key_一文详解Redis面试中常见的5种数据结构及对应使用场景
  3. linux install goolepinyin_Linux截图工具推荐(Ubuntu 18.04亲测)
  4. 福建省2013高职单招计算机类试题,13年福建-高职单招-计算机类试题及答案.doc
  5. mysql innodbmaxdirtypagespct_MySQL参数解析innodb_max_dirty_pages_pct
  6. mysql配置文件没有spring_spring cloud config使用mysql存储配置文件
  7. php 上下位置居中,CSS_css中position:fixed实现div在窗口上下左右居中,上下左右 居中 代码如下 复制 - phpStudy...
  8. TCP的协议特点(对于《Linux高性能服务器编程》的部分摘录以及自己的部分见解)
  9. bash 运行程序 下一步_怎样用 Bash 编程:语法和工具
  10. OpenCV中的二进制鲁棒独立基本特征——BRIEF