BZOJ 3630 JLOI2014 镜面通道 计算几何+最小点割集

题目大意：给定一个二维的通道，通道内有一些正方形和圆形的零件(允许重叠)，求最少删掉多少零件可以使光线通过反射通过这个通道

看到这题当时我们都吓得说日语了有木有啊！これはいったい何ですか？！計算幾何ですか？圖論ですか？やれますか？全然知らないあああああ！

我在第一第二题上耽误了太长的时间 导致第三题没时间了(第二题还没看到多组数据爆零了 令人感动不已) 最后读入全体数据后加和再加上我的生日然后对n取模+1，零分。。

全场唯一知道正解的写挂了 @Vmurder 写了个枚举射入的坐标横向射入，遇到元件就干掉，拿了50分

最可怕的是有个同样思路，只不过是斜射的怒切了 世界你逗我

这题正解是最小点割集

首先这里有个物理上的费马小定理 具体内容我不知道 但是总之放在这题上就是“水能过去光就能过去”

于是我们还考虑反射什么乱七八糟的搞毛！直接处理圆形和方形，相交就连边，上是源下是汇，求个最小点割集，这题就切了！

确实是道不错的题 只是我的代码写的有点渣 我慢しろ！

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#define c(a) memset(a,0,sizeof a)
#define M 610
#define INF (0x7fffffff)
struct abcde{int p,x,y,r,x1,y1,x2,y2;}a[M>>1];
struct abcd{int to,f,next;}table[200000];int head[M],tot=1;
inline void adddd(int x,int y,int z){table[++tot].to=y;table[tot].f=z;table[tot].next=head[x];head[x]=tot;}
inline void addd(int x,int y,int z){adddd(x,y,z),adddd(y,x,0);}
inline void add(int x,int y){addd(x<<1|1,y<<1,INF),addd(y<<1|1,x<<1,INF);}
inline void swap(int&x,int&y){int z=x;x=y;y=z;}
inline double dis(int x1,int y1,int x2,int y2){return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2))-(1e-10);}
inline int min(int x,int y){return x<y?x:y;}
int cx,cy,n,ans,s,t;
bool Judge(int x,int y)
{
    int p=a[x].p+a[y].p;
    if(p==2)return dis(a[x].x,a[x].y,a[y].x,a[y].y)<=a[x].r+a[y].r?1:0;
    if(p==3)
    {
        if(a[x].p==1)swap(x,y);
        if(dis(a[x].x1,a[x].y1,a[y].x,a[y].y)<=a[y].r)return 1;
        if(dis(a[x].x1,a[x].y2,a[y].x,a[y].y)<=a[y].r)return 1;
        if(dis(a[x].x2,a[x].y1,a[y].x,a[y].y)<=a[y].r)return 1;
        if(dis(a[x].x2,a[x].y2,a[y].x,a[y].y)<=a[y].r)return 1;
        if(a[x].x2>a[y].x&&a[x].x1<a[y].x)if(abs(a[x].y2-a[y].y)<=a[y].r||abs(a[x].y1-a[y].y)<=a[y].r)return 1;
        if(a[x].y2>a[y].y&&a[x].y1<a[y].y)if(abs(a[x].x2-a[y].x)<=a[y].r||abs(a[x].x1-a[y].x)<=a[y].r)return 1;
        if(a[x].x1<=a[y].x&&a[x].x2>=a[y].x&&a[x].y1<=a[y].y&&a[x].y2>=a[y].y)return 1;
        return 0;
    }
    if(p==4)
    {
        if((a[x].x1>=a[y].x1&&a[x].x1<=a[y].x2)||(a[x].x2>=a[y].x1&&a[x].x2<=a[y].x2))
        if((a[x].y1>=a[y].y1&&a[x].y1<=a[y].y2)||(a[x].y2>=a[y].y1&&a[x].y2<=a[y].y2))return 1;
        if(a[x].x1<=a[y].x1&&a[y].x2<=a[x].x2&&a[y].y1<=a[x].y1&&a[x].y2<=a[y].y2)return 1;
        swap(x,y);
        if((a[x].x1>=a[y].x1&&a[x].x1<=a[y].x2)||(a[x].x2>=a[y].x1&&a[x].x2<=a[y].x2))
        if((a[x].y1>=a[y].y1&&a[x].y1<=a[y].y2)||(a[x].y2>=a[y].y1&&a[x].y2<=a[y].y2))return 1;
        if(a[x].x1<=a[y].x1&&a[y].x2<=a[x].x2&&a[y].y1<=a[x].y1&&a[x].y2<=a[y].y2)return 1;
        return 0;
    }
}
int q[65540],f[M],d[M];
unsigned short r,h;
bool bfs()
{
    int i;
    memset(d,0,sizeof d);
    h=r;q[++r]=s;d[s]=1;
    while(h!=r)
    {
        h++;
        for(i=head[q[h]];i;i=table[i].next)
        if(table[i].f&&!d[table[i].to])
        {
            q[++r]=table[i].to;
            d[table[i].to]=d[q[h]]+1;
            if(table[i].to==t)return 1;
        }
    }
    return 0;
}
int dinic(int x,int flow)
{
    int temp=flow,k,i;
    if(x==t)return flow;
    for(i=head[x];i;i=table[i].next)
    if(table[i].f&&temp&&d[table[i].to]==d[x]+1)
    {
        k=dinic(table[i].to,min(temp,table[i].f) );
        if(!k)d[table[i].to]=0;
        table[i].f-=k;
        table[i^1].f+=k;
        temp-=k;
    }
    return flow-temp;
}
int main()
{
    int i,j,flow;
    scanf("%d%d%d",&cx,&cy,&n);
    s=1;t=n+1<<1;
    a[0].p=2;a[0].x1=0;a[0].y1=cy;a[0].x2=cx;a[0].y2=cy;
    a[n+1].p=2;a[n+1].x1=0;a[n+1].y1=0;a[n+1].x2=cx;a[n+1].y2=0;
    for(i=1;i<=n;i++)
    {
        scanf("%d",&a[i].p);
        if(a[i].p==1)scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].r);
        else scanf("%d%d%d%d",&a[i].x1,&a[i].y1,&a[i].x2,&a[i].y2);
        for(j=0;j<i;j++)if(Judge(i,j))add(i,j);
    }
    for(j=1;j<i;j++)if(Judge(j,i))add(j,i);
    for(i=1;i<=n;i++)adddd(i<<1,i<<1|1,1),adddd(i<<1|1,i<<1,1);
    while( bfs() )
        while( flow=dinic(s,INF) )
            ans+=flow;
    printf("%d",ans);
}