【noip2001 提高组T4】 Car的旅行路线 预处理+最短路

又到暑假了，住在城市A的Car想和朋友一起去城市B旅游。她知道每个城市都有四个飞机场，分别位于一个矩形的四个顶点上，同一个城市中两个机场之间有一条笔直的高速铁路，第I个城市中高速铁路了的单位里程价格为Ti，任意两个不同城市的机场之间均有航线，所有航线单位里程的价格均为t。
那么Car应如何安排到城市B的路线才能尽可能的节省花费呢?她发现这并不是一个简单的问题，于是她来向你请教。找出一条从城市A到B的旅游路线，出发和到达城市中的机场可以任意选取，要求总的花费最少。

输入格式
第一行有四个正整数s，t，A，B。S(0 < S<=100)表示城市的个数，t表示飞机单位里程的价格，A，B分别为城市A，B的序号，(1<=A，B<=S)。
接下来有S行，其中第I行均有7个正整数xi1，yi1，xi2，yi2，xi3，yi3，Ti，这当中的(xi1，yi1)，(xi2，yi2)，(xi3，yi3)分别是第I个城市中任意三个机场的坐标，T I为第I个城市高速铁路单位里程的价格。
输出格式
输出最小费用(结果保留两位小数)
样例输入1
3 10 1 3
1 1 1 3 3 1 30
2 5 7 4 5 2 1
8 6 8 8 11 6 3
样例输出1
47.55

题解：题目极其恶心，处理第四个点和建边非常麻烦。需要认真+仔细。预处理维护出每两个点之间的距离，然后跑一遍floyed就ok了。其实这就是一道最短路裸题。

#include<iostream>
#include<cstdio>
#include<cmath>
typedef double DD;
using namespace std;
struct node
{
    int x[5],y[5];
}lg[106];
int t,cost[106],s,A,B,n;
DD f[106][106],ans=1000000.0;
int abs(int x)
{
    if (x>=0) return x;
    else return -x;
}
DD length(int i,int u,int v)
{
    return sqrt(abs(lg[i].x[u]-lg[i].x[v])*abs(lg[i].x[u]-lg[i].x[v])*1.0+abs(lg[i].y[u]-lg[i].y[v])*abs(lg[i].y[u]-lg[i].y[v])*1.0);
}
void read(int i)    //处理第四个点，利用三角形三边长度关系处理 
{
    DD len1,len2,len3,maxx,midx,midy;
    len1=length(i,1,2);//1,2
    len2=length(i,1,3);//1,3
    len3=length(i,2,3);//2,3
    maxx=max(len1,max(len2,len3));
    if(maxx==len1)
    {
        midx=(lg[i].x[1]+lg[i].x[2])/2.0;
        midy=(lg[i].y[1]+lg[i].y[2])/2.0;
        lg[i].x[4]=midx+midx-lg[i].x[3]*1.0;
        lg[i].y[4]=midy+midy-lg[i].y[3]*1.0;
    }
    if(maxx==len2)
    {
        midx=(lg[i].x[1]+lg[i].x[3])/2.0;
        midy=(lg[i].y[1]+lg[i].y[3])/2.0;
        lg[i].x[4]=midx+midx-lg[i].x[2]*1.0;
        lg[i].y[4]=midy+midy-lg[i].y[2]*1.0;
    }
    if(maxx==len3)
    {
        midx=(lg[i].x[2]+lg[i].x[3])/2.0;
        midy=(lg[i].y[2]+lg[i].y[3])/2.0;
        lg[i].x[4]=midx+midx-lg[i].x[1]*1.0;
        lg[i].y[4]=midy+midy-lg[i].y[1]*1.0;
    }

}
DD air(int i,int j,int u,int v)    
{
    return sqrt(abs(lg[i].x[j]-lg[u].x[v])*abs(lg[i].x[j]-lg[u].x[v])*1.0+abs(lg[i].y[j]-lg[u].y[v])*abs(lg[i].y[j]-lg[u].y[v])*1.0)*t*1.0;
}
DD street(int i,int u,int v)
{
    return sqrt(abs(lg[i].x[u]-lg[i].x[v])*abs(lg[i].x[u]-lg[i].x[v])*1.0+abs(lg[i].y[u]-lg[i].y[v])*abs(lg[i].y[u]-lg[i].y[v])*1.0)*cost[i]*1.0;
}
int main()
{
    int i,j,i1,i2,j1,j2,x,y,k;
    scanf("%d",&n);
    while (n>0)
    {
        n--;
    ans=10000000000;
    scanf("%d%d%d%d",&s,&t,&A,&B);
    for (i=1;i<=s;i++)
    {
        scanf("%d%d",&x,&y);
        lg[i].x[1]=x;lg[i].y[1]=y;
        scanf("%d%d",&x,&y);
        lg[i].x[2]=x;lg[i].y[2]=y;
        scanf("%d%d",&x,&y);
        lg[i].x[3]=x;lg[i].y[3]=y;
        scanf("%d",&cost[i]);
        read(i);
    }
    for (i1=1;i1<=s;i1++)  //处理飞机线路 
        for (j1=1;j1<=4;j1++)
            for (i2=1;i2<=s;i2++)
                //if (i!=j)
                    for (j2=1;j2<=4;j2++)
                    {
                        f[(i1-1)*4+j1][(i2-1)*4+j2]=f[(i2-1)*4+j2][(i1-1)*4+j1]=air(i1,j1,i2,j2);
                    }
    for (i=1;i<=s;i++)    //处理高速路 
        for (j1=1;j1<=4;j1++)
            for (j2=1;j2<=4;j2++)
                f[(i-1)*4+j1][(i-1)*4+j2]=f[(i-1)*4+j2][(i-1)*4+j1]=street(i,j1,j2);
    for (k=1;k<=s*4;k++)
        for (i=1;i<=s*4;i++)
            for (j=1;j<=s*4;j++)
//              if ((i!=j)&&(j!=k)&&(i!=k))
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
//  for (i=1;i<=s;i++)
//  {
//      for (j=1;j<=4;j++)
//          printf("%d %d  ",lg[i].x[j],lg[i].y[j]);
//      printf("\n");
//  }
//      
//  for (i=1;i<=s*4;i++)
//  {
//      for (j=1;j<=s*4;j++)
//      {
//          printf("%lf   ",f[i][j]);
//      }
//      printf("\n");
//  }
    for (i=(A-1)*4+1;i<=(A-1)*4+4;i++)
        for (j=(B-1)*4+1;j<=(B-1)*4+4;j++)
            if (f[i][j]<ans) 
            ans=f[i][j];
    printf("%.1lf",ans);
}
}