BZOJ 4544: 椭圆上的整点
Sol
数学.
跟圆上的整点一样...TA写了个积性函数的算法...以后再说吧...
\(x^2+3y^2=r^2\)
\(3y^2=r^2-x^2\)
\(3y^2=(r-x)(r+x)\)
\(y^2=\frac{1}{3}(r-x)(r+x)\)
\(d=(r-x)(r+x)\)
\(r-x=3du^2,r+x=dv^2\) 这里 \(r-x\) 和 \(r+x\) 并没有什么区别.
\(2r=d(3u^2+v^2)\)
枚举 \(d\) 和 \(u\)
感觉复杂度是\(O(n^{\frac{3}{4}})\)
但是可以跑最大数据的说.
Code
/**************************************************************Problem: 4544User: BeiYuLanguage: C++Result: AcceptedTime:8568 msMemory:1300 kb
****************************************************************/#include<cstdio>
#include<cmath>
#include<algorithm>
#include<vector>
#include<utility>
#include<iostream>
using namespace std;typedef long long LL;
#define debug(a) cout<<#a<<"="<<a<<" "
#define mpr(a,b) make_pair(a,b)LL T,r,n,ans;inline LL in(LL x=0,char ch=getchar()){ while(ch>'9'||ch<'0') ch=getchar();while(ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();return x; } vector<pair<LL,LL> > p;LL calc(LL d){LL res=0,m=n/d;
// cout<<"*************"<<endl;
// debug(m),debug(d);cout<<endl;for(LL u=1,v;u*u*3<=m;u++){v=sqrt(m-3*u*u+0.5);
// debug(u),debug(v),debug(3*v*v+u*u),cout<<endl;
// if(u>v) break;if(v*v+u*u*3==m&&__gcd(v*v,u*u*3)==1) res++;
// cout<<"get!",debug(d*u*u*3),debug(d*v*v),debug(d*u*u*3+d*v*v)<<endl;
// p.push_back(mpr(d*u*u*3,d*v*v));}return res;
}
int main(){
// freopen("in.in","r",stdin);for(T=in();T--;){r=in(),n=r<<1,ans=0;for(LL d=1;d*d<=n;d++) if(n%d==0){if(d*d==n) ans+=calc(d);else ans+=calc(d)+calc(n/d);}cout<<ans*4+2<<endl;
// sort(p.begin(),p.end());
// for(int i=0;i<p.size();i++) cout<<p[i].first<<" "<<p[i].second<<endl;}return 0;
}
转载于:https://www.cnblogs.com/beiyuoi/p/5857781.html
BZOJ 4544: 椭圆上的整点相关推荐
- [BZOJ 1041] 圆上的整点
Description 求一个给定的圆( x2+y2=r2 x 2 + y 2 = r 2 x^2+y^2=r^2),在圆周上有多少个点的坐标是整数 Input 只有一个正整数 n n n,n≤200 ...
- BZOJ 1041 圆上的整点 数学
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=1041 题目大意: 求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是 ...
- bzoj4544: 椭圆上的整点
链接 http://www.lydsy.com/JudgeOnline/problem.php?id=4544 题解 "我真傻,真的,"wyw抬起没有神采的眼睛,接着说." ...
- [bzoj4544]椭圆上的整点 解题报告
先来看傻逼的做法: 稍微化下式子. x2+3y2=n23y2=(n+x)(n−x) x^2+3y^2=n^2\\ 3y^2=(n+x)(n-x) 令 g=(n+x,n−x)g=(n+x,n-x),则有 ...
- bzoj 1041: [HAOI2008]圆上的整点
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 2838 Solved: 1238 [Submit][S ...
- bzoj千题计划127:bzoj1041: [HAOI2008]圆上的整点
http://www.lydsy.com/JudgeOnline/problem.php?id=1041 设 X>0 ,Y>0 X^2 + Y^2 = R^2 X^2 = R^2-Y^2 ...
- 已知长短轴求椭圆上任意一点的坐标_高中数学必修2:平面解析几何——椭圆(经典习题)...
今天接着给大家分享关于高中数学必修2平面解析几何中椭圆知识点讲解,从三个方面进行讲解:基础梳理.要点整理.经典高考习题解题过程及答案. 一.基础梳理 1.椭圆的定义 2.椭圆的标准方程和几何性质 要点 ...
- OpenCv:椭圆上点的计算方程
椭圆 椭圆(Ellipse)是平面内到定点F1.F2的距离之和等于常数(大于|F1F2|)的动点P的轨迹,F1.F2称为椭圆的两个焦点.其数学表达式为: ...
- 求圆和椭圆上任意角度的点的坐标
圆上任意角度的点的坐标 如上图,给定圆心(Cx,Cy),半径为R, 求θ\thetaθ对应的点的坐标? 此处θ\thetaθ是相对于水平轴的角度. 显然我们可以使用极坐标转换来求: {px=Cx+Rc ...
最新文章
- http://blog.51cto.com/wbb827/1175634
- 如何搞垮一个技术大牛?
- 软件维护费 一般按照多少比例_2021年国家电网招聘:笔试进面一般是多少比例?...
- 使用openSSL构造一个支持https的nodejs服务器
- python网络编程自学_五分钟搞定Python网络编程实现TCP和UDP连接
- git 历史版本导出_Git基础知识(九)
- 计算机快捷键屏幕录制,屏幕录制软件哪个好用?设置电脑录屏快捷键
- Python常见陷阱
- 汪学明导师—商业模式创新与转型专家
- 计算机网络自顶向下方法 第四章 作业习题答案
- 读书笔记《计算机是怎样跑起来的》
- 腾讯云超级播放器实现无直播占位图,有直播播放视频
- 网络云存储技术Windows server 2012 (项目六 存储服务器的数据备份与还原 windows Server Backup)
- 树莓派入门(五)—— 树莓派4B之初体验
- 奥布莱恩杯尘埃落定 人工智能立功了!
- win7鼠标指针主题包_双飞燕FB35双模蓝牙鼠标评测 终于能让手放松了
- oracle减法函数mius_Oracle 基础知识习题175道
- 第八部分 项目资源管理
- docker基础(11):自定义tomcat镜像启动与访问
- 宽带信号的DOA估计学习笔记(四):空间分辨率