BZOJ1527 : [POI2005]Pun-point
求出重心,然后把所有点关于重心极角排序,极角相同的按到重心距离从大到小排序。
按极角序依次扫描,得到相邻两个向量的夹角以及长度之比,看成字符串。
若两个字符串循环同构,则两个点集相似,KMP判断即可。
时间复杂度$O(n\log n)$。
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=25010;
const double eps=1e-10;
int T,l,i,j,k,nxt[N];
inline int sgn(double x){
if(x<-eps)return -1;
if(x>eps)return 1;
return 0;
}
struct P{
double x,y,o;
P(){}
P(double _x,double _y){x=_x,y=_y;}
P operator+(const P&b){return P(x+b.x,y+b.y);}
P operator-(const P&b){return P(x-b.x,y-b.y);}
P operator/(double b){return P(x/b,y/b);}
bool operator==(const P&b){return !sgn(x-b.x)&&!sgn(y-b.y);}
bool operator!=(const P&b){return sgn(x-b.x)||sgn(y-b.y);}
void cal(){o=atan2(y,x);}
double operator*(const P&b){return x*b.x+y*b.y;}
double len(){return sqrt(x*x+y*y);}
double lensqr()const{return x*x+y*y;}
}a[N],b[N],O;
inline bool cmp(const P&a,const P&b){
if(!sgn(a.o-b.o))return a.lensqr()>b.lensqr();
return a.o<b.o;
}
struct Set{
int n,m;P f[N];
void init0(){
n=l;
for(i=0;i<n;i++)a[i]=b[i];
}
void init1(){
n=l;
for(i=0;i<n;i++)a[i]=b[i],a[i].x*=-1;
}
void init(){
O=P(0,0);
for(i=0;i<n;i++)O=O+a[i];
O=O/((double)n);
for(i=j=m=0;i<n;i++)if(a[i]==O)m++;else a[j++]=a[i];
n=j;
for(i=0;i<n;i++)a[i]=a[i]-O,a[i].cal();
if(!n)return;
sort(a,a+n,cmp);
a[n]=a[0];
for(i=0;i<n;i++)f[i]=P(a[i]*a[i+1]/a[i].len()/a[i+1].len(),a[i].len()/a[i+1].len());
for(nxt[0]=j=-1,i=1;i<n;nxt[i++]=j){
while(~j&&f[j+1]!=f[i])j=nxt[j];
if(f[j+1]==f[i])j++;
}
}
}A,B;
void init(){
scanf("%d",&l);
for(i=0;i<l;i++)scanf("%lf%lf",&b[i].x,&b[i].y);
}
bool check(){
if(A.n!=B.n||A.m!=B.m)return 0;
if(!A.n)return 1;
for(j=-1,k=0;k<2;k++)for(i=0;i<A.n;i++){
while(~j&&B.f[j+1]!=A.f[i])j=nxt[j];
if(B.f[j+1]==A.f[i]){
j++;
if(j==A.n-1)return 1;
}
}
return 0;
}
int main(){
init();
A.init0();
A.init();
scanf("%d",&T);
while(T--){
init();
B.init0();
B.init();
if(check()){puts("TAK");continue;}
B.init1();
B.init();
puts(check()?"TAK":"NIE");
}
return 0;
}
BZOJ1527 : [POI2005]Pun-point的更多相关文章
- bzoj 1537: [POI2005]Aut- The Bus 线段树
bzoj 1537: [POI2005]Aut- The Bus 先把坐标离散化 设f[i][j]表示从(1,1)走到(i,j)的最优解 这样直接dp::: f[i][j] = max{f[i-1][ ...
- [BZOJ1529][POI2005]ska Piggy banks
[BZOJ1529][POI2005]ska Piggy banks 试题描述 Byteazar 有 N 个小猪存钱罐. 每个存钱罐只能用钥匙打开或者砸开. Byteazar 已经把每个存钱罐的钥匙放 ...
- 【BZOJ】【1532】【POI2005】Kos-Dicing
网络流/二分法 最大值最小……直接做不太好做的时候就可以用二分+判定来搞. 这题我们就也可以二分最大胜场v,那么怎么来判定呢?首先我们发现:每场比赛要么A赢,要么B赢,这一点跟二分图匹配非常类似,那么 ...
- Bzoj 1532: [POI2005]Kos-Dicing 二分,网络流
1532: [POI2005]Kos-Dicing Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1373 Solved: 444[Submit][St ...
- BZOJ1533: [POI2005]Lot-A Journey to Mars
1533: [POI2005]Lot-A Journey to Mars Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 174 Solved: 76[S ...
- BZOJ1528: [POI2005]sam-Toy Cars
1528: [POI2005]sam-Toy Cars Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 282 Solved: 129[Submit][S ...
- BZOJ1537: [POI2005]Aut- The Bus
1537: [POI2005]Aut- The Bus Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 158 Solved: 100[Submit][S ...
- BZOJ1532: [POI2005]Kos-Dicing
1532: [POI2005]Kos-Dicing Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1060 Solved: 321[Submit][St ...
- BZOJ 1537: [POI2005]Aut- The Bus(dp + BIT)
对y坐标离散化, 然后按x坐标排序, dp. 一个点(x, y), 设到达这个点接到的最多乘客数为t, 那么t可以用来更新y'>=y的所有点.用树状数组维护最大值. -------------- ...
随机推荐
- centos 7配置网络 更新yum源
cd /etc/sysconfig/network-script/ 找到对应的ifcfg-entxxxx文件,然后添加网关,修改dhcp为static,静态ip,添加IPADDR ip地址.onboo ...
- socket.io简单入门(一.实现简单的图表推送)
引子:随着nodejs蓬勃发展,虽然主要业务系统因为架构健壮性不会选择nodejs座位应用服务器.但是大量的内部系统却可以使用nodejs试水,大量的前端开发人员转入全堆开发也是一个因素. 研究本例主 ...
- 高程三:Array
一:Array数组 1.Array.isArray(参数) 检测是否是数组,*不兼容IE8,兼容IE9及以上.Chrome.Firefox等,要兼容IE8,可以用 Object.prototype.t ...
- PHP cURL 使用cookie 模拟登录
cURL是什么 cURL: http://php.net/manual/zh/book.curl.php PHP 支持 Daniel Stenberg 创建的 libcurl 库,能够连接通讯各种服务 ...
- C#在二维码中添加圆角logo
public class QRCodeHelper { #region 合并用户QR图片和用户头像 /// <summary> /// 合并用户QR图片和用户头像 /// </sum ...
- Nodejs 饭店
一个Node.js饭店的发展历程 前面的一堆理论似乎不太好明白,最后讲一个关于饭店发展历程的故事作为结尾吧. 第一年 饭店开张,只有一个厨师(同时还兼任老板.服务员.打荷.收银员),当一个客 ...
- Python全栈【异常处理】
Python全栈[异常处理] 本节内容: 1.异常处理 2.什么时候用异常处理 异常处理 1.异常处理: 异常就是程序运行时发生错误的信号,异常处理是在编程过程中为了增加友好性,在程序出现bug时一般 ...
- appium 自动化测试之知乎Android客户端
appium是一个开源框架,相对来说还不算很稳定.转载请注明出处!!!! 前些日子,配置好了appium测试环境,至于环境怎么搭建,参考:http://www.cnblogs.com/tobecraz ...
- 二维码跳转不同的 app store
说道二维码 之前是用来跳转app store 然后在就是出来的 扫码付款什么的 用的很平常,其实里面也很简单 自己刚开始接触的时候 同事说要做一个二维码下载 应用 => 我=懵逼 ...
- Tomcat7.0+ web.xml问题
Tomcat7+版本的web.xml都加上 <context-param> <param-name>webAppRootKey</param-name> <p ...