【POI2007】OSI-Axes of Symmetry【计算几何】【manacher】

题意:给一个 nnn 个点的多边形,求对称轴个数。

n≤105n\leq 10^5n105

显然对称轴一定在顶点或边的中点上。

但你 n2n^2n2 枚举完全没有一点能过的样子。

冷静分析,发现有 “中点”,“对称轴”,很自然个鬼地想到了manacher。

在边的中点插入一个点,然后复制一遍断环成链。 然后跑马拉车,扩展的时候判断是否轴对称。

设点 iii 可以扩展到 [i−pi,i+pi][i-p_i,i+p_i][ipi,i+pi],如果扩展到整个多边形就是合法的对称轴,即 2pi+1≥2n2p_i+1\geq 2n2pi+12npi≥np_i\geq npin,并且只有第一圈的点会有贡献。一个对称轴会算两次,除以 222 就是答案。

方便实现的小trick:

  1. 边界的地方随机一个点就不用特判。
  2. 判轴对称可以算出中点,用叉积判在不在已知的对称轴上。如果对称轴未确定就设成 000 向量。但要注意本来就确定的时候要特判一下不要把它改回零向量。
  3. 读入的坐标都乘上 444,就可以不用 double。

复杂度 O(n)O(n)O(n)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <cstdlib>
#define MAXN 400005
using namespace std;
inline int read()
{int ans=0,f=1;char c=getchar();while (!isdigit(c)) (c=='-')&&(f=-1),c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return f*ans;
}
typedef long long ll;
int x[MAXN],y[MAXN];
int p[MAXN],maxr,mid,cnt;
int lasx,lasy;
inline bool check(int l,int i,int r)
{if ((ll)(x[l]-x[i])*(x[l]-x[i])+(ll)(y[l]-y[i])*(y[l]-y[i])!=(ll)(x[r]-x[i])*(x[r]-x[i])+(ll)(y[r]-y[i])*(y[r]-y[i]))return false;int tx=(x[l]+x[r])/2-x[i],ty=(y[l]+y[r])/2-y[i];if ((ll)tx*lasy!=(ll)ty*lasx) return false;if ((ll)tx*tx+(ll)ty*ty) lasx=tx,lasy=ty;return true;
}
int main()
{for (int T=read();T;T--){int n=read();for (int i=1;i<=2*n;i+=2) x[i]=read()*4,y[i]=read()*4;for (int i=2*n+1;i<=4*n;i+=2) x[i]=x[i-2*n],y[i]=y[i-2*n];x[4*n+1]=x[1],y[4*n+1]=y[1];for (int i=2;i<=4*n;i+=2) x[i]=(x[i-1]+x[i+1])/2,y[i]=(y[i-1]+y[i+1])/2;x[0]=rand(),y[0]=rand(),x[4*n+1]=rand(),y[4*n+1]=rand();maxr=mid=cnt=0;for (int i=1;i<=4*n;i++){if (i<maxr) p[i]=min(p[2*mid-i],maxr-i);else p[i]=0;lasx=lasy=0;while (check(i-p[i]-1,i,i+p[i]+1)) ++p[i];if (i+p[i]>maxr) maxr=i+p[i],mid=i;cnt+=(p[i]>=n);}printf("%d\n",cnt/2);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314925.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

端到端测试实践:Jenkins集成TestCafe

上一篇《对产品质量的一点思考》中说到自动化测试的重要性&#xff0c;本文简单介绍下怎样在实际项目中实现端到端测试的自动化&#xff0c;在这里我们使用的端到端测试工具是TestCafe。环境Jenkisn&#xff1a;2.183TestCafe&#xff1a;1.3.0为什么采用TestCafe做自动化测试要…

Codeforces Round #727 (Div. 2) F. Strange Array 线段树 + 区间合并 + 排序优化

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的数组&#xff0c;对每个位置iii求一个最大价值&#xff0c;价值计算方式如下&#xff1a;选择一个包含iii的[l,r][l,r][l,r]&#xff0c;让后将其拿出来排序&#xff0c;之后价值就是当前…

【CF487E】Tourists【圆方树】【树链剖分】【multiset】

题意&#xff1a;给一张 nnn 点 mmm 边的连通无向图&#xff0c;点帯权&#xff0c;qqq 次操作&#xff1a; 修改一个点的权值。询问两点间所有简单路的最小权值的最小值。 n,m,q≤105n,m,q\leq 10^5n,m,q≤105 显然建出圆方树然后询问路径最小值。多半要树链剖分了。 对于方…

再见Jenkins,从Gitlab代码提交到k8s服务持续交付只需七毛三

Gitlab runner 快速搭建CICD pipeline背景日常开发中&#xff0c;相信大家已经做了很多的自动化运维环境&#xff0c;用的最多的想必就是利用Jenkins实现代码提交到自动化测试再到自动化打包&#xff0c;部署全流水线Jenkins在devops担任了很重要的角色&#xff0c;但是另一方面…

Codeforces Round #726 (Div. 2) F. Figure Fixing 二分图 + 思维

传送门 文章目录题意&#xff1a;思路题意&#xff1a; 给你一张nnn个点mmm条边的图&#xff0c;每个点都有一个当前值aia_iai​&#xff0c;目标值bib_ibi​&#xff0c;每次可以选择一条边(i,j)(i,j)(i,j)&#xff0c;将ai,aja_i,a_jai​,aj​都加上任意一个数&#xff0c;问…

【NOI2020】命运【树形dp】【线段树合并】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;每条边需要染成黑白两种颜色中的一种。给出 mmm 个条件&#xff0c;每个条件给出 u,vu,vu,v&#xff0c;其中 uuu 是 vvv 的祖先&#xff0c;要求 uuu 到 vvv 的链上至少一条黑边。求方案数 模 998244353998244353998244353。 …

.NET Core 3.0之深入源码理解Kestrel的集成与应用(二)

前言前一篇文章主要介绍了.NET Core继承Kestrel的目的、运行方式以及相关的使用&#xff0c;接下来将进一步从源码角度探讨.NET Core 3.0中关于Kestrel的其他内容&#xff0c;该部分内容&#xff0c;我们无需掌握&#xff0c;依然可以用好Kestrel&#xff0c;本文只是将一些内部…

Codeforces Round #726 (Div. 2) E2. Erase and Extend (Hard Version) 贪心

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的串sss&#xff0c;你有两个操作可以使用&#xff1a; (1)(1)(1)从sss的结尾删除一个字母。 (2)sss(2)sss(2)sss。 让你通过若干次操作使其变成一个长度为kkk的串&#xff0c;且其字典序最…

【BZOJ3451】Normal【期望线性性】【点分治】【NTT卷积】

题意&#xff1a;随机分治中心点分治的期望操作次数 n≤3104n\leq 3\times 10^4n≤3104 即求点分树的 siz 之和的期望 即祖孙关系对数期望 考虑一有序点对 (u,v)(u,v)(u,v)&#xff0c;uuu 在点分树上是 vvv 祖先当且仅当 uuu 是 u∼vu\sim vu∼v 路径上第一个被选为分治中心…

译 | 使用Roslyn分析器高效编写更好的代码

原文&#xff1a;Mika Dumont翻译&#xff1a;汪宇杰Roslyn 是 .NET 编译器平台&#xff0c;即使在运行代码之前&#xff0c;它也能帮助您捕获 Bug。例如内置于 Visual Studio 中的Roslyn 拼写检查分析器。假设您正在创建一个静态方法&#xff0c;并将 static 单词拼错为 statc…

Codeforces Round #724 (Div. 2) F. Omkar and Akmar 组合数学 + 博弈

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先我们先来研究一下这个游戏&#xff0c;手画几个会惊奇的发现&#xff0c;后手这个b怎么怎么画都赢啊&#xff1f;&#xff1f;&#xff1f;对&#xff0c;没错&#xff0c;就是怎么画都…

Visual Studio 支持 Java?谣言止于智者

在去年的某一段时间内&#xff0c;我们经常会看到标题类似为“宇宙第一 IDE Visual Studio 支持 Java 了&#xff01;”的文章。当时&#xff0c;韩老师看到这篇文章也是一脸懵逼。公司里做 Java 的开发童鞋们就坐在韩老师附近&#xff08;嗯&#xff0c;不超过五米的距离&…

【SCOI2014】方伯伯的商场之旅【数位dp】【单峰函数】

题意&#xff1a;给定 l,r,kl,r,kl,r,k &#xff0c;对于一个 kkk 进制数&#xff0c;将数码看成这个位置的石子个数&#xff0c;每将一个石子移动 111 的距离需要 111 的代价。求 [l,r][l,r][l,r] 中的所有数在 kkk 进制下将石子集中在一个位置的最小代价之和。 l≤r≤1015,k≤…

Codeforces Round #597 (Div. 2) D. Shichikuji and Power Grid 最小生成树 + 虚拟点

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 看错题导致误入歧途&#xff0c;如果能早点看见翻译也不至于一天多也没想出来。 求联通的最小代价&#xff0c;自然的想到了能不能建边跑最小生成树。 对于两点之间比较好弄&#xff0c;直…

【NOI2007】货币兑换【任意坐标斜率优化】【CDQ分治】

题意&#xff1a;有 A&#xff0c;B 两种金券&#xff0c;给出 nnn 天内分别的单位价格和可以购买的数量的比例。开始有 SSS 元&#xff0c;求 nnn 天后最多能有多少元。 提示&#xff1a;每次操作一定全买全卖 n≤105n\leq 10^5n≤105 设 fnf_nfn​ 表示第 nnn 天结束后手上…

小白开学Asp.Net Core 《八》

小白开学Asp.Net Core 《八》— — .Net Core 数据保护组件1、背景我在搞&#xff08;https://github.com/AjuPrince/Aju.Carefree&#xff09;这个开源项目的时候&#xff0c;想做一些防止恶意攻击的小功能&#xff08;如果 我通过页面 /Dome/GetData?id123&#xff0c;那是不…

你必须知道的Docker镜像仓库的搭建

近期工作中发现用到的容器镜像越来越多&#xff08;不多的时候没考虑过镜像仓库的问题&#xff09;&#xff0c;同一个容器镜像也存在多个版本&#xff0c;那么镜像仓库的搭建需求就涌现出来&#xff0c;本文就目前的几个常用镜像仓库的搭建进行介绍&#xff0c;我们可以根据需…

nowcoder 河 我 车 题 错 天 乐 赛 倍增处理

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 倍增经典题了&#xff0c;考虑维护一个前缀和sum[i]sum[i]sum[i]&#xff0c;对于iii位置我们肯定是贪心的选能到的最远位置&#xff0c;即第一个>sum[i−1]k>sum[i-1]k>sum[i−1…

狄利克雷生成函数

一个没啥用的东西&#xff0c;权当加深对数论函数和狄利克雷卷积的理解。 定义 序列 {f1,f2,…}\{f_1,f_2,\dots\}{f1​,f2​,…} 的狄利克雷生成函数 DGF 定义为 f(s)∑n1∞fnnsf(s)\sum_{n1}^{\infin}\frac{f_n}{n^s}f(s)n1∑∞​nsfn​​ 这里不要求 fff 是积性函数。 显…

ASP.NET Core Web Api之JWT(一)

最近沉寂了一段&#xff0c;主要是上半年相当于休息和调整了一段时间&#xff0c;接下来我将开始陆续学习一些新的技能&#xff0c;比如Docker、Jenkins等&#xff0c;都会以生活实例从零开始讲解起&#xff0c;到时一并和大家分享和交流。接下来几节课的内容将会讲解JWT&#…