对RSA的攻击

 第1关：计算φ(n)

根据提示，补全右侧编辑器中 Begin-End 区间的代码，给定n和ϕ(n)，要求输出p和q，即n的两个因子。具体要求如下：

• 从后台获取两个数字n和ϕ(n)，输出n的两个因子p和q，输出要求p<q。

• #include<bits/stdc++.h>
  using namespace std;
  #define ll  long  long
  //在下面Begin和End之间补全代码，输出相应的结果
  int main()
  {
      ll n,phi;
      cin>>n>>phi;
      /*********** Begin ***********/
  
      for (ll p = 2; p * p <= n; ++p) {
          if (n % p == 0) {
              ll q = n / p;
              if ((p - 1) * (q - 1) == phi) {
                  cout << p << " " << q << endl;
                  return 0;
              }
          }
      }
  
      /*********** End ***********/
      return 0;
  }

第2关：解密指数 

根据提示，补全右侧编辑器中 Begin-End 区间的代码，给定n和1模n的四个平方根a,b,c,d，求出n的两个因子。具体要求如下：

• 从后台获取n和a，b，c，d，求出n的两个因子。 按照从小到大的顺序输出结果，即p<q。

#include<bits/stdc++.h>
using namespace std;
#define ll  long  long
//在下面Begin和End之间补全代码，输出相应的结果
ll gcd(ll a, ll b) {
    if (b == 0) return a;
    return gcd(b, a % b);
}
int main()
{
    ll n,a,b,c,d;
    cin>>n;
    cin>>a>>b>>c>>d;
    /*********** Begin ***********/
    vector<ll> roots = {a, b, c, d};

    for (int i = 0; i < 4; ++i) {
        ll x = gcd(roots[i] + 1, n);
        if (x != 1 && x != n) {
            ll y = n / x;
            if (x > y) swap(x, y);
            cout << x << " " << y << endl;
            break;
        }
    }
    /*********** End ***********/
    return 0;
}

第3关：Wiener的低解密指数攻击

根据提示，补全右侧编辑器中 Begin-End 区间的代码，给定a,b,c,d，判断c/d是否是a/b的收敛子。具体要求如下：

• 从后台获取四个整数a,b,c,d，判断c/d是否是a/b的收敛子，是的话输出YES,否则输出NO。

• #include<bits/stdc++.h>
  using namespace std;
  #define ll  long  long
  //在下面Begin和End之间补全代码，输出相应的结果
  // 辗转相除法求最大公约数
  int gcd(int a, int b) {
      if (b == 0)
          return a;
      return gcd(b, a % b);
  }
  
  // 计算分数的连分数展开
  std::vector<int> continuedFractionExpansion(double a, double b) {
      std::vector<int> cont_frac;
      while (std::fabs(b) > 1e-9) {
          int q = static_cast<int>(a / b);
          cont_frac.push_back(q);
          double temp = a;
          a = b;
          b = temp - q * b;
      }
      return cont_frac;
  }
  
  // 检查c/d是否是a/b的收敛子
  bool isConvergent(double a, double b, double c, double d) {
      std::vector<int> cont_frac_a = continuedFractionExpansion(a, b);
      std::vector<int> cont_frac_c = continuedFractionExpansion(c, d);
  
      int len_a = cont_frac_a.size();
      int len_c = cont_frac_c.size();
  
      if (len_c > len_a)
          return false;
  
      for (int i = 0; i < len_c; ++i) {
          if (cont_frac_a[i] != cont_frac_c[i])
              return false;
      }
  
      return true;
  }
  
  int main()
  {
      double a,b,c,d;
      cin>>a>>b>>c>>d;
      /*********** Begin ***********/
      if (isConvergent(a, b, c, d))
          std::cout << "YES";
      else
          std::cout << "NO";
  
      /*********** End ***********/
      return 0;
  }