bzoj1706: [Usaco2007 Nov]relays 奶牛接力跑 (Floyd+新姿势)

题目大意:有t(t<=100)条无向边连接两点,求s到e刚好经过n(n<=10^7)条路径的最小距离。

第一反应分层图,但是一看n就懵逼了,不会写。看了题解之后才知道可以这么玩。。。

首先有100条边最多200个点,但点编号到1000,所以离散化一下。

任何一个正整数都能用2的幂相加得到,所以先把n转变成2进制来看,按位考虑。dist[i][j][k]表示刚好经过2^i条边从j到k的最短距离,则dist[i,j,k]=min{dist[i-1][j][l]+dist[i-1][l][k]}。用类似快速幂的方法,若是2进制的n这一位是1的话,就把答案数组g和dist跑floyd(min{g[i-1][j][l]+dist[i-1][l][k]}),否则dist自己和自己跑(min{dist[i-1][j][l]+dist[i-1][l][k]}),这样子g就是跑n条路径得出的最小距离了。

代码如下:

type
  map=array[1..1000,1..1000]of longint;
var
  n,t,s,e,i,j,k,len,x,y,tot:longint;
  c,d,g:array[1..1000,1..1000]of longint;
  num:array[1..1000]of longint;

procedure merge(var a,b:map);
var
  i,j,k:longint;
begin
  fillchar(c,sizeof(c),63);
  for k:=1 to tot do
  for i:=1 to tot do
  for j:=1 to tot do
  if c[i,j]>a[i,k]+b[k,j] then
  c[i,j]:=a[i,k]+b[k,j];
  a:=c;
end;

procedure work;
begin
  fillchar(g,sizeof(g),63);
  for i:=1 to tot do g[i,i]:=0;
  while n>0 do
  begin
    if n and 1=1 then merge(g,d);
    merge(d,d);
    n:=n>>1;
  end;
  writeln(g[num[s],num[e]]);
end;

begin
  readln(n,t,s,e);
  fillchar(d,sizeof(d),63);
  for i:=1 to t do
  begin
    readln(len,x,y);
    if num[x]=0 then
    begin
      inc(tot);num[x]:=tot;
    end;
    if num[y]=0 then
    begin
      inc(tot);num[y]:=tot;
    end;
    d[num[x],num[y]]:=len;d[num[y],num[x]]:=len;
  end;
  work;
end.

时间: 2024-10-21 21:00:43

bzoj1706: [Usaco2007 Nov]relays 奶牛接力跑 (Floyd+新姿势)的相关文章

[bzoj1706] [usaco2007 Nov]relays 奶牛接力跑

大概是叫倍增Floyd? 显然最多200个点...f[i][j][k]表示从j到k,走2^i步的最小路程.就随便转移了.. 查询的话就是把n二进制位上是1的那些都并起来. 1 #include<cstdio> 2 #include<iostream> 3 #include<cstring> 4 #include<algorithm> 5 #include<cstdlib> 6 #define ll long long 7 #define ull

[BZOJ] 1706: [usaco2007 Nov]relays 奶牛接力跑

1706: [usaco2007 Nov]relays 奶牛接力跑 Time Limit: 5 Sec  Memory Limit: 64 MBSubmit: 707  Solved: 367[Submit][Status][Discuss] Description FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项目.至于进行接力跑的地点 自然是在牧场中现有的T(2 <= T <= 100)条跑道上. 农场上的跑道有一些交汇点,每条跑道都连结

矩阵乘法专题2——bzoj 1706 [usaco2007 Nov] relays 奶牛接力跑 题解

转载请注明:http://blog.csdn.net/jiangshibiao/article/details/24960651 [原题] 1706: [usaco2007 Nov]relays 奶牛接力跑 Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 340  Solved: 162 [Submit][Status] Description FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项目.至于进行接力

BZOJ 1706 usaco2007 Nov relays 奶牛接力跑 倍增Floyd

题目大意:给定一张无向图,求从s出发恰好经过n条边到达e的最短路 倍增Floyd--为何大家都管这个叫做矩阵乘法- - 算了为何要纠结这种事- - 令f[p][i][j]表示走2^p步从i到达j的最短路 有f[p][i][j]=min{f[p-1][i][k]+f[p-1][k][j]} 将n进行二进制拆分 用矩阵g记录答案矩阵 对于每一位p 用f[p]和g两个矩阵搞出h 再将h的值赋给g 切忌直接用f[p]更新g 这样可能导致g的上一次的值被继承到下一次里 这样相当于这一步没跑 100条边,因

bzoj 1706: [usaco2007 Nov]relays 奶牛接力跑【矩阵乘法+Floyd】

唔不知道怎么说--大概核心是把矩阵快速幂的乘法部分变成了Floyd一样的东西,非常之神 首先把点离散一下,最多有200个,然后建立邻接矩阵,a[u][v]为(u,v)之间的距离,没路就是inf 然后注意重载乘号的部分,注意到这样乘一次就相当于把本来存了经过k条路的最短路的邻接矩阵变成存了经过k+1条路的最短路的邻接矩阵 这样看来乘n次就行了,这里用矩阵快速幂 妙啊 #include<iostream> #include<cstdio> using namespace std; co

zjoj1706: [usaco2007 Nov]relays 奶牛接力跑

矩阵乘法(快速幂) 为说明方便,这里让\(k\)为点数,\(n\)为路径长度. 先将点都离散化,这样最后的点只有\(2k\)个. 先考虑一种暴力,每次用\(O(k^3)\)的复杂度来暴力更新,设当前长度\(l\)点的两两最短路矩阵为\(S\),现在要转移到\(l+1\)时的最短路矩阵\(T\).我们考虑用每条边更新,对于某条从\(x\)连向\(y\)的长度为\(z\)的边,对于任一点\(i\),有: \[T[i][y]=min(T[i][y],T[i][x]+z)\] 另外,每次更新时,\(T\

【Floyd矩阵乘法】BZOJ1706- [usaco2007 Nov]relays 奶牛接力跑

[题目大意] 给出一张无向图,求出恰巧经过n条边的最短路. [思路] 首先题目中只有100条边,却给出了10000个点(实际上最多只能有200个),离散化一下. 后面就是Floyd的新姿势,以前看过的集训队论文里面有:D 一开始的邻接矩阵是经过一条边的最短路,把这个邻接矩阵记作f[0] f[1]=f[0]*f[0]=f[0]^2(这里的乘法是矩阵乘法),就可以表示恰巧经过两条边的啦. f[2]=f[1]*f[0]=f[0]^3,恰巧表示经过两条边. -- 所以恰巧经过n条边的最短路是f[n-1]

BZOJ 1706 usaco 2007 Nov relays 奶牛接力跑/POJ 3613 Cow Relays 倍增Floyd

题目大意:求恰好走k步从S到T的最短路. 思路:设f[p][i][j]为从i到j恰好走2^p步的最短路,DP方程十分简单:f[p][i][j] = min(f[p][i][j],f[p - 1][i][k] + f[p - 1][k][j]); 对总步数T进行二进制拆分,在T有1的位置上,假如这个位置为p,那么就用f[p][][]来更新答案g[][],最后得到的g[][]就是答案矩阵. 注意要离散化一下.. CODE: #include <cstdio> #include <cstrin

【BZOJ】【1046】/【POJ】【3613】【USACO 2007 Nov】Cow Relays 奶牛接力跑

倍增+Floyd 题解:http://www.cnblogs.com/lmnx/archive/2012/05/03/2481217.html 神题啊= =Floyd真是博大精深…… 题目大意为求S到E,恰好经过N条边的最短路径(姑且称为路径吧,虽然好像已经不是了……) 总共只有大约200个点(很多点根本没走到,离散化一下即可)所以可以考虑Floyd算最短路. 引用下题解: 题目求i,j之间边数恰为N的最短路径(边可以重复走),我们知道线性代数中有:01邻接矩阵A的K次方C=A^K,C[i][j