最小环flody hdu6080

来源:互联网 发布:淘宝助理如何在线发货 编辑:程序博客网 时间:2024/05/16 13:47

度度熊保护村庄

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 118    Accepted Submission(s): 51


Problem Description
哗啦啦村袭击了喵哈哈村!

度度熊为了拯救喵哈哈村,带着自己的伙伴去救援喵哈哈村去了!度度熊与伙伴们很快的就过来占据了喵哈哈村的各个军事要地,牢牢的守住了喵哈哈村。

但是度度熊发现,这是一场旷日持久的战斗,所以度度熊决定要以逸待劳,保存尽量多的体力,去迎战哗啦啦村的战士。

于是度度熊决定派尽量多的人去休息,但是同时也不能松懈对喵哈哈村的保护。

换句话而言,度度熊希望尽量多的人休息,而且存在一个包围圈由剩下的人组成,且能够恰好的包围住喵哈哈村的所有住房(包括边界)。

请问最多能让多少个人休息呢?
 

Input
本题包含若干组测试数据。

第一行一个整数n,表示喵哈哈村的住房数量。

接下来n行,每行两个整数(x1[i],y1[i]),表示喵哈哈村的住房坐标。

第n+1行一个整数m,表示度度熊的士兵数量。

接下来m行,每行两个整数(x2[i],y2[i]),表示度度熊伙伴的坐标。

满足:

1<=n,m<=500

-10000<=x1[i],x2[i],y1[i],y2[i]<=10000
 

Output
请输出最多的人员休息的数目。

如果无法保护整个村庄的话,输出"ToT"
 

Sample Input
21 12 240 00 44 24 011 120 00 1
 

Sample Output
1

ToT

//后面的点集包围前面输入的点集#include <iostream>#include <cstdio>#include <cstdlib>#include <cstring>#include <algorithm>#include <cmath>#define N 1010#define EPS 1e-10#define INF 0x3f3f3f3fusing namespace std;struct PO {double x,y;} p[N],s[N];int n,m,dis[N][N];inline int dc(double x) {if(x>EPS) return 1;else if(x<-EPS) return -1;return 0;}inline double cross(const PO &a,const PO &b,const PO &c) {return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);}inline double dot(const PO &a,const PO &b,const PO &c) {return (b.x-a.x)*(c.x-a.x)+(b.y-a.y)*(c.y-a.y);}inline bool check(const PO &a,const PO &b) {for(int i=1; i<=n; i++)if(dc(cross(a,b,s[i]))<0) return false;return true;}inline void prep() {memset(dis,0x3f,sizeof dis);for(int i=1; i<=m; i++)for(int j=1; j<=m; j++) {if(i==j) continue;if(check(p[i],p[j])) dis[i][j]=1;}}inline bool onpoint() {for(int i=1; i<=m; i++) {int cnt=0;for(int j=1; j<=n; j++) {if(dc(p[i].x-s[j].x)==0&&dc(p[i].y-s[j].y)==0) cnt++;else break;}if(cnt==n) return true;}return false;}inline bool onseg() { //如果所有点共线,且没有一条线段包含他们所有点for(int i=3; i<=n; i++)if(dc(cross(s[1],s[2],s[i]))!=0) return false;for(int i=1; i<=m; i++)for(int j=i+1; j<=m; j++) {int cnt=0;for(int k=1; k<=n; k++)if(dc(dot(s[k],p[i],p[j]))<=0) cnt++;if(cnt==n) return false;}return true;}inline void go() {if(n==0) {printf("%d\n",m);return;}if(onpoint()) {printf("%d\n",m-1);return;}if(onseg()) {puts("ToT");return;}int ans=INF;prep();for(int k=1; k<=m; k++)for(int i=1; i<=m; i++)if(dis[i][k]<INF)for(int j=1; j<=m; j++)dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);for(int i=1; i<=m; i++) ans=min(ans,dis[i][i]);if(ans==INF) puts("ToT");elseprintf("%d\n",m-ans);}void read() {while(~scanf("%d",&n)) {double tmp;for(int i=1; i<=n; i++)scanf("%lf%lf",&s[i].x,&s[i].y);scanf("%d",&m);for(int i=1; i<=m; i++)scanf("%lf%lf",&p[i].x,&p[i].y);go();}}int main() {read();}