2

我对密文和最后的正方形进行了一系列计算。问题是即使有足够的噪声预算来执行平方和重新线性化（在操作之前和之后），当我解密它时，我得到的结果不正确。

奇怪的是，如果我解密和解码，然后在一个新的密文中再次编码和加密，数字，在平方之前，计算是正确的。

我不明白问题出在哪里，因为如果设置的参数有误，那么无论如何我都会得到不正确的结果。

我指定我使用分数编码器，其中积分部分的多项式系数为 64，小数部分的精度为 32 位（以 3 为底），并且我在平方之前的所有计算都是在密文和明文之间执行的。

这是我的意思的一个例子：

#include <vector>
#include "seal/seal.h"

using namespace std;
using namespace seal;



int main(int argc, char const *argv[])
{
    //Set parameters
    EncryptionParameters parms;
    parms.set_poly_modulus("1x^4096 + 1");
    parms.set_coeff_modulus(coeff_modulus_128(4096));
    parms.set_plain_modulus(1<<20);
    SEALContext context(parms);
    KeyGenerator keygen(context);
    PublicKey public_key = keygen.public_key();
    SecretKey secret_key = keygen.secret_key();
    EvaluationKeys ev_keys16;
    keygen.generate_evaluation_keys(16, ev_keys16);
    Encryptor encryptor(context, public_key);
    Evaluator evaluator(context);
    Decryptor decryptor(context, secret_key);
    FractionalEncoder fraencoder(context.plain_modulus(), context.poly_modulus(), 64, 32, 3);




    float data=0.5;
    float weight= 1.5;

    //encrypt data and encode weight

    Ciphertext encrypted_data;
    encryptor.encrypt(fraencoder.encode(data),encrypted_data);
    cout<<"Noise budget in a freshly encrypted: "<<decryptor.invariant_noise_budget(encrypted_data)<<endl;
    Plaintext encoded_weight=fraencoder.encode(weight);

    //Operation: (0.5*1.5)*5=3.75
    vector<Ciphertext> mul_vector(5);
    for(int i=0;i<5;i++){
        evaluator.multiply_plain(encrypted_data,encoded_weight,mul_vector[i]);
    }

    evaluator.add_many(mul_vector,encrypted_data);
    cout<<"Noise budget after 5 plain multiplications and 4 additions: "<<decryptor.invariant_noise_budget(encrypted_data)<<endl;

    //Operation: 3.75*4=15
    vector<Ciphertext> add_vector(4);
    for(int i=0;i<4;i++){
        add_vector[i]= Ciphertext(encrypted_data);
    }
    evaluator.add_many(add_vector,encrypted_data);

    cout<<"Noise budget after 4 additions: "<<decryptor.invariant_noise_budget(encrypted_data)<<endl;

    //Operation: (15-1.5)*1.5=20.25
    evaluator.sub_plain(encrypted_data,encoded_weight);
    evaluator.multiply_plain(encrypted_data,encoded_weight);

    cout<<"Noise budget after 1 plain sub and 1 plain multiplication: "<<decryptor.invariant_noise_budget(encrypted_data)<<endl;

    //Operation: (20.25*1.5)*6=182.25
    vector<Ciphertext> mul_vector2(6);
    for(int i=0;i<6;i++){
        evaluator.multiply_plain(encrypted_data,encoded_weight,mul_vector2[i]);
    }

    evaluator.add_many(mul_vector2,encrypted_data);
    cout<<"Noise budget after 6 plain multiplications and 5 additions: "<<decryptor.invariant_noise_budget(encrypted_data)<<endl;

    // here I decrypt, decode and encrypt again, to obtain the right result
    Plaintext tmp;
    float res;
    //If I remove the following 4 lines the final result is incorrect
    decryptor.decrypt(encrypted_data,tmp);
    res = fraencoder.decode(tmp);
    cout<<"Decrypted result before square: "<<res<<endl;
    encryptor.encrypt(fraencoder.encode(res),encrypted_data);

    //182.25^2=33215.1
    evaluator.square(encrypted_data);
    evaluator.relinearize(encrypted_data,ev_keys16);

    cout<<"Noise budget after square and relianearization: "<<decryptor.invariant_noise_budget(encrypted_data)<<endl;


    decryptor.decrypt(encrypted_data,tmp);
    res= fraencoder.decode(tmp);
    cout<<res<<endl;


    return 0;
}

我错过了什么？

4

1 回答 1

3

问题是你plain_modulus对于这个计算来说太小了。

如果你tmp在最后打印出来 ( tmp.to_string()) 你会看到它有非常大的系数。请注意，这些系数是模数plain_modulus，因此预计它们中的许多看起来很大。尽管如此，res 作为整数的无穷范数（即，如果您使用足够大的plain_modulus）是 266441508，刚好低于 2 29。由于您plain_modulus只有 2 20，因此您最终会得到不正确的结果。您的重新编码会有所帮助，因为在最后一次计算之前，系数仍然不是太大，并且重新编码多项式会将系数降低到无穷范数 1（在您的情况下，base=3）。

plain_modulus解决方案是至少增加到2 29。当然，这会导致更多的噪音增长，并且您已经非常接近噪音溢出，但它仍然可以工作。

于 2018-08-28T22:40:31.430 回答