首页 > 代码库 > 【POJ3613】【USACO 2007 Nov Gold】 2.Cow Relays 矩阵乘法?
【POJ3613】【USACO 2007 Nov Gold】 2.Cow Relays 矩阵乘法?
题意:给你一个m条边的图,求s到t的正好用k条边的最短路。
(输入k,t,s,t)
题解:
先说说暴力。
动规f[k][i][j]表示i到j经过k条边的最短路,然后外层循环k一遍遍跑最后出解。
显然大概率T。
然后有一种思路:
我们可以动规求得f[k][i][j]表示i到j经过k条边的最短路,然后再求g[i]表示从终点走i步回到终点的最短路。
这样我们就可以乱搞过了。(没写,也没细想)
再之后是正解:
我们可以利用类似于快速幂的方法求f[i][j]表示i到j正好用多少多少步。
然后思想是使用边数的累加,或者叫倍增。
贴代码:
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define N 205 #define M 1005 #define inf 0x3f3f3f3f using namespace std; int n,m,s,t; int hash[M],dis[N][N],temp[N][N],f[N][N]; void power_floyd(int K) { int i,j,k; while(K) { if(K&1) { memset(temp,0x3f,sizeof temp); for(k=1;k<=n;k++)for(i=1;i<=n;i++)for(j=1;j<=n;j++)temp[i][j]=min(temp[i][j],f[i][k]+dis[k][j]); for(i=1;i<=n;i++)for(j=1;j<=n;j++)f[i][j]=temp[i][j]; } memset(temp,0x3f,sizeof temp); for(k=1;k<=n;k++)for(i=1;i<=n;i++)for(j=1;j<=n;j++)temp[i][j]=min(temp[i][j],dis[i][k]+dis[k][j]); for(i=1;i<=n;i++)for(j=1;j<=n;j++)dis[i][j]=temp[i][j]; K>>=1; } return ; } int main() { // freopen("test.in","r",stdin); int i,j,k; int a,b,c; memset(f,0x3f,sizeof f); memset(dis,0x3f,sizeof dis); for(i=1;i<=200;i++)f[i][i]=0; scanf("%d%d%d%d",&k,&m,&s,&t); while(m--) { scanf("%d%d%d",&c,&a,&b); if(!hash[a])hash[a]=++n; if(!hash[b])hash[b]=++n; dis[hash[a]][hash[b]]=dis[hash[b]][hash[a]]=c; } s=hash[s],t=hash[t]; power_floyd(k); printf("%d\n",f[s][t]); return 0; }
【POJ3613】【USACO 2007 Nov Gold】 2.Cow Relays 矩阵乘法?
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。