14992: [Usaco2017 Feb]Why Did the Cow Cross the Road

#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long LL;
const int N=100+5,M=N*N*16;
const int dx[16]={-3, 0, 3, 0,-2,-2, 2, 2,-1,-1, 1, 1, 1,-1, 0, 0};
const int dy[16]={ 0,-3, 0, 3, 1,-1, 1,-1, 2,-2, 2,-2, 0, 0, 1,-1};
int n;
LL a[N][N],t;
struct Edge{
    int cnt,y[M],nxt[M],fst[N*N];
    LL z[M];
    void set(){
        cnt=0;
        memset(y,0,sizeof y);
        memset(z,0,sizeof z);
        memset(nxt,0,sizeof nxt);
        memset(fst,0,sizeof fst);
    }
    void add(int a,int b,LL c){
        cnt++;
        y[cnt]=b,z[cnt]=c;
        nxt[cnt]=fst[a],fst[a]=cnt;
    }
}e;
bool check(int x,int y){
    return 1<=x&&x<=n&&1<=y&&y<=n;
}
int hash(int x,int y){
    return (x-1)*n+y-1;
}
void adden(int x,int y){
    if (!check(x,y))
        return;
    e.add(hash(x,y),hash(n,n),t*(abs(n-x)+abs(n-y)));
}
LL dis[N*N];
bool f[N*N];
int q[N*N],qmod,head,tail;
void spfa(){
    for (int i=0;i<n*n;i++)
        dis[i]=1LL<<55;
    memset(f,0,sizeof f);
    qmod=n*n+1;
    head=0,tail=0;
    q[tail=(tail+1)%qmod]=0;
    dis[0]=0;
    f[0]=1;
    while (head!=tail){
        int x=q[head=(head+1)%qmod];
        f[x]=0;
        for (int i=e.fst[x];i;i=e.nxt[i]){
            int y=e.y[i];
            LL z=e.z[i];
            if (dis[y]>dis[x]+z){
                dis[y]=dis[x]+z;
                if (!f[y]){
                    f[y]=1;
                    q[tail=(tail+1)%qmod]=y;
                }
            }
        }
    }
}
int main(){
    scanf("%d%lld",&n,&t);
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            scanf("%lld",&a[i][j]);
    e.set();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++){
            int x=hash(i,j);
            for (int k=0;k<16;k++){
                int i_=i+dx[k],j_=j+dy[k];
                if (check(i_,j_))
                    e.add(x,hash(i_,j_),t*3LL+a[i_][j_]);
            }
        }
    adden(n-2,n),adden(n-1,n-1),adden(n,n-2),adden(n-1,n),adden(n,n-1);
    spfa();
    printf("%lld",dis[hash(n,n)]);
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值