Sherwood算法求解离散对数问题

最新推荐文章于 2024-07-10 17:13:21 发布

南七行者

最新推荐文章于 2024-07-10 17:13:21 发布

阅读量858

点赞数

分类专栏：数据结构与算法文章标签：算法数据结构算法导论

本文链接：https://blog.csdn.net/weixin_37477009/article/details/109013101

版权

数据结构与算法专栏收录该内容

11 篇文章 0 订阅

订阅专栏

一、离散对数问题
1.定义
设 $a = g ^ x\;mod\;p$ ，记 $log_{g,\;p} a = x$ ，称x为a的对数(以g为底，模p)。
给定p, g, a计算x称为离散对数问题。
2.应用之一：加密
若g，p已知，x为私钥，求a公钥很容易，即解密；
若g，p已知，a为公钥，求x公钥很困难(p是很大的质数更难)，即破解困难(下文将描述)；

二、简单算法
根据费马小定理( $a^{p - 1}\;mod\;p = 1$ ，p是质数，整数a不是p的倍数)， $x = 0$ 和 $x = p - 1$ 计算 $g ^ x\;mod\;p$ 重复，最坏的情况下，x从0~p-2都不重复；所以，简单算法的思想就是遍历x，找到满足 $a = g ^ x\;mod\;p$ 的x(一般a不取0，所以遍历x从1~p-1)。

代码如下：

#include <cstdio>

int dlog(int g, int a, int p) {
    int x = 0, y = 1;
    do{
        x++;
        y *= g;
    }while((a != y % p) && (x != p));
    return x;
}

int main() {
    //测试
    int g = 2, a = 6, p = 19;
    int ans = dlog(g, a, p);
    if(ans == p) printf("无解");
    else printf("x = %d", ans);
    return 0;
}

分析：
当a = 14时， x = 7；意味着尝试了7次；
当a = 7时，x = 14；意味着尝试了14次；
也就是说，执行时间与a的取值相关。

动机：用sherwood算法使得执行时间与a无关。

三、sherwood算法
1.一般算法f(x)改造成Sherwood算法：

RH(x) {
     // 用Sherwood算法计算f(x)
     n ← length[x];             // x的size为n
     r ← uniform(An);           // 随机取一元素
     y ← u(x, r);                 //将原实例x转化为随机实例y
     s ← f(y);                  // 用确定算法求y的解s
     return v(r, s);                 // 将s的解变换为x的解
}

Sherwood算法的核是找到u，v方法
这个伪码看不明白，看下面实例

2.计算离散对数的改造思路
(1)构造 : $y=log_{g,\;p}c$
其中，
$\;mod\; p$ ( $b = g^r\;mod\;p$ )
r对应于uniform(An)，实际是从0~p-2中随机选一个数。

$log_{g,\;p} b = r$
$log_{g,\;p} a = x$

(2)推导过程
$y =log_{g,\;p} c=log_{g,\;p}(b·a\;mod\;p)=(log_{g,\;p} b + log_{g,\;p} a)\;mod\:(p - 1)=(r+x)\;mod\;(p-1)$

推导依据的定理
$log_{g,\;p}(s·t\;mod\;p)=(log_{g,\;p} s + log_{g,\;p} t)\;mod\:(p - 1)$

总结：
$y=(r+x)\;mod\;(p-1)$ ，得 $\;mod \;(p - 1)$

推导过程：
$r + x = k \cdot (p - 1) + y (k 为一常系数)$ ， $x = k \cdot (p - 1) + y - r$ , 故 $(y-r)\;mod\;p$

上述提到的：
u指的便是，把求x转化为求y，即 $y=g^c\;mod\;p$ ， $\;mod\; p$
v指的便是，已知y后求x，即 $\;mod \;(p - 1)$

3.代码

#include <cstdio>
#include <random>
#include <ctime>

int dlog(int g, int a, int p) {
    int x = 0, y = 1;
    do{
        x++;
        y *= g;
    }while((a != y % p) && (x != p));
    return x;
}

int uniform(int maxN) {
    return rand() % maxN;
}

int ModularExponent(int g, int r, int p) {
    long long sum = 1;
    for(int i = 1; i<= r; i++) sum *= g;
    return sum % p;
}


int dlogRH(int g, int a, int p) {
    int r = uniform(p - 1);
    int b = ModularExponent(g, r, p);
    int c = (b * a) % p;
    int y = dlog(g, c, p);  //这里y可能已经无解了；
    if(y == p) return -1;
    else return (y - r + p -1) % (p - 1);   //C++,y-r可能是负数，所以要加p - 1;
}


int main() {
    //测试
    srand((unsigned)time(NULL));
    int g = 2, a = 6, p = 19;
    printf("确定算法的结果：\n");
    int ans1 = dlog(g, a, p);
    if(ans1 == p) printf("无解\n");
    else printf("x = %d\n", ans1);
     printf("Sherwood算法的结果：\n");
    int ans2 = dlogRH(g, a, p);
    if(ans2 == -1) printf("无解\n");
    else printf("x = %d\n", ans2);
    return 0;
}