2021牛客暑期多校训练营3-E

2021牛客暑期多校训练营3-E

题意

大致题意为 满足x^2 + y^2 整除 x*y+1的x与y的对数

思考

要满足题意 可以列出方程 x^2 + y^2 =k(xy+1)
首先我们把x进行固定 即视为一个常数 方程变为关于y的二元方程
韦达定理 有
y1+y2=-b/a=kx ; y1*y2=c/a=x^2 + k
设 y=kx 代入原方程解出k=x^2
然后一直递推下去 讲y放进vector

代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<iomanip>
#include<cstring>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<vector>
#define maxn 100005
#define N 1000005
#define ninf 0xc0c0c0c0
#define ll long long
#define ull __int128
#define endl '\n'
#define read(a) {char wmc;while((wmc=getchar())>47) a=a*10+(wmc^48);}
using namespace std;
typedef pair<int,int> P;
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	
	vector<ll>q;
	q.push_back(1);
	
	
	for(ll i=2;i<=1e6;i++){
		ull a=i;
		ull b=i*i*i;
		ull k=i*i;
		ull res;
		while(b<=1e18){
			q.push_back(b);
			res=k*b-a;
			a=b;
			b=res;
		}
	}
	sort(q.begin(),q.end());
	int t;
	cin>>t;
	while(t--)
	{
		
		ll n;
		cin>>n;
		ll ans=upper_bound(q.begin(),q.end(),n)-q.begin();
		cout<<ans<<endl;
	}
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值