#CSPJ202202. 解密(AC代码)

2022年 CSP_J 第二轮 第二题 

#CSPJ202202. [CSP-J 2022] 解密

题目描述

给定一个正整数 k,有 k 次询问,每次给定三个正整数 n_i, e_i, d_i,求两个正整数 p_i, q_i,使 n_i = p_i \times q_ie_i \times d_i = (p_i - 1)(q_i - 1) + 1

输入格式

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

接下来 k 行,第 i 行三个正整数 n_i, e_i, d_i

输出格式

输出 k 行,每行两个正整数 p_i, q_i 表示答案。

为使输出统一,你应当保证 p_i \leq q_i

如果无解,请输出 NO

样例 #1

样例输入 #1

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

样例输出 #1

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

提示

【样例 #2】

见附件中的 decode/decode2.in 与 decode/decode2.ans

【样例 #3】

见附件中的 decode/decode3.in 与 decode/decode3.ans

【样例 #4】

见附件中的 decode/decode4.in 与 decode/decode4.ans

【数据范围】

以下记 m = n - e \times d + 2

保证对于 100% 的数据,1 \leq k \leq {10}^5,对于任意的 1 \leq i \leq k1 \leq n_i \leq {10}^{18}1 \leq e_i \times d_i \leq {10}^{18}1 \leq m \leq {10}^9

测试点编号k≤n≤m≤特殊性质
110^310^3保证有解
2
310^96\times 10^4保证有解
4
510^9保证有解
6
710^510^{18}保证若有解则 p=q
8保证有解
9
10

【题意】

给三个正整数,n,e,d,求两个正整数p,q,使n=p\times qe\times d=(p-1)\times (q-1) +1

【题解】

即解二元一次方程:

e\times d=(p-1)\times (q-1) +1展开一下:

e\times d = (p-1)\times (q-1)+1 =p\times q-(p+q) +2

已知p\times q=n,是个常数,那么对于p+q=n-e\times d+2,这个等式右边都是已知的数,那么我们令m=n-e\times d +2

也就是到现在,我们有:m=p+qn=p\times q,可得一元二次方程:

p^2-m\times p+n=0

然后我们用一元二次方程的求根公式即可求得pq

p=\frac{m-\sqrt{m^2-4n}}{2}, q=\frac{m+\sqrt{m^2-4n}}{2}

另外需要注意的是:输入数据中保证了1\leqslant m\leqslant 10^9 ,那么根号底下的数会超过int但不会超过long long,我们需要自己编写一个二分开根号的函数。
并不是所有数据都保证有解,如果要拿满分的话,别忘了判断解的合法性(可以在求得pq的时候,来再去校验一下,来简单判断)。

【代码】 

#include <bits/stdc++.h>
using namespace std;
long long solve(long long x)
{
    long long left=0,right=x;
    while(left<right){
        long long mid=(left+right+1)/2;
        if((double)mid*mid<=x){
            left=mid;
        }else{
            right=mid-1;
        }
    }
    return left;
}
int main()
{
    //freopen("decode.in","r",stdin);
    //freopen("decode.out","w",stdout);
    long long t;
    cin>>t;
    while(t--){
        long long n,e,d;
        cin>>n>>e>>d;
        long long m=n-e*d+2;
        long long f=solve(m*m-4*n);
        long long p=0,q=0;
        if (f*f==(m*m-4*n)){
            p=(m-f)/2;
            q=(m+f)/2;
        }
        if((p+q)!=m || p*q!=n){
            cout<<"NO"<<endl;
        }else{
            cout<<p<<' '<<q<<endl;
        }
    }
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值