【BZOJ4016】[FJOI2014]最短路径树问题(点分治,最短路)
【BZOJ4016】[FJOI2014]最短路径树问题(点分治,最短路)
题面
题解
首先把最短路径树给构建出来,然后直接点分治就行了。
这个东西似乎也可以长链剖分,然而没有必要。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define MAX 30300
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int n,m,K;
struct Line{int v,next,w;}e[MAX<<2];
int h[MAX],cnt=1;
inline void Add(int u,int v,int w){e[cnt]=(Line){v,h[u],w};h[u]=cnt++;}
vector<pair<int,int> > E[MAX<<1];
int dis[MAX];
bool vis[MAX];
void Dijkstra()
{
priority_queue<pair<int,int> >Q;
memset(dis,63,sizeof(dis));
Q.push(make_pair(0,1));dis[1]=0;
while(!Q.empty())
{
int u=Q.top().second;Q.pop();
if(vis[u])continue;vis[u]=true;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;
if(dis[v]>dis[u]+e[i].w)
{
dis[v]=dis[u]+e[i].w;
Q.push(make_pair(-dis[v],v));
}
}
}
}
void Build()
{
priority_queue<int,vector<int>,greater<int> > Q;Q.push(1);
memset(vis,0,sizeof(vis));
while(!Q.empty())
{
int u=Q.top();Q.pop();vis[u]=true;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(vis[v])continue;
if(dis[v]==dis[u]+e[i].w)
{
E[u].push_back(make_pair(v,e[i].w));
E[v].push_back(make_pair(u,e[i].w));
if(!vis[v])Q.push(v),vis[v]=true;
}
}
}
}
int rt,Size,mx,sz[MAX];
void getroot(int u,int ff)
{
sz[u]=1;int ret=0;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(v==ff||vis[v])continue;
getroot(v,u);sz[u]+=sz[v];
ret=max(ret,sz[v]);
}
ret=max(ret,Size-sz[u]);
if(mx>ret)mx=ret,rt=u;
}
int len[MAX],num[MAX],mxl;
int Len[MAX],Num[MAX],Mxl;
int av,as;
void chk(int v,int s)
{
if(av<v)av=v,as=0;
if(av==v)as+=s;
}
void dfs(int u,int ff,int l,int v)
{
if(l>K)return;mxl=max(mxl,l);
if(len[l]==v)++num[l];
else if(len[l]<v)len[l]=v,num[l]=1;
if(l==K)chk(v,1);
else if(l<K&&Len[K-l])chk(Len[K-l]+v,Num[K-l]);
for(int i=h[u];i;i=e[i].next)
{
int vv=e[i].v;if(vv==ff||vis[vv])continue;
dfs(vv,u,l+1,v+e[i].w);
}
}
void Divide(int u)
{
vis[u]=true;Mxl=0;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(vis[v])continue;
mxl=0;dfs(v,u,1,e[i].w);
for(int j=1;j<=mxl;++j)
{
if(Len[j]==len[j])Num[j]+=num[j];
else if(Len[j]<len[j])Len[j]=len[j],Num[j]=num[j];
num[j]=len[j]=0;
}
Mxl=max(Mxl,mxl);
}
for(int i=1;i<=Mxl;++i)Len[i]=Num[i]=0;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(vis[v])continue;
Size=mx=sz[v];getroot(v,u);
Divide(rt);
}
}
int main()
{
n=read();m=read();K=read()-1;
for(int i=1;i<=m;++i)
{
int u=read(),v=read(),w=read();
Add(u,v,w);Add(v,u,w);
}
Dijkstra();
Build();
memset(h,0,sizeof(h));cnt=1;
memset(vis,0,sizeof(vis));
for(int u=1;u<=n;++u)
for(int i=0,l=E[u].size();i<l;++i)
Add(u,E[u][i].first,E[u][i].second);
Size=mx=n;getroot(1,0);
Divide(rt);
printf("%d %d\n",av,as);
return 0;
}
【BZOJ4016】[FJOI2014]最短路径树问题(点分治,最短路)的更多相关文章
- P2993 [FJOI2014]最短路径树问题 点分治+最短路
这道题还是非常简单的,由于我们要保证最小字典序,因此我们需要把边进行排序,然后从大到小插入,因为链式前向星是倒着存的.我们只需要先跑一个最短路,然后查询边是不是在最短路上,这个可以通过枚举边并用 di ...
- [BZOJ4016][FJOI2014]最短路径树问题
[BZOJ4016][FJOI2014]最短路径树问题 试题描述 给一个包含n个点,m条边的无向连通图.从顶点1出发,往其余所有点分别走一次并返回. 往某一个点走时,选择总长度最短的路径走.若有多条长 ...
- [BZOJ4016][FJOI2014]最短路径树问题(dijkstra+点分治)
4016: [FJOI2014]最短路径树问题 Time Limit: 5 Sec Memory Limit: 512 MBSubmit: 1796 Solved: 625[Submit][Sta ...
- 【BZOJ4016】[FJOI2014]最短路径树问题 最短路径树+点分治
[BZOJ4016][FJOI2014]最短路径树问题 Description 给一个包含n个点,m条边的无向连通图.从顶点1出发,往其余所有点分别走一次并返回. 往某一个点走时,选择总长度最短的路径 ...
- 【BZOJ4016】[FJOI2014]最短路径树问题
[BZOJ4016][FJOI2014]最短路径树问题 题面 bzoj 洛谷 题解 虽然调了蛮久,但是思路还是蛮简单的2333 把最短路径树构出来,然后点分治就好啦 ps:如果树构萎了,这组数据可以卡 ...
- bzoj 4016 [FJOI2014]最短路径树问题(最短路径树+树分治)
4016: [FJOI2014]最短路径树问题 Time Limit: 5 Sec Memory Limit: 512 MBSubmit: 426 Solved: 147[Submit][Stat ...
- BZOJ_4016_[FJOI2014]最短路径树问题_最短路+点分治
BZOJ_4016_[FJOI2014]最短路径树问题_最短路+点分治 Description 给一个包含n个点,m条边的无向连通图.从顶点1出发,往其余所有点分别走一次并返回. 往某一个点走时,选择 ...
- 【BZOJ-4016】最短路径树问题 Dijkstra + 点分治
4016: [FJOI2014]最短路径树问题 Time Limit: 5 Sec Memory Limit: 512 MBSubmit: 1092 Solved: 383[Submit][Sta ...
- [FJOI2014]最短路径树问题 长链剖分
[FJOI2014]最短路径树问题 LG传送门 B站传送门 长链剖分练手好题. 如果你还不会长链剖分的基本操作,可以看看我的总结. 这题本来出的很没水平,就是dijkstra(反正我是不用SPFA)的 ...
随机推荐
- 扫描不同域下的AD账户进行删除
public ResultModel GetEntryOneToDel(string sAMAccountName) { bool del=false; ResultModel result = ne ...
- Windows+Qt使用gRPC
上篇文章<Windows+VS2017使用gRPC>编译出了Windows下可用的gRPC静态lib库文件,在此基础上要想在Qt上使用,需要使用MSVC2017 64bit构建组件进行构建 ...
- maven 学习---Maven添加远程仓库
默认情况下,Maven从Maven中央仓库下载所有依赖关系.但是,有些库丢失在中央存储库,只有在Java.net或JBoss的储存库远程仓库中能找到. 1. Java.net资源库 添加Java.ne ...
- Ubuntu 18.04上安装 phpMyAdmin
我们将安装 phpMyAdmin 在 Ubuntu18.04 上配合 Apache 一起工作. 在安装 phpMyAdmin 之前需要已经安装了LAMP栈并提供了web页面. 如果没有安装可以参照 U ...
- 2016 ACM-ICPC Asia Regional Dalian Online HDU 5875 Function(线段树)
题意 求区间l~r的a[l]%a[l+1]%--%a[r]的值 思路 因为取模的变化是很快的,所以线段树查找区间内第一个小于等于a[l]的数的位置,更新ans后继续查找即可. 注意查询满足某种条件的位 ...
- Mysql 数据库部署(解压免安装)
下载压缩包,解压. 在根目录下添加my.ini, 复制一下内容,保存. [mysqld] skip-grant-tables #绑定IPv4和3306端口 bind-address = 0.0.0.0 ...
- HDU4815 Little Tiger vs. Deep Monkey——0-1背包
题目描述 对于n道题目,每道题目有一个分值,答对加分,答错不得分,你要和一个叫深猴的比赛,题目你可以假设成判断题(不是对就是错),深猴对于所有的题目都是随机选择一个答案,而你是有脑子的,求为了不输掉比 ...
- vue-cli 3 和 vue-cli 2的区别
分享文章: 浅谈vue-cli 3 和 vue-cli 2的区别 https://blog.csdn.net/weixin_42080056/article/details/81631661 vue- ...
- 【Spring AOP】AOP的实现(三)
一.Spring 对AOP的支持 Spring中AOP代理由Spring的IOC容器负责生成.管理,其依赖关系也由IOC容器负责管理.因此,AOP代理可以直接使用容器中的其它bean实例作为目标,这种 ...
- LeetCode 654. Maximum Binary Tree最大二叉树 (C++)
题目: Given an integer array with no duplicates. A maximum tree building on this array is defined as f ...