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
概率算法课算法课的小作业。