BZOJ 4289: PA2012 Tax(最短路)
2024-10-01 14:06:35
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 755 Solved: 240
[Submit][Status][Discuss]
Description
给出一个N个点M条边的无向图,经过一个点的代价是进入和离开这个点的两条边的边权的较大值,求从起点1到点N的最小代价。起点的代价是离开起点的边的边权,终点的代价是进入终点的边的边权
N<=100000
M<=200000
Input
Output
Sample Input
4 5
1 2 5
1 3 2
2 3 1
2 4 4
3 4 8
1 2 5
1 3 2
2 3 1
2 4 4
3 4 8
Sample Output
12
HINT
Source
这题居然卡long long,也是没谁了
首先一个很显然的思路是暴力拆边
即把每个点每一条入边和每一条出边的两两看做一个点,权值为两边的较大值
但是这样很显然是$O(m^2)$,肯定会GG
所以我们考虑一种神仙操作。
对于一条无向边,我们把它看成两条有向边
然后我们这样构图
1.对于一个点,我们把它的出边从小到大排序
2.枚举每一条边,如果这条边连接着1或者N,那么我们从S连向这条边或者从这条边连向T,权值为该边的权值
3.从改边所对应的入边向该边连一条边,边权为它们的权值
4.枚举每一条出边,从权值较小的向权值较大的连权值为两边差值的边,从权值较大的向权值较小的连权值为0的边
可能这样说不是很清楚,借鉴一下这位大佬的图
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
#define Pair pair<long long,int>
#define F first
#define S second
const int MAXN=*1e6+;
using namespace std;
inline int read()
{
char c=getchar();int x=,f=;
while(c<''||c>''){if(c=='-')f=-;c=getchar();}
while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
return x*f;
}
struct Edge
{
int u,v,w,nxt;
}E[MAXN];
int headE[MAXN],numE=;
inline void add_edge(int x,int y,int z)
{
E[numE].u=x;
E[numE].v=y;
E[numE].w=z;
E[numE].nxt=headE[x];
headE[x]=numE++;
}
struct node
{
int u,v,w,nxt;
}edge[MAXN];
int head[MAXN],num=;
inline void AddEdge(int x,int y,int z)
{
edge[num].u=x;
edge[num].v=y;
edge[num].w=z;
edge[num].nxt=head[x];
head[x]=num++;
}
int N,M,S,T;
int temp[MAXN];
long long dis[MAXN];
bool vis[MAXN];
void Dijstra()
{
memset(dis,0xf,sizeof(dis));dis[S]=;
priority_queue<Pair>q;
q.push(make_pair(,S));
while(q.size()!=)
{
while(vis[q.top().second]&&q.size()>) q.pop();
long long p=q.top().second;
vis[p]=;
for(int i=head[p];i!=-;i=edge[i].nxt)
if(dis[edge[i].v]>dis[p]+edge[i].w)
dis[edge[i].v]=dis[p]+edge[i].w,
q.push(make_pair(-dis[edge[i].v],edge[i].v));
}
printf("%lld\n",dis[T]);
}
int comp(const int &a,const int &b)
{
return E[a].w<E[b].w;
}
int main()
{
#ifdef WIN32
freopen("a.in","r",stdin);
#else
#endif
memset(headE,-,sizeof(headE));
memset(head,-,sizeof(head));
N=read();M=read();S=,T=*(M+);
for(int i=;i<=M;i++)
{
int x=read(),y=read(),z=read();
add_edge(x,y,z);
add_edge(y,x,z);
}
for(int i=;i<=N;i++)
{
int tempnum=;
for(int j=headE[i];j!=-;j=E[j].nxt)
temp[++tempnum]=j;
sort(temp+,temp+tempnum+,comp);
for(int j=;j<=tempnum;j++)
{
int x=temp[j],y=temp[j+];
if(E[x].u==)
AddEdge(S,x,E[x].w);
if(E[x].v==N)
AddEdge(x,T,E[x].w);
AddEdge(x^,x,E[x].w);
if(j!=tempnum)
AddEdge(x,y,E[y].w-E[x].w),
AddEdge(y,x,);
}
}
Dijstra();
return ;
}
最新文章
- Ajax1
- CDN系统对网站的性能有极大的提升
- Java笔记5-修饰符,重载,递归,数组
- 【OpenGL】VS2010环境配置 [转]
- paip.自适应网页设计 跟 响应式 设计方法与工具补充(2)o54
- Uploadify v3.2.1 属性、事件、方法说明
- 一个Java方法覆盖的小问题
- shell括号操作符
- C++类型引用浅析
- 全国省市区Json文件 ,做省市区联动很轻松
- go - 内置基础类型
- java控件之树形结构JTree
- [WinForm]委托应用①——窗口之间方法/控件调用
- 基于MMSE的预测
- Variable number of arguments (Varargs)
- Docker Registry V2 Garbage Collection
- 使用 CSS3 的 box-sizing 属性设置元素大小包含 border 与 padding
- Bzoj5332: [Sdoi2018]旧试题
- python 安装pip setuptools
- React Native 组建之IOS和Android通用抽屉
热门文章
- [实例]ROS使用OpenCV读取图像并发布图像消息在rviz中显示
- POJ 2392 DP
- Android local socket学习总结
- 比较两个时间的大小 举例:CompareDate(";12:00";,";11:15";)
- IIS7.0与AP.NET
- js通过经纬度计算两点之间的距离
- DataBaseFactory基础了解
- Ubuntu下快速配置Caffe
- 全文检索lucene6.1的检索方式
- NuSOAP笔记:如何创建复杂数据类型