矩阵乘法(快速幂)

为说明方便,这里让\(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\)矩阵的初始值为无限大。

然后我们就可以用\(O(nk^3)\)的复杂度去做这道题了。但这明显不行。

我们设没有直接连通的两个点距离为无限大,构建出邻接矩阵\(D\),就可以魔改一下上面的式子,改成:

\[T[i][j]=min(T[i][x]+D[x][j])
\]

其中\(x\)为自己枚举的中间节点,然后就出现的如下的代码:

for(int i=0;i<k;++i){
for(int j=0;j<k;++j){
for(int l=0;l<k;++l){
ret.a[i][j]=min(ret.a[i][j],a.a[i][l]+b.a[l][j]);
}
}
}

发现,这不是就是矩阵乘法吗?

因为取最小值满足可加性,所以使用矩阵快速幂是可行的。这样,我们就能把复杂度优化为\(O(lognk^3)\)

然后,我就不开O2过不了了。

我们发现从源点能到达的点数最多只有\(k+1\)(因为即使走过每条边都发现一个新节点,也只能发现这么多点。)所以我们可以只用源点能到的点进行离散化,可以将点数从\(2k\)到\(k\),从而在矩阵乘法时省掉8倍常数,然后就可以不开O2AC了。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1000010,M=200;
int n,k,s,t;
struct data{
ll a[M][M];
data(){memset(a,0,sizeof a);}
}a;
data operator*(const data&a,const data&b){
data ret;
memset(ret.a,0x3f,sizeof ret.a);
for(int i=0;i<k;++i){
for(int j=0;j<k;++j){
for(int l=0;l<k;++l){
ret.a[i][j]=min(ret.a[i][j],a.a[i][l]+b.a[l][j]);
}
}
}
return ret;
}
data mpow(data a,int n){
data ret=a;
n--;
while(n){
if(n&1)ret=ret*a;
n/=2;
a=a*a;
}
return ret;
}
int tot,bian[N],nxt[N],head[N];
void add(int x,int y){
tot++,bian[tot]=y,nxt[tot]=head[x],head[x]=tot;
}
struct edge{
int x,y;
ll z;
}e[M];
int vis[N];
vector<int>v;
void dfs(int x){
if(vis[x])return;
v.push_back(x);
vis[x]=1;
for(int i=head[x];i;i=nxt[i]){
int y=bian[i];
dfs(y);
}
}
int main(){
cin>>n>>k>>s>>t;
memset(a.a,0x3f,sizeof a.a);
for(int i=1;i<=k;++i){
scanf("%lld%d%d",&e[i].z,&e[i].x,&e[i].y);
add(e[i].x,e[i].y);
add(e[i].y,e[i].x);
}
dfs(s);
sort(v.begin(),v.end());
for(int i=1;i<=k;++i){
if(!vis[e[i].x])continue;
int x=lower_bound(v.begin(),v.end(),e[i].x)-v.begin(),
y=lower_bound(v.begin(),v.end(),e[i].y)-v.begin();
a.a[y][x]=a.a[x][y]=min(a.a[x][y],e[i].z); }
data ret=mpow(a,n);
s=lower_bound(v.begin(),v.end(),s)-v.begin();
t=lower_bound(v.begin(),v.end(),t)-v.begin();
cout<<ret.a[s][t]<<endl;
}

zjoj1706: [usaco2007 Nov]relays 奶牛接力跑的更多相关文章

  1. BZOJ_[usaco2007 Nov]relays 奶牛接力跑_离散化+倍增弗洛伊德

    BZOJ_[usaco2007 Nov]relays 奶牛接力跑_离散化+倍增弗洛伊德 Description FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们 ...

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

    [BZOJ1706][usaco2007 Nov]relays 奶牛接力跑 Description FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项 ...

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

    Description FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项目.至于进行接力跑的地点 自然是在牧场中现有的T(2 <= T < ...

  4. bzoj 1706: [usaco2007 Nov]relays 奶牛接力跑——倍增floyd

    Description FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项目.至于进行接力跑的地点 自然是在牧场中现有的T(2 <= T < ...

  5. 【bzoj1706】[usaco2007 Nov]relays 奶牛接力跑 离散化+倍增Floyd

    题目描述 FJ的N(2 <= N <= 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项目.至于进行接力跑的地点 自然是在牧场中现有的T(2 <= T <= 100) ...

  6. 【BZOJ】1706: [usaco2007 Nov]relays 奶牛接力跑

    [题意]给定m条边的无向图,起点s,终点t,要求找出s到t恰好经过n条边的最短路径.n<=10^6,m<=100. [算法]floyd+矩阵快速幂 [题解] 先对点离散化,得到点数N. 对 ...

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

    题目大意:有t(t<=100)条无向边连接两点,求s到e刚好经过n(n<=10^7)条路径的最小距离. 第一反应分层图,但是一看n就懵逼了,不会写.看了题解之后才知道可以这么玩... 首先 ...

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

    大概是叫倍增Floyd? 显然最多200个点...f[i][j][k]表示从j到k,走2^i步的最小路程.就随便转移了.. 查询的话就是把n二进制位上是1的那些都并起来. #include<cs ...

  9. 【bzoj1706】[usaco2007 Nov]relays 奶牛接力跑

    题意 给出一张无向图,求出恰巧经过n条边的最短路. 题解 考虑先离散化,那么点的个数只会有202个最多.于是复杂度里面就可以有一个\(n^3\).考虑构造矩阵\(d^1\)表示经过一条边的最短路,那么 ...

随机推荐

  1. springboot整合mybatis,mongodb,redis

    springboot整合常用的第三方框架,mybatis,mongodb,redis mybatis,采用xml编写sql语句 mongodb,对MongoTemplate进行了封装 redis,对r ...

  2. mysql中数据表记录的增删查改(2)

    select `数据表.字段1`, group_concat(`数据表.字段2`) from `数据表` group by `数据表.字段1` order by `数据表.字段1` desc; sel ...

  3. 通过names.index()方法找到第2个eva值 ,并将其改成EVA

    names= ['alex','rain','peiqi','eva','mac','jack','eva','kangkang','jain']first_index=names.index('ev ...

  4. Debian 下忘记root密码的特殊修改方式

    当系统开机进入 grub页面时,在启动条目上按下键盘的 e 进入编辑 找到 linux 开头的一行,类似下面这样 linux /boot/vmlinux-4.9.0.8-amd64 root=/dev ...

  5. 【转】webpack4安装过程遇到的问题及处理方法

    随便百度一下,安装使用webpack的教程铺天盖地,安装一步步来,最后的最后打包没反应......,浪费了不少的时间. 这里我要提醒一下,如果安装webpack1,2,3按照百度上的教程应该不会有问题 ...

  6. java之struts2之类型转换

    在使用servlet开发中,表单中提交的数据到servlet后都是字符串类型,需要程序员手动进行类型转换. 但是到struts2后,基本数据类型struts2都可以转换.但是如果是自定义类型,stru ...

  7. 一张图看懂SharpImage

    通过下面的图片,可以瞬间看懂整个类库的脉络.图片比较大,如果看不清,可

  8. showModalDialog的使用方法

    基本介绍: showModalDialog()         (IE 4+ 支持) showModelessDialog()      (IE 5+ 支持) window.showModalDial ...

  9. springboot笔记05——profile多环境配置切换

    前言 一个应用程序从开发到上线,往往需要经历几个阶段,例如开发.测试.上线.每个阶段所用到的环境的配置可能都是不一样的,Springboot 应用可以很方便地在各个环境中对配置进行切换.所以,今天主要 ...

  10. 5_PHP数组_3_数组处理函数及其应用_7_数组排列函数

    以下为学习孔祥盛主编的<PHP编程基础与实例教程>(第二版)所做的笔记. 数组排列函数 1. sort() 函数 程序: <?php $array = array("img ...