AcWing寒假每日一题2023-1.20解密

宣传一下算法提高课整理 <—

csdn

题目描述

给定一个正整数 k k k,有 k k k 次询问,每次给定三个正整数 n i , e i , d i n_i,e_i,d_i ni,ei,di,求两个正整数 p i , q i p_i,q_i pi,qi,使 n i = p i × q i n_i=p_i×q_i ni=pi×qi e i × d i = ( p i − 1 ) ( q i − 1 ) + 1 e_i×d_i=(p_i−1)(q_i−1)+1 ei×di=(pi1)(qi1)+1

输入格式

第一行一个正整数 k k k,表示有 k k k 次询问。

接下来 k k k 行,第 i i i 行三个正整数 n i , d i , e i n_i,d_i,e_i ni,di,ei

输出格式

输出 k k k 行,每行两个正整数 p i , q i p_i,q_i pi,qi 表示答案。

为使输出统一,你应当保证 p i ≤ q i p_i≤q_i piqi

如果无解,请输出 NO

数据范围

以下记 m = n − e × d + 2 m=n−e×d+2 m=ne×d+2

保证对于 100 100 100% 的数据, 1 ≤ k ≤ 1 0 5 1≤k≤10^5 1k105,对于任意的 1 ≤ i ≤ k 1≤i≤k 1ik 1 ≤ n i ≤ 1 0 18 1≤n_i≤10^{18} 1ni1018 1 ≤ e i × d i ≤ 1 0 18 1≤e_i×d_i≤10^{18} 1ei×di1018 1 ≤ m ≤ 1 0 9 1≤m≤10^9 1m109
a

样例输入

10
770 77 5
633 1 211
545 1 499
683 3 227
858 3 257
723 37 13
572 26 11
867 17 17
829 3 263
528 4 109

样例输出

2 385
NO
NO
NO
11 78
3 241
2 286
NO
NO
6 88

思路1:二分

先看数据范围: 1 ≤ k ≤ 1 0 5 1≤k≤10^5 1k105 1 ≤ n i ≤ 1 0 18 1≤n_i≤10^{18} 1ni1018 1 ≤ e i × d i ≤ 1 0 18 1≤e_i×d_i≤10^{18} 1ei×di1018 1 ≤ m ≤ 1 0 9 1≤m≤10^9 1m109
那么1e18的数据一定是要开long long的。
仔细读题可以发现,数据范围处给了一个 令 m = n − e × d + 2 令m=n-e×d+2 m=ne×d+2,我们思考:这是否会是解题的突破口呢?作者在比赛的时候没有看到m,这道题用暴力做的,只拿了20分

对于此等式: e i × d i = ( p i − 1 ) ( q i − 1 ) + 1 e_i×d_i=(p_i−1)(q_i−1)+1 ei×di=(pi1)(qi1)+1

化简后: e i × d i = p q − ( p + q ) + 2 e_i×d_i=pq-(p+q)+2 ei×di=pq(p+q)+2

∵ n = p × q ∵n=p×q n=p×q

∴ e i × d i = n − ( p + q ) + 2 ∴e_i×d_i=n-(p+q)+2 ei×di=n(p+q)+2

移项,得 p + q = n − e i × d i + 2 p+q=n-e_i×d_i+2 p+q=nei×di+2

m = p + q m=p+q m=p+q

联立,得
{ n = p × q m = p + q \left\{\begin{matrix} n=p×q \\ m=p+q \end{matrix}\right. {n=p×qm=p+q
a
a
(图片来自@JcWing)

思路1时间复杂度

共k组数据,每组数据二分范围1~ n / 2 n/2 n/2,因此时间复杂度为 O ( k log ⁡ m ) O(k\log m) O(klogm)

可以AC

思路1-AC Code

C + + C++ C++

#include <iostream>
#include <cstring>

using namespace std;

typedef long long LL;

LL k, n, m, e, d;
LL p, q;

int main()
{
	scanf("%lld", &k);
	
	while (k -- )
	{
		scanf("%lld%lld%lld", &n, &d, &e);
		m = n - e * d + 2;
		
		LL l = 1, r = m >> 1;
		while (l < r)
		{
			LL mid = l + r >> 1;
			if (mid * (m - mid) >= n) r = mid;
			else l = mid + 1;
		}
		
		p = l, q = m - l;
		if (p * q == n) printf("%lld %lld\n", p, q);
		else puts("NO");
	}
	
	return 0;
}

思路2:逆造方程

{ n = p × q m = p + q \left\{\begin{matrix} n=p×q \\ m=p+q \end{matrix}\right. {n=p×qm=p+q

运用韦达逆定理逆造方程,则

p 2 + m p + n = 0 p^2+mp+n=0 p2+mp+n=0

运用求根公式 − b ± b 2 − 4 a c 2 a \frac{-b\pm \sqrt{b^2-4ac} }{2a} 2ab±b24ac

得到 △ = m 2 − 4 n △=m^2-4n =m24n

p , q p,q p,q分别为 m − △ m-\sqrt{\triangle } m m + △ m+\sqrt{\triangle } m+

思路2时间复杂度

k k k次循环,因此是 O ( k ) O(k) O(k)

可真挺ok的

思路2-AC Code

C + + C++ C++

#include <iostream>
#include <cmath>

using namespace std;

typedef long long LL;

int main()
{
    int k;
    scanf("%d", &k);

    while (k -- )
    {
        LL n, d, e;
        scanf("%lld%lld%lld", &n, &d, &e);
        LL m = n - e * d + 2;
        LL dt = m * m - 4 * n;
        LL r = sqrt(dt);

        if (dt < 0 || r * r != dt) puts("NO");
        else printf("%lld %lld\n", (m - r) / 2, (m + r) / 2);
    }

    return 0;
}

最后,如果觉得对您有帮助的话,点个赞再走吧!

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

星河依旧长明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值