bzoj 1041 圆上的整点

最新推荐文章于 2022-08-21 02:26:27 发布

转载最新推荐文章于 2022-08-21 02:26:27 发布 · 614 阅读

Brush Mode 专栏收录该内容

17 篇文章

订阅专栏

本文介绍了一种计算单位圆内整数坐标点数量的方法，通过数学推导得出核心公式，并利用C++编程实现算法。该算法首先考虑x、y大于0的情况，通过枚举2r的约数和变量u来找出所有符合条件的整数坐标，最终结果乘以4再加上坐标轴上的点数。

这里先只考虑x，y都大于0的情况

如果x^2+y^2=r^2，则(r-x)(r+x)=y*y

令d=gcd(r-x,r+x),r-x=d*u^2,r+x=d*v^2,显然有gcd(u,v)=1且u<v

有2r=d*(u^2+v^2)，y=d*u*v，x=d(v^2-u^2)/2

枚举2r的约数d，再花费sqrt(2r/d)的时间枚举u，求出v=sqrt(2r/d-u^2)然后判断gcd(u,v)=1

最后结果乘以4（四个象限）+4（坐标轴上）即可

/**************************************************************
    Problem: 1041
    User: ******
    Language: C++
    Result: Accepted
    Time:92 ms
    Memory:1284 kb
****************************************************************/
 
#include <vector>
#include <list>
#include <limits.h>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <string.h>
#include <stdlib.h>
#include <cassert>
 
using namespace std;
 
int main() {
    long long r;
    cin >> r;
    long long d, a, b;
    long long ans = 0;
    for (d = 1; d * d <= 2*r; ++d) {
        if ((2 * r) % d) continue;
        long long t = (2 * r) / d;
        for (a = 1; a * a <= t; ++a) {
            b = sqrt(t - a * a);
            if (b * b != t - a * a) continue;
            if (a >= b) continue;
            if (__gcd(a, b) > 1) continue;
            // cout << d << " " << a << " " << b << endl;
            // cout << "X: " << - a * a * d + r << endl;
            ++ans;
        }
        if (d * d == 2 * r) continue;
        t = d;      
        for (a = 1; a * a <= t; ++a) {
            b = sqrt(t - a * a);
            if (b * b != t - a * a) continue;
            if (a >= b) continue;
            if (__gcd(a, b) > 1) continue;
            // cout << t << " " << a << " " << b << endl;
            // cout << "X: " << - a * a * (2*r/d) + r << endl;
            ++ans;
        }
    }
    ans = (ans + 1) * 4;
    cout << ans << endl;
    return 0;
}

转摘至： http://blog.youkuaiyun.com/lwfcgz/article/details/39927801