L3-018 森森美图 （30 分)
森森最近想让自己的朋友圈熠熠生辉，所以他决定自己写个美化照片的软件，并起名为森森美图。众所周知，在合照中美化自己的面部而不美化合照者的面部是让自己占据朋友圈高点的绝好方法，因此森森美图里当然得有这个功能。 这个功能的第一步是将自己的面部选中。森森首先计算出了一个图像中所有像素点与周围点的相似程度的分数，分数越低表示某个像素点越“像”一个轮廓边缘上的点。 森森认为，任意连续像素点的得分之和越低，表示它们组成的曲线和轮廓边缘的重合程度越高。为了选择出一个完整的面部，森森决定让用户选择面部上的两个像素点A和B，则连接这两个点的直线就将图像分为两部分，然后在这两部分中分别寻找一条从A到B且与轮廓重合程度最高的曲线，就可以拼出用户的面部了。 然而森森计算出来得分矩阵后，突然发现自己不知道怎么找到这两条曲线了，你能帮森森当上朋友圈的小王子吗？
为了解题方便，我们做出以下补充说明：
- 图像的左上角是坐标原点(0,0)，我们假设所有像素按矩阵格式排列，其坐标均为非负整数（即横轴向右为正，纵轴向下为正）。
- 忽略正好位于连接A和B的直线（注意不是线段）上的像素点，即不认为这部分像素点在任何一个划分部分上，因此曲线也不能经过这部分像素点。
- 曲线是八连通的（即任一像素点可与其周围的8个像素连通），但为了计算准确，某像素连接对角相邻的斜向像素时，得分额外增加两个像素分数和的√2倍减一。例如样例中，经过坐标为(3,1)和(4,2)的两个像素点的曲线，其得分应该是这两个像素点的分数和(2+2)，再加上额外的(2+2)乘以(√2−1)，即约为5.66。
输入格式：
输入在第一行给出两个正整数N和M（5≤N,M≤100），表示像素得分矩阵的行数和列数。
接下来N行，每行M个不大于1000的非负整数，即为像素点的分值。
最后一行给出用户选择的起始和结束像素点的坐标(Xstart,Ystart)和(Xend,Yend)。4个整数用空格分隔。
输出格式：
在一行中输出划分图片后找到的轮廓曲线的得分和，保留小数点后两位。注意起点和终点的得分不要重复计算。
输入样例：
6 6
9 0 1 9 9 9
9 9 1 2 2 9
9 9 2 0 2 9
9 9 1 1 2 9
9 9 3 3 1 1
9 9 9 9 9 9
2 1 5 4
输出样例：
27.04
#include<set>
#include<map>
#include<list>
#include<queue>
#include<stack>
#include<math.h>
#include<vector>
#include<bitset>
#include<iomanip>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#define eps (1e-8)
#define MAX 0x3f3f3f3f
#define u_max 1844674407370955161
#define l_max 9223372036854775807
#define i_max 2147483647
#define re register
#define pushup() tree[rt]=tree[rt<<1]+tree[rt<<1|1]
#define nth(k,n) nth_element(a,a+k,a+n); // 将 第K大的放在k位
#define ko() for(int i=2;i<=n;i++) s=(s+k)%i // 约瑟夫
#define ok() v.erase(unique(v.begin(),v.end()),v.end()) // 排序,离散化
#define Catalan C(2n,n)-C(2n,n-1) (1,2,5,14,42,132,429...) // 卡特兰数
using namespace std;
inline int read(){
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' & c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
typedef long long ll;
const double pi = atan(1.)*4.;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3fLL;
const int M=63;
const int N=1e5+5;
int n,m,sx,sy,ex,ey;
double s[105][105];
struct fun{
int x,y,px,py;
double ss;
}f[105][105];
struct cmp{ // 优先队列，用普通队列也可以
bool operator()(fun &a,fun &b)const{
return a.ss>b.ss;
}
};
// 判断点p在直线的位置（上方 ， 上边 ，下方）
int J(int x1,int y1,int x2,int y2,int x,int y){
return (x2-x1)*(y-y1)-(x-x1)*(y2-y1);
}
/*void print(int x,int y){ // 由于结果老是错，我打印了一下路径
if(x==-1&&y==-1) return ;
print(f[x][y].px,f[x][y].py);
printf("%d %d\n",f[x][y].x,f[x][y].y);
}*/
double mapp[105][105]; // 保存到达此点的最优结果
double bfs(int x,int y,int leap){ // 标记 leap为1，是只走直线上部分，为0，下部分
for(int i=0;i<n;i++) // 初始化
for(int j=0;j<m;j++)
mapp[i][j]=inf;
//v[x][y][=1;
// fun t;
f[x][y].x=x; f[x][y].y=y; f[x][y].ss=s[x][y];
f[x][y].px=-1,f[x][y].py=-1;
mapp[x][y]=s[x][y];
// t.x=x; t.y=y; t.ss=s[x][y];
priority_queue<fun,vector<fun>,cmp>qq;
//queue<fun>qq;
qq.push(f[x][y]);
while(!qq.empty()){
fun g=qq.top();
qq.pop();
if(g.x==ex&&g.y==ey){
//print(g.x,g.y);
return g.ss;
}
int xx,yy;
for(int i=0;i<8;i++){
int flag=0;
if(i==0){
xx=g.x+1;
yy=g.y;
}
else if(i==1){
xx=g.x-1;
yy=g.y;
}
else if(i==2){
xx=g.x;
yy=g.y+1;
}
else if(i==3){
xx=g.x;
yy=g.y-1;
}
else if(i==4){
xx=g.x-1;
yy=g.y-1;
flag=1;
}
else if(i==5){
xx=g.x-1;
yy=g.y+1;
flag=1;
}
else if(i==6){
xx=g.x+1;
yy=g.y-1;
flag=1;
}
else if(i==7){
xx=g.x+1;
yy=g.y+1;
flag=1;
}
if(xx<0||xx>=n||yy<0||yy>=m) continue;
int h=J(sx,sy,ex,ey,xx,yy);
if(h==0){
if(xx==ex&&yy==ey){
f[xx][yy].x=xx; f[xx][yy].y=yy;
f[xx][yy].px=g.x; f[xx][yy].py=g.y;
if(flag){
f[xx][yy].ss=g.ss+s[xx][yy];
f[xx][yy].ss+=(s[g.x][g.y]+s[xx][yy])*(sqrt(2.0)-1);
}
else
f[xx][yy].ss=g.ss+s[xx][yy];
//v[xx][yy]=1;
if(mapp[xx][yy]>f[xx][yy].ss){
mapp[xx][yy]=f[xx][yy].ss;
qq.push(f[xx][yy]);
}
}
}
else if(h>0){
if(leap){
f[xx][yy].x=xx; f[xx][yy].y=yy;
f[xx][yy].px=g.x; f[xx][yy].py=g.y;
if(flag){
f[xx][yy].ss=g.ss+s[xx][yy];
f[xx][yy].ss+=(s[g.x][g.y]+s[xx][yy])*(sqrt(2.0)-1);
}
else
f[xx][yy].ss=g.ss+s[xx][yy];
if(mapp[xx][yy]>f[xx][yy].ss){
mapp[xx][yy]=f[xx][yy].ss;
qq.push(f[xx][yy]);
}
}
}
else{
if(!leap){
f[xx][yy].x=xx; f[xx][yy].y=yy;
f[xx][yy].px=g.x; f[xx][yy].py=g.y;
if(flag){
f[xx][yy].ss=g.ss+s[xx][yy];
f[xx][yy].ss+=(s[g.x][g.y]+s[xx][yy])*(sqrt(2.0)-1);
}
else
f[xx][yy].ss=g.ss+s[xx][yy];
if(mapp[xx][yy]>f[xx][yy].ss){
mapp[xx][yy]=f[xx][yy].ss;
qq.push(f[xx][yy]);
}
}
}
}
}
}
int main(){
scanf("%d %d",&n,&m);
for(int i=0;i<n;i++)
for(int j=0;j<m;j++)
scanf("%lf",&s[i][j]);
scanf("%d %d %d %d",&sy,&sx,&ey,&ex);
double ans=bfs(sx,sy,1);
ans+=bfs(sx,sy,0);
printf("%.2lf\n",ans-s[sx][sy]-s[ex][ey]); // 起点与终点都加了两次
return 0;
}