首页 > 代码库 > BZOJ1041 [HAOI2008]圆上的整点
BZOJ1041 [HAOI2008]圆上的整点
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。
本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!
题目链接:BZOJ1041
正解:数学
解题报告:
参考博客:http://www.cppblog.com/zxb/archive/2010/10/18/130330.html。
讲的比较清楚,公式推导就不赘述了。
//It is made by ljh2000 //有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。 #include <iostream> #include <cstdlib> #include <cstring> #include <cstdio> #include <cmath> #include <algorithm> #include <ctime> #include <vector> #include <queue> #include <map> #include <set> #include <string> #include <complex> #include <bitset> using namespace std; typedef long long LL; typedef long double LB; typedef complex<double> C; const double pi = acos(-1); LL r,d,ans,rr; inline LL gcd(LL x,LL y){ if(y==0) return x; return gcd(y,x%y); } inline int getint(){ int w=0,q=0; char c=getchar(); while((c<‘0‘||c>‘9‘) && c!=‘-‘) c=getchar(); if(c==‘-‘) q=1,c=getchar(); while (c>=‘0‘&&c<=‘9‘) w=w*10+c-‘0‘,c=getchar(); return q?-w:w; } inline void solve(LL x){ LL u,v,lin; for(u=1;u*u<=x;u++) { lin=x-u*u; v=sqrt(lin); if(u>=v) break;//前提 if(v*v==lin && gcd(u,v)==1/*!!!*/) ans++; } } inline void work(){ r=getint(); rr=r*2; ans=1; for(d=1;d*d<=rr;d++) { if(rr%d!=0) continue;//!!! solve(rr/d); if(d*d!=rr) solve(d); } cout<<ans*4; } int main() { work(); return 0; } //有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
BZOJ1041 [HAOI2008]圆上的整点
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。