第四平方和定理,用c语言实现

1.实验题目

1.7【问题描述】

四平方和定理,又称为拉格朗日定理:每个正整数都可以表示为至多4个正整数的平方和。  如果把0包括进去,就正好可以表示为4个数的平方和。比如:5 = 0^2 + 0^2 + 1^2 + 2^2 7 = 1^2 + 1^2 + 1^2 + 2^2 ^符号表示乘方的意思)对于一个给定的正整数,可能存在多种平方和的表示法。 要求你对4个数排序: 0 <= a <= b <= c <= d  并对所有的可能表示法按 a,b,c,d 为联合主键升序排列,最后输出第一个表示法  

【输入】  一个正整数N (N<5000000)  

【输出】  4个非负整数,按从小到大排序,中间用空格分开  

输入范例1 5  

输出范例10 0 1 2  

输入范例2】773535

输出范例2】1 1 267 838 

2.需求分析

本演示程序用VC++6.0编写,完成一个数经过拉格朗日定理所得结果
 输入的形式和输入值的范围:一个正整数N (N<5000000)。在所有输入中,元素的值都是整数
 输出的形式:4个非负整数,按从小到大排序,中间用空格分开

以下是C语言的代码实现: 1. RSA算法 加密: ```c #include <stdio.h> #include <stdlib.h> #include <string.h> #include <openssl/rsa.h> #include <openssl/pem.h> int main() { // 生成RSA公私钥 RSA *key = RSA_generate_key(2048, RSA_F4, NULL, NULL); // 保存公私钥 FILE *f; f = fopen("private_key.pem", "w"); PEM_write_RSAPrivateKey(f, key, NULL, NULL, 0, NULL, NULL); fclose(f); f = fopen("public_key.pem", "w"); PEM_write_RSAPublicKey(f, key); fclose(f); // 使用公钥加密 char plaintext[] = "hello world"; size_t plaintext_len = strlen(plaintext); unsigned char ciphertext[256]; RSA_public_encrypt(plaintext_len, (unsigned char *)plaintext, ciphertext, key, RSA_PKCS1_OAEP_PADDING); return 0; } ``` 解密: ```c #include <stdio.h> #include <stdlib.h> #include <string.h> #include <openssl/rsa.h> #include <openssl/pem.h> int main() { // 读取私钥 FILE *f; f = fopen("private_key.pem", "r"); RSA *key = RSA_new(); key = PEM_read_RSAPrivateKey(f, &key, NULL, NULL); fclose(f); // 使用私钥解密 unsigned char plaintext[256]; size_t plaintext_len = RSA_private_decrypt(sizeof(ciphertext), ciphertext, plaintext, key, RSA_PKCS1_OAEP_PADDING); return 0; } ``` 2. 模重复平方算法 ```c #include <stdio.h> #include <stdint.h> uint64_t mod_exp(uint64_t base, uint64_t exponent, uint64_t modulus) { uint64_t result = 1; while (exponent > 0) { if (exponent & 1) { result = (result * base) % modulus; } exponent >>= 1; base = (base * base) % modulus; } return result; } ``` 3. 蒙哥马利算法 ```c #include <stdio.h> #include <stdint.h> uint64_t montgomery_reduction(uint64_t x, uint64_t r, uint64_t n, uint64_t n_inv) { // Step 1 uint64_t m = (x * r) % n; // Step 2 uint64_t t = (x + m * n) * n_inv % r; // Step 3 if (t >= r) { return t - r; } else { return t; } } uint64_t mod_exp_montgomery(uint64_t base, uint64_t exponent, uint64_t modulus) { // Step 1 uint64_t r = 1; while (r < modulus) { r <<= 1; } // Step 2 uint64_t n_inv = -modulus % r; // Step 3 base = (base * r) % modulus; // Step 4 uint64_t x = r % modulus; while (exponent > 0) { if (exponent & 1) { x = montgomery_reduction(x * base, r, modulus, n_inv); } base = montgomery_reduction(base * base, r, modulus, n_inv); exponent >>= 1; } // Step 5 return montgomery_reduction(x, r, modulus, n_inv); } ``` 4. 中国剩余定理 ```c #include <stdio.h> #include <stdint.h> uint64_t chinese_remainder_theorem(uint64_t *c, uint64_t *n, uint64_t len) { // Step 1 uint64_t N = 1; for (int i = 0; i < len; i++) { N *= n[i]; } // Step 2 uint64_t x = 0; for (int i = 0; i < len; i++) { uint64_t Ni = N / n[i]; uint64_t Mi = mod_exp_montgomery(Ni, n[i] - 2, n[i]); x += c[i] * Ni * Mi % N; } // Step 3 return x % N; } ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值