P2685 [TJOI2012]桥
xcxcli题解
下面思路仿照上述题解,代码基本照抄上述题解
u⇝vu\leadsto vu⇝v表示uuu到vvv的最短路
u→vu\to vu→v表示uuu和vvv直接相连的边
d1ud1_ud1u表示1⇝u1\leadsto u1⇝u的最短路
dnudn_udnu表示n⇝vn\leadsto vn⇝v的最短路
题意化简一下就是让你求删除一条边使得1⇝n1\leadsto n1⇝n最大化,并求删边方案数。
首先不难发现,如果我们选择删除的边不在1⇝n1\leadsto n1⇝n的路径上,并不会对1⇝n1\leadsto n1⇝n的大小产生影响,也就是如果要产生影响必须删除该路径上的一条边。
不妨假设1⇝n1\leadsto n1⇝n路径上的点如下
1→⋯→a→b→c→d→⋯→n1\to \dots\to a\to b{\color{red}\to} c \to d \to \dots \to n1→⋯→a→b→c→d→⋯→n
如果我们现在选择删除b→cb\to cb→c这条边,发现这条路径被切成两部分1→⋯→a→b1\to \dots\to a\to b1→⋯→a→b和c→d→⋯→nc \to d \to \dots \to nc→d→⋯→n,由于我们现在已经删除b→cb\to cb→c这条边,新的1⇝n1\leadsto n1⇝n一定是1⇝x⇝u→v⇝y⇝n1\leadsto x\leadsto u\to v\leadsto y\leadsto n1⇝x⇝u→v⇝y⇝n,并且x∈1→⋯→a→b,y∈c→d→⋯→nx\in 1\to \dots\to a\to b,y\in c \to d \to \dots \to nx∈1→⋯→a→b,y∈c→d→⋯→n
也就是我们现在让其强制走一条不在1⇝n1\leadsto n1⇝n上的路径即u→vu\to vu→v,并且走此路径能够避开b→cb\to cb→c,1⇝x⇝u→v⇝y⇝n1\leadsto x\leadsto u\to v\leadsto y\leadsto n1⇝x⇝u→v⇝y⇝n路径的权值显然是d1u+wu→v+dnvd1_u+w_{u\to v}+dn_vd1u+wu→v+dnv,让其更新删除b→cb\to cb→c后1⇝n1\leadsto n1⇝n的最短路。
显然我们不能枚举x,yx,yx,y然后在枚举u→vu \to vu→v去更新删除b→cb\to cb→c的结果,有一个直接的想法就是1⇝x⇝u→v⇝y⇝n1\leadsto x\leadsto u\to v\leadsto y\leadsto n1⇝x⇝u→v⇝y⇝n路径能够避开哪些1→⋯→a→b→c→d→⋯→n1\to \dots\to a\to b{\color{red}\to} c \to d \to \dots \to n1→⋯→a→b→c→d→⋯→n的边,假设它能够避开a→b→c→da\to b\to c\to da→b→c→d显然这个结果能够更新删除a→ba\to ba→b或者b→cb\to cb→c或者c→dc\to dc→d的答案。
于是我们只需要枚举u→vu\to vu→v这条边(需要保证此边不在1⇝n1\leadsto n1⇝n的路径上),然后找到1⇝x⇝u→v⇝y⇝n1\leadsto x\leadsto u\to v\leadsto y\leadsto n1⇝x⇝u→v⇝y⇝n能够避开哪些边(仔细想想可知避开的路径一定是连续的),进行区间更新(线段树)即可。
这也是上述博客预处理LuL_uLu和RuR_uRu的原因!!!
有了上述两个数组那么1⇝x⇝u→v⇝y⇝n1\leadsto x\leadsto u\to v\leadsto y\leadsto n1⇝x⇝u→v⇝y⇝n避开的路径则是Lu⇝RvL_u \leadsto R_vLu⇝Rv之间的边。
dengyaotriangle详细证明
#include<queue>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
using pii=pair<int,int>;
constexpr int N=100010,M=400010;
constexpr int INF=0x3f3f3f3f;
int h[N],e[M],ne[M],w[M],idx;
void add(int a,int b,int c){e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;}
int n,m;
bool st[N];
int d1[N],dn[N];
int path[N],cnt,id[N];
bool in[M];
int L[N],R[N],ans[N];
void dijkstra(int S,int d[])
{memset(d,0x3f,(n+1)*sizeof(int));memset(st,0,sizeof st);d[S]=0;priority_queue<pii,vector<pii>,greater<pii> >q;q.push({0,S});while(q.size()){int u=q.top().second;q.pop();if(st[u]) continue;st[u]=1;for(int i=h[u];i!=-1;i=ne[i]){int v=e[i];if(d[v]<=d[u]+w[i]) continue;d[v]=d[u]+w[i];q.push({d[v],v});}}
}
void bfs(int k,int d[],int f[])
{queue<int> q;q.push(path[k]);f[path[k]]=k;while(q.size()){int u=q.front();q.pop();for(int i=h[u];i!=-1;i=ne[i]){int v=e[i];if(!id[v]&&!f[v]&&d[u]+w[i]==d[v]) f[v]=k,q.push(v);}}
}
struct node
{int l,r,v;
}tree[N<<2];
void build(int u,int l,int r)
{tree[u]={l,r,INF};if(l==r) return;int mid=l+r>>1;build(u<<1,l,mid),build(u<<1|1,mid+1,r);
}
void modify(int u,int l,int r,int x)
{if(tree[u].l>=l&&tree[u].r<=r) return tree[u].v=min(tree[u].v,x),void();int mid=tree[u].l+tree[u].r>>1;if(l<=mid) modify(u<<1,l,r,x);if(r>mid) modify(u<<1|1,l,r,x);
}
void pushdown(int u)
{if(tree[u].l==tree[u].r){ans[tree[u].l]=tree[u].v;return;}tree[u<<1].v=min(tree[u<<1].v,tree[u].v);tree[u<<1|1].v=min(tree[u<<1|1].v,tree[u].v);int mid=tree[u].l+tree[u].r>>1;pushdown(u<<1),pushdown(u<<1|1);}
int main()
{cin>>n>>m;memset(h,-1,sizeof h);for(int i=1;i<=m;i++){int a,b,c;cin>>a>>b>>c;add(a,b,c),add(b,a,c);}dijkstra(1,d1),dijkstra(n,dn);// 预处理1->n路径上的点int u=1;while(u!=n){path[++cnt]=u;id[u]=cnt;for(int i=h[u];i!=-1;i=ne[i]){int v=e[i];if(dn[v]+w[i]==dn[u]) {in[i]=1;u=v;break;}}}path[++cnt]=n;id[n]=cnt;//1->n路径上的点 for(int i=1;i<=cnt;i++) bfs(i,d1,L),bfs(i,dn,R);--cnt;build(1,1,cnt);//我们将1->n路径上的边编号for(int u=1;u<=n;u++)for(int i=h[u];i!=-1;i=ne[i]){int v=e[i];if(!in[i]&&L[u]<R[v])// L[u]->R[v]边的编号是L[u]->R[v]-1modify(1,L[u],R[v]-1,d1[u]+w[i]+dn[v]);}pushdown(1);int mx=0,res=0;for(int i=1;i<=cnt;i++){if(mx<ans[i]) mx=ans[i],res=1;else if(mx==ans[i]) res++;}if(mx==d1[n]) res=m;cout<<mx<<' '<<res<<'\n';
}