PazeSSH/tests/dbg_p2.c

/* dbg_p2.c —— 单测 p^2 / p 的 Knuth D 过程 */
#include "paze/bignum.h"
#include "paze/p256.h"
#include <stdio.h>

int main(void) {
    paze_bn_t p, a, q, r;
    paze_p256_load_p(&p);
    paze_bn_mul(&a, &p, &p);
    fprintf(stderr, "=== p^2/p ===\n");
    paze_bn_divmod(&q, &r, &a, &p);
    {
        char h[800];
        paze_bn_to_hex(&q, h, sizeof(h)); printf("q = %s\n", h);
        paze_bn_to_hex(&r, h, sizeof(h)); printf("r = %s\n", h);
    }
    /* 验证 q*p + r == a */
    paze_bn_t prod, chk;
    paze_bn_mul(&prod, &q, &p);
    paze_bn_add(&chk, &prod, &r);
    printf("q*p+r==a: %d\n", paze_bn_cmp(&chk, &a) == 0);
    return 0;
}