Codeforces-1487 D. Pythagorean Triples（数学）

最新推荐文章于 2021-11-10 15:46:19 发布

tomjobs

最新推荐文章于 2021-11-10 15:46:19 发布

阅读量382

点赞数 1

分类专栏： # codeforces

本文链接：https://blog.csdn.net/tomjobs/article/details/113825440

版权

codeforces 专栏收录该内容

553 篇文章 5 订阅

订阅专栏

在这里插入图片描述
思路：
数学题，打表可以找出规律，不过我们也可以推一下。

$c^2 = a^2 + b^2$
$c = a^2 - b$
$\Rightarrow$ $c^2 = b^2 + c + b$
$\Rightarrow$ $c (c - 1) = b (b + 1)$
因为 $c, c - 1$ 互质， $b, b + 1$ 互质，所以要么 $c = b$ ，要么 $c = b + 1$
而 $c > b$
所以 $c = b + 1$

推出 $b=(a^2-1)/2$
推出 $c=(a^2-1)/2+1$

易得 $a$ 必须得为大于3的奇数， $c$ 才能有意义。而我们已经知道 $c$ 的最大范围 $n$ ，所以只需要算出存在多少 $a^2-1)/2+1≤n$ ， $a = 2 k + 1, k \geq 1$ 。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <map>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 7;
const int mod = 1e9 + 7;
typedef long long ll;
ll c[maxn];
int cnt;
void init() {
    cnt = 0;
    for(int a = 3;a <= 100000;a += 2) {
        c[++cnt] = (1ll * a * a - 1) / 2 + 1;
    }
}
int main() {
    init();
    int T;scanf("%d",&T);
    while(T--) {
        int n;scanf("%d",&n);
        int ans = upper_bound(c + 1,c + 1 + cnt,n) - (c + 1);
        printf("%d\n",ans);
    }
    return 0;
}