2021牛客多校第三场E

最新推荐文章于 2022-07-04 18:53:12 发布

Gone_Summer

最新推荐文章于 2022-07-04 18:53:12 发布

阅读量159

点赞数

分类专栏： 2021牛客多校第三场文章标签：算法

本文链接：https://blog.csdn.net/Gone_Summer/article/details/120499510

版权

2021牛客多校第三场专栏收录该内容

3 篇文章 0 订阅

订阅专栏

本文探讨了一道涉及数学规律和编程的题目，目标是找到满足x² + y² = k(xy + 1)的整数解(x, y)，其中1 < x <= y <= n。通过观察和推导，作者使用了动态规划的方法来计算每个x对应的解，并利用了lower_bound和upper_bound函数进行区间查找。最后，代码示例展示了如何在1到n的范围内计算解的数量。

摘要由CSDN通过智能技术生成

题意：找出满足x²+y²=k（xy+1）的（x，y）在1<=x<=y<=n范围内共有多少个

思路：看到范围，显然这应该是一道打表规律题

于是我们开始推公式：把x视为常量，（x，y1）为一个点，（x，y2）又是一个点，由韦达定理可以推出下面结论：

起始条件y1=0（这个点不算在内）则k=x²（这个是以后也是这样的，推出k恒为x²），y2=kx

然后把y2视为新的x，x视为y1进行下一轮推算，见下图：

因此：

之后就判断就行

代码用到了一个知识：

参考文章：关于lower_bound( )和upper_bound( )的常见用法_brandong-CSDN博客_lower_bound

代码：

#include <bits/stdc++.h>
using namespace std;
//(x,kx-y);k=x2，这是带入原方程得到的，所以一直是这个数，而其他的x和y是根据上一个点推出来的
//(x,x^3)(x^3,x)(x^3,x^5-x),原来是这样算的，k=x^2是定的
typedef long long ll;
const ll M = 1e18 + 1, N = 1e6 + 10;

ll n;
int t;
ll ans[N];
vector<ll> v;//放答案,我也不知道有多少个，正常这里可以二分，但是既然vector中有这个函数，为什么不用呢

int main()
{
	cin >> t;
	//因为第一个点是（x，x^3)，所以是x*x*x,所以数组只需要开6次方
	for (int x = 2; 1ll * x * x * x <= M; x++)x=1是无解的,但是从1开始会死循环
	{//嘶，我就想不明白了，为啥加上=M就能过，没有+就不过？？？？？、
		ans[0] = 0, ans[1] = x;
		for (int i = 2;; i++)
		{
			__int128 tem = 1ll * x * x;
			tem = tem * ans[i - 1] - ans[i - 2];
			if (tem > M)
			{
				break;
			}
			ans[i] = tem;
			v.push_back(tem);
		}
	}
	sort(v.begin(), v.end());
	while (t--)
	{
		scanf("%lld", &n);
		int cnt = upper_bound(v.begin(), v.end(), n) - v.begin() + 1;
		printf("%d\n", cnt);
	}
	return 0;
}