【BZOJ】1041: [HAOI2008]圆上的整点(几何)

http://www.lydsy.com:808/JudgeOnline/problem.php?id=1041

【BZOJ】1041: [HAOI2008]圆上的整点(几何)

所谓的神题,我不会,直接题解。。看了半天看懂题解了。详见hzwer博客

这题呢,我只能吸收些思想,即,当我们要找合法解的时候,我们可以深究它的性质,然后用性质来判定是否存在合法解。

此神题直接看题解打码。

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define read(a) a=getnum()
#define print(a) printf("%d", a)
#define debug(a) printf("%lld\n", a)
inline int getnum() { int ret=0; char c; for(c=getchar(); c<'0' || c>'9'; c=getchar()); for(; c>='0' && c<='9'; c=getchar()) ret=ret*10+c-'0'; return ret; }
typedef long long ll;
ll gcd(ll a, ll b) { return b?gcd(b, a%b):a; } inline bool check(ll A, ll B) {
if(((ll)sqrt(B)*(ll)sqrt(B))==B && A!=B)
if(gcd(A, B)==1) return true;
return false;
} int main() {
int ans=0;
ll d, d2, r, r2;
scanf("%lld", &r);
r2=r<<1;
ll m=sqrt(r2);
ll a;
for(d=1; d<=m; ++d) {
if(!(r2%d)) {
d2=d<<1;
for(a=1; a<=(ll)sqrt(r2/d2); ++a)
if(check(a*a, r2/d-a*a)) ++ans;
if(d!=r2/d) {
for(a=1; a<=(ll)sqrt(d/2); ++a)
if(check(a*a, d-a*a)) ++ans;
}
}
}
printf("%lld\n", (ll)(ans*4+4));
return 0;
}

Description

求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数。

Input

r

Output

整点个数

Sample Input

4

Sample Output

4

HINT

n<=2000 000 000

Source

上一篇:Java 并发专题 : Executor详细介绍 打造基于Executor的Web服务器


下一篇:ASP.NET MVC应用程序执行过程分析