抽代的成分远远大于OI的成分
首先把一个点定为原点,然后我们发现如果我们不旋转此时答案就是所有位置的gcd
如果要旋转怎么办,我们考虑把我们选定的网格边连同方向和大小看做单位向量→e
那么此时我们把坐标系变成复平面,每个点都可以表示成(a+bi)→e的形式
当a,b均为整数时,它其实是个高斯整数的形式,那么我们可以把带余除法推广到高斯整数环
一些具体的姿势详见:高斯整数余数的一个问题
然后直接推广出辗转相除即可,利用不等式放缩容易得出此时一定是最优的
#include<cstdio>
#include<iostream>
#define int long long
#define RI register int
#define CI const int&
using namespace std;
const int N=100005;
inline int Fix(int x,int y)
{
int f=1; if (x<0) x=-x,f=-1; return (x%y*2<y)?((x/y)*f):((x/y+1)*f);
}
struct Complex
{
int x,y;
inline Complex(CI X=0,CI Y=0)
{
x=X; y=Y;
}
friend inline Complex operator + (const Complex& A,const Complex& B)
{
return Complex(A.x+B.x,A.y+B.y);
}
friend inline Complex operator - (const Complex& A,const Complex& B)
{
return Complex(A.x-B.x,A.y-B.y);
}
friend inline Complex operator * (const Complex& A,const Complex& B)
{
return Complex(A.x*B.x-A.y*B.y,A.y*B.x+A.x*B.y);
}
friend inline Complex operator / (const Complex& A,const Complex& B)
{
Complex iB=Complex(B.x,-B.y),ret=A*iB; int dv=B.x*B.x+B.y*B.y;
return Complex(Fix(ret.x,dv),Fix(ret.y,dv));
}
friend inline Complex operator % (const Complex& A,const Complex& B)
{
return A-((A/B)*B);
}
}a[N],g; int n,mix,mxx,miy,mxy;
inline Complex gcd(const Complex& A,const Complex& B)
{
if (!B.x&&!B.y) return A; return gcd(B,A%B);
}
signed main()
{
RI i; for (scanf("%lld",&n),i=1;i<=n;++i) scanf("%lld%lld",&a[i].x,&a[i].y);
for (i=2;i<=n;++i) a[i]=a[i]-a[1]; for (g=a[2],i=3;i<=n;++i) g=gcd(g,a[i]);
for (i=2;i<=n;++i) a[i]=a[i]/g,mix=min(mix,a[i].x),mxx=max(mxx,a[i].x),
miy=min(miy,a[i].y),mxy=max(mxy,a[i].y); return printf("%lld",max(mxx-mix,mxy-miy)),0;
}
__EOF__
【推荐】还在用 ECharts 开发大屏?试试这款永久免费的开源 BI 工具!
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· AI与.NET技术实操系列:使用Catalyst进行自然语言处理
· 分享一个我遇到过的“量子力学”级别的BUG。
· Linux系列:如何调试 malloc 的底层源码
· AI与.NET技术实操系列:基于图像分类模型对图像进行分类
· go语言实现终端里的倒计时
· 物流快递公司核心技术能力-海量大数据处理技术
· 几个技巧,教你去除文章的 AI 味!
· 系统高可用的 10 条军规
· 如何在 Github 上获得 1000 star?
· 关于普通程序员该如何参与AI学习的三个建议以及自己的实践
2017-12-22 Luogu P1082 同余方程
2017-12-22 EZ 2017 12 17初二初三第一次膜你赛