L3-018 森森美图 - bfs

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的非负整数，即为像素点的分值。

最后一行给出用户选择的起始和结束像素点的坐标(X​start​​,Y​start​​)和(X​end​​,Y​end​​)。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

思路：

QwQ读题读好久嘤嘤嘤

有几个坑：

（1）横轴向右为正，纵轴向下为正，即题上输入的坐标是以横为x轴，竖为y轴

（2）判断条件在不在线的上方的方法：

• AB和AC的叉积大于0，则AC在AB上方

• AB和AC的叉积小于0，则AC在AB下方

（3）bfs的时候别忘了终止条件啊qwq，因为终点也是在线上的

代码入下:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>
using namespace std;
#define ll long long
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
typedef pair<int,int>P;
const double INF=0x3f3f3f3f;
const int N=115,mod=32767;
double mape[N][N];
int sx,sy,ex,ey,n,m;
int dx[8]={0,1,0,-1,-1,-1,1,1};
int dy[8]={-1,0,1,0,-1,1,1,-1};
double d[N][N];

map<P,int>vis;
int cross(P a,P b,P c){
    return ((b.first-a.first)*(c.second-a.second)-(b.second-a.second)*(c.first-a.first));
}
bool judge(int flag,int x,int y){
    if(x==ex&&y==ey)return true;
    if(x<0||x>=n||y<0||y>=m)return false;
    int tmp=cross(P(sx,sy),P(ex,ey),P(x,y));
    if(tmp>0&&flag||tmp<0&&!flag)return true;
    return false;
}

void bfs(int flag){
    vis.clear();
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++){
            d[i][j]=INF;
        }
    }
    queue<P>Q;
    Q.push(P(sx,sy));
    vis[P(sx,sy)]=1;
    d[sx][sy]=mape[sx][sy];
    while(Q.size()){
        P p=Q.front();
        vis[p]=0;
        Q.pop();
        for(int i=0;i<8;i++){
            int nx=p.first+dx[i];
            int ny=p.second+dy[i];
            if(!judge(flag,nx,ny))continue;
            double tmp=d[p.first][p.second]+mape[nx][ny];
            if(i>3)tmp+=(mape[nx][ny]+mape[p.first][p.second])*(sqrt(2)-1);
            if(d[nx][ny]>tmp){
                d[nx][ny]=tmp;
                if(!vis[P(nx,ny)]){
                    vis[P(nx,ny)]=1;
                    Q.push(P(nx,ny));
                }
            }
        }
    }
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++){
            scanf("%lf",&mape[i][j]);
        }
    }
    scanf("%d%d%d%d",&sy,&sx,&ey,&ex);
    bfs(1);
    double ans=d[ex][ey];
    bfs(0);
    ans+=d[ex][ey];
    printf("%.2f\n",ans-mape[sx][sy]-mape[ex][ey]+1e-6);
}