打篮球 ( pass ) 图论

【题目描述】

信息组的同学是比较喜欢打篮球的。

喜欢打篮球的同学很容易发现传球在篮球运动中是非常重要的，球传的好也就意味着球

打的顺。经常传球虽然不一定保证

100%

概率进球，但是总归是使球场气氛更加活跃了。

现在我们来看一场篮球比赛，肉蟹队vs宇宙大总统队。每支队伍都有N个(肉蟹)或(宇宙大总统)，现在轮到肉蟹队进攻。比赛可以看成在一个无限大的二维平面上进行，每个球员都有一个坐标篮筐也有一个坐标(50,0)。肉蟹们进攻的方式就是先传球若干次，再投篮。

一旦投篮，这一轮进攻就立即结束。

由于有宇宙大总统们的防御，所以肉蟹们传球和投篮都是有成功率的。

传球和投篮的轨迹可以看成一条线段。

关于传球和投篮的成功率如下：

传球：

Cp*(1-(ls/150)^2)*dr/(dr+1)

投篮：

(Cs*dr/(dr+1))^ln(ls)

其中

Cp和Cs

是给定的常数，ls是球的运动轨迹的长度，dr是防御队员的干扰系数，

由于宇宙大总统非常有智慧，所以大总统们知道他们既可以干扰传球的，也可以干扰接球的，既可以盖帽，也可以干扰球（这里干扰球不算犯规），因此dr就等于所有宇宙大总统到球运动轨迹两端点距离的最小值。

肉蟹们知道直接与宇宙大总统对抗是没有胜算的，所以她们对球场施了一个神秘的魔法，使得场上的所有人都瞬间变肉，以至于无法移动。

现在球在第一个肉蟹手上，请你设计一个传球和投篮的方案。求肉蟹们在这一轮就进球的概率最大是多少。

【输入数据】

第一行三个数，分别表示

N,Cp,Cs

接下来N行，每行两个数，表示肉蟹们的坐标

接下来N行，每行两个数，表示宇宙大总统们的坐标

【输出数据】

一行，一个数，即最大的概率。

精确到小数点后6位。

【输入样例1】

3 1 1

50 50

35 60

70 15

75 5

72 25

45 17

【输出样例1】

0.806732

【输入样例2】

1 0.5 0.5

50 4

50 3

【输出样例2】

0.146342

【输入样例3】

5 0.999999 0.8

50 50

40 50

40 40

40 30

50 20

50 41

44 29

48

27

45 41

48 64

【输出样例3】0.294171

【数据约定】

30%：1≤N≤ 5100%：1≤N≤ 500≤坐标≤ 1000≤Cp,Cs≤ 1

---

这题其实可以用ＤＰ和最短路做，我选择的是最短路。

我们只需预处理一些数据，然后就可以跑一遍ｓｐｆａ，然后在一一和篮筐判定选出答案就可以了。

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define ll long long
#define il inline
#define db double
using namespace std;
int n;
db cp,cs;
struct rou
{
    db x,y;
}r[145];
struct tong
{
    db x,y;
}t[145];
il db distance(db x1,db y1,db x2,db y2)
{
    db r1=x2-x1,r2=y2-y1;
    return (db)sqrt(r1*r1+r2*r2);
}
db dist1[145][145];
db dist2[145][145];
db dist[145];
db lan[145];
db gan[145];
db minx=9999;
bool vis[145];
il void spfa()
{
    int t[145];
    int head=0,tail=1;
    t[0]=1;
    dist[1]=1;
    while(head!=tail)
    {
        int h=t[head++];
        vis[h]=0;
        for(int i=1;i<=n;i++)
        {
            db ls=dist1[h][i],dr=min(gan[h],gan[i]);
            if(dist[i]<dist[h]*(cp*(1-(ls/150)*(ls/150))*(dr/(dr+1))))
            {
                dist[i]=dist[h]*(cp*(1-(ls/150)*(ls/150))*(dr/(dr+1)));
                if(!vis[i])
                {
                    vis[i]=1;
                    t[tail++]=i;
                }
            }
        }
    }
}
int main()
{
    freopen("pass.in","r",stdin);
    freopen("pass.out","w",stdout);
    cin>>n>>cp>>cs;
    for(int i=1;i<=n;i++)
    cin>>r[i].x>>r[i].y;
    for(int i=1;i<=n;i++)
    cin>>t[i].x>>t[i].y;
    for(int i=1;i<=n;i++)
    lan[i]=distance(r[i].x,r[i].y,50.0,0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            dist1[i][j]=dist1[j][i]=distance(r[i].x,r[i].y,r[j].x,r[j].y);
            dist2[i][j]=distance(r[i].x,r[i].y,t[j].x,t[j].y);
        }
    for(int i=1;i<=n;i++)
    gan[i]=9999;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            gan[i]=min(gan[i],dist2[i][j]);
    for(int i=1;i<=n;i++)
        minx=min(minx,distance(t[i].x,t[i].y,50.0,0));
    spfa();
    db ans=0;
    for(int i=1;i<=n;i++)
    {
        db dr=min(minx,gan[i]);
        ans=max(ans,dist[i]*(pow((cs*dr/(dr+1)),(log(lan[i])))));
    }
    printf("%f\n",ans);
    return 0;
}