1

clf
m = 317;
a = 151;
r = ones(1,m-1);
for i=2:m-1
r(i) = mod(a*r(i-1), m);
end
N = 1:m-1;
k = zeros(1,m-1);
for i=1:m-1
k(i) = nnz(r==i);
end
plot(N,k,'.')
title('Nk vs k')
xlabel('k')
ylabel('N(k)')
fprintf(['Question 1: Note every value is 1, because it only appears once\n']);