x=randint(1,10)
pos=randint(1,2,[1,10])
x(1,pos)=~x(1,pos)
如果妳的碼流長度為N
pos=randint(1,2,[1,10])改成pos=randint(1,ceil(N*p),[1,N])即可