Implementar a reconstrução do compartilhamento secreto de Shamir


11

O esquema de compartilhamento secreto de Shamir é uma maneira simples de proteger um segredo, dividindo-o em várias partes necessárias para reconstruí-lo.

Sua tarefa é implementar a reconstrução do Compartilhamento Secreto de Shamir sobre o Campo Finito definido pelos primos 1928049029. Se você tiver alguma dúvida sobre o que isso significa, basta perguntar ou ver Campo Finito e Aritmética de Campos Finitos na wikipedia (mais recursos abaixo).

Entrada

A entrada é feita usando stdin. Primeiro vem um número inteiro k, depois as linhas k. Cada uma dessas linhas contém um par de números inteiros x yque representam um segredo. Em outras palavras, f(x) = yno polinômio original usado para construir os segredos.

O número de segredos dados é sempre suficiente para construir o segredo correspondente.

Resultado

Saída para mostrar o segredo reconstruído.

Exemplos

Entrada:

5         
1 564797566
2 804114535
4 1354242660
6 1818201132
7 503769263

Resultado:

1234

Entrada:

7
1 819016192
2 1888749673
3 1737609270
4 365594983
5 1628804870
6 1671140873
7 492602992

Resultado:

456457856

Recursos

Artigo da Wikipedia

Papel

Campo finito Fonte: Wikipedia

Aritmética de campos finitos Fonte: Wikipedia

Polinômio de Lagrange Fonte: Wikipedia

Capítulo sobre aritmética de campos finitos

Respostas:


4

bash, 271 caracteres

r () {
[$ {1/0 /}] && {r $ (($ 2% $ 1)) $ 1; ((t = u, u = v- $ 2 / $ 1 * u, v = t));}
}
ler
((N = 1928049029, n = 0))
enquanto lê x [$ n] y [$ n]
do ((n ++))
feito
para ((i = n; z = (z + l)% N, i -;)) faça
para ((j = n, l = y [i]; j -;)) faça
((u = 0, v = 1, d = x [j] -x [i], M = N + d))
r MN
[$ {d / 0 /}] && ((l = l * x [j]% N * (u + N)% N))
feito
feito
eco $ z

As novas linhas podem ser substituídas na maioria dos casos por ponto e vírgula, mas não acho que exista espaço em branco desnecessário.

(Eu não tinha percebido antes hoje que os números inteiros do bash são de 64 bits - muito útil).

Para o bash, o GCD recursivo (explorando o estado global) parece ser mais compacto do que o iterativo. Isso é principalmente direto; o truque interessante é o [ ${d/0/} ]&&fooque é efetivamenteif [ $d -ne 0 ];then foo;fi


Agradável! Eu nunca esperava ver uma resposta básica para esse problema. 1
Juan

@ Juan, comecei a fazê-lo em Perl e me cansei de ter que forçá-lo a fazer divisão inteira em vez de flutuar. E eu sei o bash melhor de qualquer maneira, por isso envolve menos batidas da cabeça contra a parede.
Peter Taylor

3

199 caracteres no Octave:

m=@(x)mod(x,1928049029);[d,l]=scanf('%d');c=d(1);e=repmat(int64(d(2:2:l)),1,c);[_,b]=gcd(e-e',1928049029*ones(c));b=eye(c)+m(e.*b);x=b(1,:);for i=2:c;x=m(x.*b(i,:));end;disp(m(sum(m(x'.*d(3:2:l)))))

3

Golfscript, 114 112 111 110 109 65 (86) caracteres

Se você não se importa em obter resultados esta semana, basta 65 caracteres:

~](;2/0\:X{~\.X{0=}%^\{\.@- 1928049029:P.,\@{@*\%(!}++?**}+/+P%}/

Mas se você está procurando eficiência, é um pouco mais longo com 86 caracteres:

~](;2/0\:X{~\.X{0=}%^\{\[.0](@-[1928049029:P%P]{.~/{\.(;@@~@*-+\}+2*.1=}do;0=*}+/+P%}/

Isso é dissecado com muito mais detalhes do que desejo repetir aqui no meu blog .


Principalmente não é o meu trabalho, mas o uso pesado do Nabb fornece 47 caracteres:

n%(!\:A{~A{~;.3$- 1928049029:N((?1or**}/\/+N%}/

Nota: Eu apenas raciocinei sobre esse código: tentar executá-lo seria inútil, dado o período de tempo e a quantidade de memória que ele usaria.


3

Golfe - 52 46 (67)

Uma abordagem de força bruta para inversos modulares em 46 caracteres. Calcula repetidamente um ^ (N-2) com números inteiros de precisão arbitrários.

n%(!\:A{~A{~;.3$-.!+1928049029:N((?**}/\/+N%}/

A implementação do algoritmo euclidiano estendido custa apenas 15 caracteres adicionais.

n%(!\:A{~A{~;.3$-:|!1\1928049029:N{@2$|3$/*-\|\:|%.}do;;**}/\/+N%}/

Esse código é totalmente detalhado no meu blog , incluindo algumas alternativas para calcular a inversa multiplicativa modular.


11
Bom, mas acho que ainda há pelo menos dois caracteres a serem salvos. Substitua {*N%2<}por {*N%1=}como no blog e você pode abandonar o (;depois N,. Mas então, para a entrada de desempenho é irrelevante, você pode usar o pequeno teorema de Fermat sem se preocupar com o lado modular da exponenciação - apenas deixe isso para a arrumação final - para que o destinatário se torne N((?.
Peter Taylor #

11
@ Peter: {*N%1=}+vai faltar o caso com o denominador zero, que levaria pelo menos 3 caracteres para lidar. Porém, é bom pegar simplesmente x ^ (N-2), e podemos obter 46 caracteres usando isso.
Nabb

2

Lua 444 Chars

Funciona para o exemplo na página wiki

3
2 1942
4 3402
5 4414

Mas de alguma forma não funciona para os exemplos aqui nesta página. Se alguém puder encontrar o erro?

Versão sem golfe:

-- Reconstruct shamir secret
-- convention, poly = {[0]=a0,a1,...,an}
i=io.read
f=math.fmod
w=1928049029
k=i():match"%d+"
x={} -- Will contain X values
y={} -- Will contain Y values
p={} -- will contain lagrange polynomials

-- Read data
for j=0,k-1 do
    x[j],y[j]=i():match("(%d+) (%d+)")
    print(j,x[j],y[j])
end
-- Multiplication and scaling function
function mul(p,q,s)
    -- multiply polies
    r={} -- poly to be returned
    for k=0,#p do 
        for l=0,#q do
            r[l+k]=r[l+k] or 0 -- if the coeff for degree l+k of x doesn't exist, put 0
            p[k]=p[k] or 0 -- if p hasn't got a coeff for x^k
            q[l]=q[l] or 0 -- idem for q
            r[l+k]=(r[l+k]+s*p[k]*q[l]%w -- calculate increment for coeff for x^(l+k) 
        end
    end
    -- Debugging
    io.write"Multiplied "
    printPoly(p)
    io.write"With       "
    printPoly(q)
    io.write("And scaling factor ",tostring(s),"\n")
    io.write"Yielding   "
    printPoly(r)
    return r
end

function printPoly(p) -- "Pretty" printing of the polynomial
    for k=#p,1,-1 do
        io.write(tostring(p[k] or 0),"x^",tostring(k),"+")
    end
    io.write(p[0])
    io.write"\n"
end
function egcd(a,b)
    if a == 0 then
        return b, 0, 1
    else
        local g, y, x = egcd(b % a, a)
        return g, x - math.floor(b / a) * y, y
    end
end

function inv(a,m)
    a=a>=0 and a or a+m
    local g,x,y = egcd(a,m)
    if g== 1 then
        return x%m
    else
        print(a,"has no inverse mod",m)
    end
end


-- generate lagrange polynomials
for j=0,#x do
    print("j=",j,"*********")
    for m=0,k-1 do
        if m~=j then -- if m==j, continue
            p[j]=p[j]or{[0]=1} -- if this poly doesn't exist, take 1
            p[j]=mul( p[j], {[0]=-x[m],1},inv(x[j]-x[m],w))-- multiply with (x-x_m)/(x_j-x_m)
            io.write"---------------------------------\n"
        end
    end
end
r=0 -- Result for x^0
for k=0,#p do
    print("l_"..k)
    printPoly(p[k]) -- print l_k
    r=r+f(y[k]*p[k][0],w) -- add coeff for x^0 to result
end
print("Secret was",f(r,w)) -- display result

Golfe (sem campo finito), 444 caracteres:

i=io.read f=math.fmod w=1928049029 k=i():match"%d+"x={}y={}p={}for j=0,k-1 do x[j],y[j]=i():match("(%d+) (%d+)")end
function mul(p,q,s)r={}for k=0,#p do for l=0,#q do r[l+k]=r[l+k]or 0 p[k]=p[k]or 0 q[l]=q[l]or 0 r[l+k]=f(r[l+k]+s*p[k]*q[l],w)end end return r end
for j=0,#x do for m=0,k-1 do if m~=j then p[j]=p[j]or{[0]=1}p[j]=mul(p[j],{[0]=-x[m],1},1/(x[j]-x[m]))end end end r=0 for k=0,#p do r=r+f(y[k]*p[k][0],w)end
print(f(r,w))

O exemplo da Wikipedia não usa um campo finito, o que é realmente uma pena, que teria sido muito mais instrutivo. Essa é provavelmente a fonte do seu erro.
Aaaaaaaaaaaa

2

Java, 435407 caracteres

import java.util.*;public class G{public static void main(String[]args){Scanner s=new Scanner(System.in);int i,k,n=s.nextInt();long N=1928049029L,x[]=new long[n],y[]=new long[n],z=0,l,c;for(i=n;i-->0;){x[i]=s.nextInt();y[i]=s.nextInt();}for(i=n;i-->0;){l=y[i];for(long j:x)if(x[i]!=j){c=1;for(long a=N+j-x[i],b=N,d=0,t;b>0;){t=d;d=c-a/b*d;c=t;t=b;b=a%b;a=t;}l=l*j%N*(c+N)%N;}z+=l;}System.out.println(z%N);}}

Ungolfed:

import java.util.*;
public class G {
    public static void main(String[] args) {
        Scanner s=new Scanner(System.in);
        int i,k,n=s.nextInt();
        long N=1928049029L,x[]=new long[n],y[]=new long[n],z=0,l,c;
        for (i=n; i-->0;) {
            x[i]=s.nextInt();
            y[i]=s.nextInt();
        }
        for (i=n; i-->0;) {
            l=y[i];
            for (long j:x)
                if (x[i]!=j) {
                    // Extended Euclid algorithm - iterative version -
                    // to find the reciprocal of j-x[i] (mod N)
                    c=1;
                    for (long a=N+j-x[i], b=N, d=0, t; b>0;) {
                        t=d; d=c-a/b*d; c=t;
                        t=b; b=a%b; a=t;
                    }
                    l = l*j%N;
                    l = l*(c+N)%N;
                }
                z+=l;
        }
        System.out.println(z%N);
    }
}

2

Haskell, 183

p=1928049029
a#0=(1,0)
a#b=let(s,t)=b#mod a b in(t,s-div a b*t)
s d=sum[y*product[z*fst((z-x)#p)|[z,_]<-d,z/=x]|[x,y]<-d]
main=interact$show.(`mod`p).s.map(map read.words).tail.lines
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.