Cow Relays POJ - 3613

Posted MeowMeowMeow

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Cow Relays POJ - 3613相关的知识,希望对你有一定的参考价值。

这题面是外星人写的吧

定义一个群,其元素为一个矩阵,定义一个该群的元素M,并且M存着一个图。

定义一个二元运算符*,运算结果仍为该群的元素。

如M*M=R,则R(i,j)=Min(Rij,Mik+Mkj)

这个操作满足交换律和结合律。

其意义代表不言而喻。(手动滑稽)

因此可利用快速幂求解。

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
const int MOD=1000;
const int sz=201;
int N,t,s,e,n;
struct mat{
    int a[sz][sz];
    mat(){memset(a,0,sizeof a);}
    void initE(){for(int i=0;i<n;i++)a[i][i]=1;}
    void initfib(){a[1][0]=a[0][1]=a[0][0]=1;a[1][1]=0;}
    mat operator*(const mat&t)const{
        mat q;memset(q.a,0x3f,sizeof q.a);
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
            for(int k=0;k<n;k++)
                q.a[i][j]=min(q.a[i][j],a[i][k]+t.a[k][j]);
        return q;
    }
    void out(){
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++)printf("%d ",a[i][j]);printf("\n");
        }
    }
    mat operator^(int k)const{
        mat res=*this;
        mat m=*this;
        //res.out();m.out();
        while(k>0){
            if(k&1)res=res*m;
            m=m*m;k>>=1;
        }
        return res;
    }

};
int mp[1010];
int main(){ 
    scanf("%d%d%d%d",&N,&t,&s,&e);
    n=0;
    for(int i=0;i<1010;i++)mp[i]=-1;
    mat m;memset(m.a,0x3f,sizeof m.a);
    while(t--){
        int l,x,y;
        scanf("%d%d%d",&l,&x,&y);
        int sx=mp[x],sy=mp[y];
        if(sx==-1)sx=mp[x]=n++;
        if(sy==-1)sy=mp[y]=n++;
        m.a[sx][sy]=m.a[sy][sx]=l;
    }
   // m.out();
    m=m^(N-1);
  //  m.out();
    printf("%d\n",m.a[mp[s]][mp[e]]);
    return 0;
}

 

以上是关于Cow Relays POJ - 3613的主要内容,如果未能解决你的问题,请参考以下文章

poj3613Cow Relays——k边最短路

POJ3613Cow Relays

Cow Relays POJ - 3613

poj3613:Cow Relays(倍增优化+矩阵乘法floyd+快速幂)

poj 3613 Cow Relays矩阵快速幂+Floyd

[POJ3613] Cow Relays(Floyd+矩阵快速幂)