//一些点一些圆,过圆不能连线,相邻点不能连线,问最多连几条线
//计算几何模板+区间dp
//关键是判断圆和线段是否相交
#include <cstdio>
#include <cmath>
#include <algorithm>
#define N 500
#define ll long long
#define sqr(x) ((x)*(x))
using namespace std;
bool ok[N][N];
ll n,m,r,f[N][N];
struct poi
{
ll x,y;
poi(ll X,ll Y)
{
x=X;y=Y;
}
poi()
{
}
} p[N],q[N];
ll cro(poi a,poi b)
{
return (a.x*b.y-a.y*b.x);
}
ll dot(poi a,poi b)
{
return a.x*b.x+a.y*b.y;
}
poi tr(poi a,poi b)
{
return poi(a.x-b.x,a.y-b.y);
}
double dis(poi a,poi b)
{
return sqrt((double)(sqr(a.x-b.x)+sqr(a.y-b.y)));
}
bool cmp(poi a,poi b)
{
return (a.x<b.x)||(a.x==b.x && a.y<b.y);
}
bool cmp2(poi a,poi b)
{
return cro(tr(a,p[]),tr(b,p[]))>;
}
bool check(poi a,poi b)
{
double C=dis(a,b);
for(int i=;i<=m;i++)
{
double A=dis(q[i],a);
double B=dis(q[i],b);
double COSA=(B*B+C*C-A*A)/(*B*C);
double COSB=(A*A+C*C-B*B)/(*A*C);
if(COSA<1e-||COSB<1e-){
if(A<1.0*r+1e-||B<1.0*r+1e-)return ;
else continue;
}
double SINA=sqrt(-COSA*COSA);
double H=B*SINA;
if(H<1.0*r+1e-)return ;
}
return ;
}
int main()
{
scanf("%d%d%d",&n,&m,&r);
for(int i=;i<=n;i++)
scanf("%d%d",&p[i].x,&p[i].y);
sort(p+,p++n,cmp);
sort(p+,p++n,cmp2);
for(int i=;i<=m;i++)
scanf("%d%d",&q[i].x,&q[i].y);
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
if(i!= || j!=n)
ok[i][j]=check(p[i],p[j]);
for(int len=;len<=n;len++)
for(int start=;start<=n;start++){
ll end=start+len;
for(int k=start+;k<min(end,n);k++)
f[start][end]=max(f[start][end],f[start][k]+f[k][end]);
f[start][end]+=ok[start][end];
}
printf("%d\n",f[][n]);
}

POJ3178 计算几何+DP的更多相关文章

  1. sdut 2153:Clockwise(第一届山东省省赛原题,计算几何+DP)

    Clockwise Time Limit: 1000ms   Memory limit: 65536K  有疑问?点这里^_^ 题目描述 Saya have a long necklace with ...

  2. poj3178 Roping the Field (计算几何 + dp)

    Roping the Field Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 858   Accepted: 250 De ...

  3. 『HGOI 20190917』Cruise 题解 (计算几何+DP)

    题目概述 在平面直角坐标系的第$1$象限和第$4$象限有$n$个点,其中第$i$个点的坐标为$(x_i,y_i)$,有一个权值$p_i$ 从原点$O(0,0)$出发,不重复的经过一些点,最终走到原点, ...

  4. HDU 4562 守护雅典娜 (计算几何+DP)

    守护雅典娜 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total Submi ...

  5. Code Chef MINPOLY(计算几何+dp)

    题面 传送门 题解 我们枚举这个凸多边形\(y\)坐标最小的点\(p_i\),然后对于所有\(y\)坐标大于等于它的点极角排序 我们预处理出\(s_{j,k}\)表示三角形\(p_i,p_j,p_k\ ...

  6. bzoj 3778: 共鸣【计算几何+dp】

    枚举起点,然后设f[i][j]为上凸壳上一个点是i当前点是j的最大面积,g是下凸壳,然后合并的时候枚举结束点t合并上下凸壳即可 这样的好处是每次转移都是往凸多边形里加一个三角形(s,i,j),所以判断 ...

  7. SDOI2019 省选前模板整理

    目录 计算几何✔ DP 斜率优化✔ 四边形不等式✔ 轮廓线DP✘ 各种分治 CDQ分治✔ 点分治✔ 整体二分✔ 数据结构 线段树合并✔ 分块✔ K-D Tree LCT 可持久化Trie✔ Splay ...

  8. 2019.03.15 ZJOI2019模拟赛 解题报告

    得分: \(20+45+15=80\)(三题暴力全写挂...) \(T1\):Lyk Love painting 首先,不难想到二分答案然后\(DP\)验证. 设当前需验证的答案为\(x\),则一个暴 ...

  9. HZNU Training 4 for Zhejiang Provincial Collegiate Programming Contest 2019

    今日这场比赛我们准备的题比较全面,二分+数论+最短路+计算几何+dp+思维+签到题等.有较难的防AK题,也有简单的签到题.为大家准备了一份题解和AC代码. A - Meeting with Alien ...

随机推荐

  1. matlab之sortrows()函数

    sortrows()函数的格式: sortrows(A,column) A是一个矩阵,如果没有第二个参数column,则默认按照第一列升序排列,如果遇到重复数字,则按照第二列升序排列,依次类推... ...

  2. CSS样式命名整理

    CSS样式命名整理 页面结构 容器: container/wrap 整体宽度:wrapper 页头:header 内容:content 页面主体:main 页尾:footer 导航:nav 侧栏:si ...

  3. github添加ssh公钥

    使用git作为版本维护工具非常方便,而且一般个人用github作为远端库就够用了.而一般git连接github的方式采用ssh的方法,http的会略微慢一些.所以为了方便一般会在github设置中添加 ...

  4. [原创]Java在线编辑word文档调用PageOffice实现并发控制

    1.功能介绍 PageOffice的并发控制功能用来解决多个用户在线编辑同一篇文档可能造成的互相覆盖修改结果的技术难题. B/S架构下用户访问都是并发的,也就是说经常会出现同时N个用户对一个服务器页面 ...

  5. hdu-5804 Price List(水题)

    题目链接: Price List Time Limit: 2000/1000 MS (Java/Others)     Memory Limit: 262144/131072 K (Java/Othe ...

  6. VijosP1180:选课

    描述 学校实行学分制.每门的必修课都有固定的学分,同时还必须获得相应的选修课程学分.学校开设了N(N<300)门的选修课程,每个学生可选课程的数量M是给定的.学生选修了这M门课并考核通过就能获得 ...

  7. xml解析中的sax解析

    title: xml解析中的sax解析 tags: grammar_cjkRuby: true --- SAXPasser 类: parser(File file, DefaultHandler ha ...

  8. Python3中内置类型bytes和str用法及byte和string之间各种编码转换,python--列表,元组,字符串互相转换

    Python3中内置类型bytes和str用法及byte和string之间各种编码转换 python--列表,元组,字符串互相转换 列表,元组和字符串python中有三个内建函数:,他们之间的互相转换 ...

  9. 小程序[邮箱提取器-EmailSplider]总结

    1.背景情况     学东西做快的是付诸实践,写这个小程序的目的就是为了综合运用各个知识点,从而提升学习的效果.   2.涉及知识     A.Swing 的布局     B.Swing中,线程访问U ...

  10. 2018上海大都会 J-Beautiful Numbers(数位dp-模未知)

    J-Beautiful Numbers 链接:https://www.nowcoder.com/acm/contest/163/J 来源:牛客网 时间限制:C/C++ 8秒,其他语言16秒 空间限制: ...