gcd_1.c 文件参考

#include "../../../include/lammp/numth.h"
#include "../../../include/lammp/lmmpn.h"

gcd_1.c 的引用(Include)关系图:

浏览源代码.

函数
mp_limb_t	lmmp_gcd_11_ (mp_limb_t u, mp_limb_t v)
	计算 [numa,na] 在B^n 下的逆元
mp_limb_t	lmmp_gcd_1_ (mp_srcptr up, mp_size_t un, mp_limb_t vlimb)
	计算两个无符号整数的最大公约数

函数说明

lmmp_gcd_11_()

mp_limb_t lmmp_gcd_11_	(	mp_limb_t	u,
		mp_limb_t	v
	)

计算 [numa,na] 在B^n 下的逆元

参数

dst	结果指针（长度为 n 个limb）
numa	待求逆元指针（长度为 na 个limb）
na	待求逆元的 limb 长度
n	结果的 limb 长度

警告

n>=na>0, numa!=NULL, dst!=NULL, numa[0]%2==1, sep(dst,numa)

返回

dst 的实际 limb 长度

计算两个无符号整数的最大公约数

参数

u	第一个无符号整数
v	第二个无符号整数

返回

最大公约数

警告

u!=0, v!=0

在文件 gcd_1.c 第 10 行定义.

                                                 {
    lmmp_param_assert(u > 0 && v > 0);
    int count, k;
    k = lmmp_tailing_zeros_(u | v);
    u >>= lmmp_tailing_zeros_(u);
    v >>= lmmp_tailing_zeros_(v);
    while (u != v) {
        if (u > v) {
            u -= v;
            count = lmmp_tailing_zeros_(u); 
            u >>= count;
        } else {
            v -= u;
            count = lmmp_tailing_zeros_(v);
            v >>= count;
        }
    }
    return u << k;
}

引用了 k, lmmp_param_assert , 以及 lmmp_tailing_zeros_().

被这些函数引用 lmmp_gcd_1_(), lmmp_gcd_22_() , 以及 lmmp_gcd_lehmer_().

函数调用图:

这是这个函数的调用关系图:

lmmp_gcd_1_()

mp_limb_t lmmp_gcd_1_	(	mp_srcptr	up,
		mp_size_t	un,
		mp_limb_t	vlimb
	)

计算两个无符号整数的最大公约数

参数

up	第一个无符号整数指针
un	第一个无符号整数的 limb 长度
v	第二个无符号整数

警告

v!=0, up!=NULL, un>0

返回

最大公约数

在文件 gcd_1.c 第 30 行定义.

                                                                   {
    lmmp_param_assert(un > 0);
    lmmp_param_assert(vlimb > 0);
    mp_limb_t ulimb;
    if (un == 1) {
        ulimb = up[0];
    }
    else {
        ulimb = lmmp_mod_1_(up, un, vlimb);
    }
    if (ulimb == 0)
        return vlimb;
    else
        return lmmp_gcd_11_(ulimb, vlimb);
}

引用了 lmmp_gcd_11_(), lmmp_mod_1_() , 以及 lmmp_param_assert.

被这些函数引用 lmmp_gcd_lehmer_().

函数调用图:

这是这个函数的调用关系图: