L3-018 森森美图

解题思路：这道题用叉乘（这个我感觉比较有创意，我一开始想到的是用gcd来求分界的）将矩阵分成两部分，然后就是用spfa求最短路。

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
using namespace std;
int mp[106][105],v[105][105];
int dirx[]={1,0,-1,1,-1,1,0,-1};
int diry[]={1,1,1,0,0,-1,-1,-1};
int x11,y11,x2,y2,xs,ys,xe,ye;
double k=sqrt(2.0)-1;
int n,m;
int vis[106][106];
double dis[106][106];
int cross(int x,int y)
{
	x2=x-xs;
	y2=y-ys;
	int t=x11*y2-x2*y11;
	if(t==0) return 0;
	if(t>0) return 1;
	return -1;
}
struct node
{
	int x,y;
	node(){
	}
	node(int x,int y):x(x),y(y){
	}
};
queue<node> q;
double spfa(int f)
{
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
		{
			dis[i][j]=1e15;
			v[i][j]=0;
		}
	}
	dis[xs][ys]=mp[xs][ys];
	q.push(node(xs,ys));
	int x3,y3,x,y;
	double tp;
	while(!q.empty())
	{
		node t=q.front();
		q.pop();
		x=t.x;y=t.y;
		v[x][y]=0;
		for(int i=0;i<8;i++)
		{
			x3=x+dirx[i];
			y3=y+diry[i];
			if(x3<0||x3>=n||y3<0||y3>=m) continue;
			if(vis[x3][y3]!=f) continue;
			tp=dis[x][y]+mp[x3][y3];
			if(dirx[i]!=0&&diry[i]!=0)
			tp+=k*(mp[x3][y3]+mp[x][y]);
			if(dis[x3][y3]>tp)
			{
				dis[x3][y3]=tp;
				if(!v[x3][y3])
				{
					v[x3][y3]=1;
					q.push(node(x3,y3));
				}
			}
		} 
	}
	return dis[xe][ye];
}
int main()
{
	//freopen("t.txt","r",stdin);
	scanf("%d%d",&n,&m);
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
		{
			scanf("%d",&mp[i][j]);
		}
	}
	scanf("%d%d%d%d",&ys,&xs,&ye,&xe);
	x11=xe-xs,y11=ye-ys;
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
		{
			vis[i][j]=cross(i,j);
		}
	}
	double ans = 0 - mp[xs][ys] - mp[xe][ye];
	vis[xs][ys]=vis[xe][ye]=1;
	ans+=spfa(1);
	vis[xs][ys]=vis[xe][ye]=-1;
	ans+=spfa(-1);
	printf("%.2lf\n",ans);
	return 0;
}