【spfa】bzoj3921 Mimori与树海

考虑“删除后图仍连通”,即其不是无向图的桥(bridge),可以用Tarjan算法预处理,这里不赘述。

【算法一】

枚举删除的是哪条边,然后枚举起点,暴搜,统计答案。

可以通过0、1号测试点。

预计得分:20分。

【算法二】

枚举删除的是哪条边,接着,

①枚举起点,∵op=0,∴裸最短路,Heap-dijkstra/SPFA。时间复杂度: O(m2*n*logn)/O(m2*n*k)。

②Floyd,时间复杂度O(m*n3)。

可以通过0、1、2、3号测试点。

预计得分:40分。

【算法三】

枚举删的是哪条边,接着,枚举终点,逆着推回去,考虑op=1的边,则有w[U][V]=(dis[U]+18)/19。然后跑Heap-dijkstra/SPFA即可。

时间复杂度:O(m2*n*logn)/O(m2*n*k)。

可以通过0、1、2、3、4、5号测试点。

预计得分:60分。

【算法四】

考虑删除某条边E后,以某个点U为终点的最短路是否发生变化。

我们发现,在某个无向图中,以U为根节点,到其他所有点的最短路构成了一棵树,即“最短路树(shortest-path tree)”。我们删除的边如果不在最短路树上,对当前图中的以U为起点的所有最短路不发生变化。

我们可以在Heap-dijkstra/SPFA的同时处理出最短路树:

if(dis[V]>dis[U]+w[U][V]) TreeEdge[S][V]=Edge[U][V];

//TreeEdge[S][V]表示以S为起点的单源最短路树中,V结点的父边。

于是算法得出:

首先枚举终点,然后枚举删除的边,仅仅需要重复求删除n-1条在最短路树上的边的答案即可,其他的边可以直接用之前预处理的答案就好了。

时间复杂度:O(n2*m*logn)/O(n2*k*m)。Heap-dijkstra可能需要卡常数。

预计得分:100分。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
#define N 111
#define M 1101
queue<int>q;
int dis[N],ans,ans2[M<<1],sumv[N],need[N];
int n,m,first[N],next[M<<1],v[M<<1],en,TEs[N][N]/*树边*/,x,y,es[N][N];
bool w[M<<1]/*道路类型*/,z,Del[M<<1]/*当前删除的边*/,ontree[N][M<<1]
/*ontree[i][j]表示边j是否在以i为汇点的最短路树上*/,inq[N];
void AddEdge(const int &U,const int &V,const bool &W)
{
	v[en]=V; w[en]=W;
	next[en]=first[U];
	first[U]=en;
	es[U][V]=en++;
}
int calc(const int &U,const bool &op){return op ? (dis[U]+18)/19 : 1;}
void spfa(const int &s,const bool &op)
{
    memset(dis,0x7f,sizeof(dis));
    memset(inq,0,sizeof(inq));
    dis[s]=need[s]; inq[s]=1; q.push(s);
    while(!q.empty())
      {
        int U=q.front();
        for(int i=first[U];i!=-1;i=next[i]) if(!Del[i])
          {
          	int C=calc(U,w[i]);
          	if(dis[v[i]]>dis[U]+C)
              {
              	if(!op) TEs[s][v[i]]=es[U][v[i]];
            	dis[v[i]]=dis[U]+C;
            	if(!inq[v[i]])
            	  q.push(v[i]),inq[v[i]]=1;
              }
          }
        q.pop(); inq[U]=0;
      }
    if(!op)
      {
      	for(int i=1;i<=n;++i)
      	  sumv[s]+=dis[i];
      	ans+=sumv[s];
      	for(int i=1;i<=n;++i)
		  if(i!=s)
      	    ontree[s][TEs[s][i]]=ontree[s][TEs[s][i]^1]=1;
      }
}
int dfn[N],dep;
bool bridge[M<<1];
int Tarjan(int U,int fa)
{
    int lowU=dfn[U]=++dep;
    for(int i=first[U];i!=-1;i=next[i])
      if(!dfn[v[i]])
		{
          int lowV=Tarjan(v[i],U);
          lowU=min(lowU,lowV);
          if(lowV>dfn[U])
            bridge[i]=bridge[i^1]=1;
        }
      else if(v[i]!=fa)
	    lowU=min(lowU,dfn[v[i]]);
    return lowU;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;++i) scanf("%d",&need[i]);
	memset(first,-1,sizeof(first));
	for(int i=1;i<=m;++i)
	  {
	  	scanf("%d%d%d",&x,&y,&z);
	  	AddEdge(x,y,z);
	  	AddEdge(y,x,z);
	  }
	Tarjan(1,0);
	for(int i=1;i<=n;++i)
	  spfa(i,0);
	for(int i=1;i<=n;++i)
	  {
	  	for(int j=1;j<=n;++j)
	  	  if((!bridge[TEs[i][j]])&&(j!=i))
	  	    {
	  	      Del[TEs[i][j]]=Del[TEs[i][j]^1]=1;
	  	      spfa(i,1);
	  	      Del[TEs[i][j]]=Del[TEs[i][j]^1]=0;
	  	      int tmp=ans2[TEs[i][j]];
	  	      for(int k=1;k<=n;++k)
	  	        ans2[TEs[i][j]]+=dis[k];
	  	      ans2[TEs[i][j]^1]+=(ans2[TEs[i][j]]-tmp);
	  	    }
	  	for(int j=0;j<en;++j)
		  if((!bridge[j])&&(!ontree[i][j]))
	  	    ans2[j]+=sumv[i];
	  }
	printf("%d %d\n",ans,*max_element(ans2,ans2+en));
	return 0;
}
时间: 2024-10-09 02:31:02

【spfa】bzoj3921 Mimori与树海的相关文章

UESTC30-最短路-Floyd最短路、spfa+链式前向星建图

最短路 Time Limit: 3000/1000MS (Java/Others) Memory Limit: 65535/65535KB (Java/Others) 在每年的校赛里,所有进入决赛的同学都会获得一件很漂亮的T-shirt.但是每当我们的工作人员把上百件的衣服从商店运回到赛场的时候,却是非常累的!所以现在他们想要寻找最短的从商店到赛场的路线,你可以帮助他们吗? Input 输入包括多组数据. 每组数据第一行是两个整数NN ,MM (N≤100N≤100 ,M≤10000M≤1000

畅通project续HDU杭电1874【dijkstra算法 || SPFA】

http://acm.hdu.edu.cn/showproblem.php?pid=1874 Problem Description 某省自从实行了非常多年的畅通project计划后.最终修建了非常多路.只是路多了也不好,每次要从一个城镇到还有一个城镇时,都有很多种道路方案能够选择,而某些方案要比还有一些方案行走的距离要短非常多.这让行人非常困扰. 如今,已知起点和终点,请你计算出要从起点到终点.最短须要行走多少距离. Input 本题目包括多组数据.请处理到文件结束. 每组数据第一行包括两个正

HDU 2722 Here We Go(relians) Again (spfa)

Here We Go(relians) Again Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 32768/32768K (Java/Other) Total Submission(s) : 1   Accepted Submission(s) : 1 Font: Times New Roman | Verdana | Georgia Font Size: ← → Problem Description The Gorelians

[BZOJ 1295][SCOI2009]最长距离(SPFA+暴力)

题目:http://www.lydsy.com/JudgeOnline/problem.php?id=1295 分析:很巧妙的一道spfa从搜索的角度是搜索在所有1中搜索删除哪T个1,对整个图询问,这样肯定TLE 不妨反过来想一想:对于两个点,弄出联通这两个点所需删除的最少的1,那么就知道这两个点是否可以作为题目要求的起点和终点,如果满足算一下结果和ans比较一下就可以. 所以算法就出来了: 枚举起点(S,T),用SPFA跑出图上的所有点到起点这条路径联通的最少删除的1,那么ans=max(di

poj3268 Silver Cow Party (SPFA求最短路)

其实还是从一个x点出发到所有点的最短路问题.来和回只需分别处理一下逆图和原图,两次SPFA就行了. #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<string> #include<cmath> #include<map> #include<set> #include<vector> #

POJ--3259--Wormholes【SPFA判负权值回路】

题意:有n个点,之间有m条双向路径,还有w个虫洞,单向,从一点到另一点需要花费时间,但是有虫洞的话会减少时间,一个人想要走某一条路使得他能碰到过去的自己,问这个图是否能让他实现他的想法. 其实就是判一个图是否存在负权值回路,SPFA可以实现,原理是:如果存在负权值回路,那么从源点到某个顶点的距离就可以无限缩短,因此就会无限入队,所以在SPFA中统计每个顶点的入队次数,如果超过了n个(顶点个数)则说明存在负权值回路. 我把输出yes和输出no写反了,WA了两发,看了半天都没发现... #inclu

POJ 3259 Wormholes SPFA算法题解

本题其实也可以使用SPFA算法来求解的,不过就一个关键点,就是当某个顶点入列的次数超过所有顶点的总数的时候,就可以判断是有负环出现了. SPFA原来也是可以处理负环的. 不过SPFA这种处理负环的方法自然比一般的Bellman Ford算法要慢点了. #include <stdio.h> #include <string.h> #include <limits.h> const int MAX_N = 501; const int MAX_M = 2501; const

ZOJ 3794 Greedy Driver spfa

题意: 给定n个点,m条有向边,邮箱容量. 起点在1,终点在n,开始邮箱满油. 下面m行表示起点终点和这条边的耗油量(就是长度) 再下面给出一个数字m表示有P个加油站,可以免费加满油. 下面一行P个数字表示加油站的点标. 再下面一个整数Q 下面Q行 u v 表示在u点有销售站,可以卖掉邮箱里的任意数量的油,每以单位v元. 问跑到终点能获得最多多少元. 先求个每个点的最大剩余油量 f[i], 再把边反向,求每个点距离终点的最短路 dis[i]. 然后枚举一下每个销售点即可,( f[i] - dis

【最短路】【spfa】hdu6071 Lazy Running

给你一个4个点的环,问你从2号点出发, 再回到2号点,长度>=K的最短路是多少.环上的边长度不超过30000. 跑出来所有dis(2,j)以后,然后for一遍j,根据dis(2,j)+t*2*w>=K,解出来对于每个j而言最小的t,然后尝试更新答案即可.如果dis(2,j)已经大于等于K了,那直接用其尝试更新答案. 跟CDOJ1633很像. #include<cstdio> #include<queue> #include<algorithm> #inclu