Protocol LV10 算法的实现

clc;close all;clear all;
%R1 has SIZE strings
SIZE=10;
n=100;
range=2^n;


tag=1;


%生成10个在n比特范围内的十进制数
X=round(range*rand(1,SIZE));
% Y=X;
Y=round(range*rand(1,SIZE));
% X(SIZE/2:SIZE)=Y(SIZE/2:SIZE);
% 测试，能否accept，Y=X;
%从0到n平方的素数序列


TempPrimes=primes(n^2);


count=0;


while(tag)
    count=count+1;
    P=zeros(1,SIZE);
    P=randPrimes(SIZE,TempPrimes);
    S=mod(X,P);
    Q=mod(Y,P);
    J=S==Q;
    if(J==zeros(1,SIZE))
        disp('reject');
        disp('计算次数为： ');
        disp(count);
        tag=0;
    else
        for i=1:SIZE
            if(J(i)==1)
                if(X(i)==Y(i))
                    disp('accept');
                    disp('计算次数为： ');
                    disp(count);
                    tag=0;
                else
                    disp('?');
                    tag=1;
                end
                break;
            end
        end
    end
end

function [V] = randPrimes(Size,PV)
   V=zeros(1,Size);
   for i=1:Size
       temp = round((length(PV)-1)*rand(1,1))+1;
       V(i) = PV(temp); 
   end
end

概率算法课算法课的小作业。