3 * $Id: share.c,v 1.4 2000/10/08 12:16:17 mdw Exp $
5 * Shamir's secret sharing
7 * (c) 2000 Straylight/Edgeware
10 /*----- Licensing notice --------------------------------------------------*
12 * This file is part of Catacomb.
14 * Catacomb is free software; you can redistribute it and/or modify
15 * it under the terms of the GNU Library General Public License as
16 * published by the Free Software Foundation; either version 2 of the
17 * License, or (at your option) any later version.
19 * Catacomb is distributed in the hope that it will be useful,
20 * but WITHOUT ANY WARRANTY; without even the implied warranty of
21 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
22 * GNU Library General Public License for more details.
24 * You should have received a copy of the GNU Library General Public
25 * License along with Catacomb; if not, write to the Free
26 * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
30 /*----- Revision history --------------------------------------------------*
33 * Revision 1.4 2000/10/08 12:16:17 mdw
34 * Use @MP_EQ@ instead of @MP_CMP@.
36 * Revision 1.3 2000/06/24 18:29:05 mdw
37 * Interface change: allow shares to be extracted from a context on demand,
38 * rather than building them all up-front.
40 * Revision 1.2 2000/06/18 23:05:19 mdw
41 * Minor performance tweak: use Barrett reduction rather than Montgomery.
42 * Fast secret sharing isn't done here, though: see `gfshare' instead.
44 * Revision 1.1 2000/06/17 12:09:38 mdw
45 * Shamir's secret sharing system.
49 /*----- Header files ------------------------------------------------------*/
58 #include "mpbarrett.h"
64 /*----- Main code ---------------------------------------------------------*/
66 /* --- @share_create@ --- *
68 * Arguments: @share *s@ = pointer to share context to initialize
69 * @unsigned t@ = threshold for the system
73 * Use: Initializes a sharing context.
76 void share_create(share
*s
, unsigned t
)
85 /* --- @share_destroy@ --- *
87 * Arguments: @share *s@ = pointer to share context to destroy
91 * Use: Disposes of a sharing context. All memory is freed, all
92 * integers are dropped.
95 void share_destroy(share
*s
)
99 /* --- Dispose of the share vector --- */
102 for (i
= 0; i
< s
->t
; i
++) {
109 /* --- Other stuff --- */
117 /* --- @share_mkshares@ --- *
119 * Arguments: @share *s@ = pointer to share context to fill in
120 * @grand *r@ = pointer to random number source
124 * Use: Initializes a sharing context to be able to create shares.
125 * The context structure is expected to be mostly filled in. In
126 * particular, @t@ and @s@ must be initialized. If @p@ is zero,
127 * a prime number of appropriate size is generated
128 * automatically. If @v@ is zero, a vector of appropriate size
129 * is allocated. You should use the macro @SHARE_INIT@ or
130 * @share_create@ to construct sharing contexts.
133 void share_mkshares(share
*s
, grand
*r
)
137 /* --- If there's no prime, construct one --- */
143 unsigned bits
= (mp_octets(s
->s
) + 1) * 8;
146 p
= mprand(MP_NEW
, bits
, r
, 1);
147 s
->p
= pgen("p", p
, p
, 0, 0, 0, pgen_filter
, &pf
,
148 rabin_iters(bits
), pgen_test
, &pr
);
151 /* --- Construct the polynomial --- */
154 s
->v
= xmalloc(s
->t
* sizeof(share_pt
));
155 for (i
= 0; i
< s
->t
- 1; i
++)
156 s
->v
[i
].y
= mprand_range(MP_NEWSEC
, s
->p
, r
, 0);
157 s
->v
[s
->t
- 1].y
= mp_copy(s
->s
);
160 /* --- @share_get@ --- *
162 * Arguments: @share *s@ = pointer to share conext
163 * @mp *d@ = destination for the share
164 * @unsigned x@ = share index to fetch
166 * Returns: The share, as requested.
168 * Use: Extracts a share from the system. You may extract @MPW_MAX@
169 * shares, or @s->p@ shares from the system, whichever is
170 * smaller. Shares are indexed from 0.
173 mp
*share_get(share
*s
, mp
*d
, unsigned x
)
180 /* --- Various bits of initialization --- */
182 mp_build(&u
, &uw
, &uw
+ 1);
186 /* --- Evaluate the polynomial at %$x = i + 1$% --- */
189 mpbarrett_create(&mb
, s
->p
);
190 for (i
= 0; i
< s
->t
; i
++) {
191 d
= mp_mul(d
, d
, &u
);
192 d
= mp_add(d
, d
, s
->v
[i
].y
);
193 d
= mpbarrett_reduce(&mb
, d
, d
);
195 mpbarrett_destroy(&mb
);
200 /* --- @share_add@ --- *
202 * Arguments: @share *s@ = pointer to sharing context
203 * @unsigned x@ = which share number this is
204 * @mp *y@ = the share value
206 * Returns: Number of shares required before recovery may be performed.
208 * Use: Adds a share to the context. The context must have been
209 * initialized with the correct prime @p@ and threshold @t@.
212 unsigned share_add(share
*s
, unsigned x
, mp
*y
)
214 /* --- If no vector has been allocated, create one --- */
218 s
->v
= xmalloc(s
->t
* sizeof(share_pt
));
220 for (i
= 0; i
< s
->t
; i
++)
224 assert(((void)"Share context is full", s
->i
< s
->t
));
226 /* --- Store the share in the vector --- */
228 s
->v
[s
->i
].x
= x
+ 1;
229 s
->v
[s
->i
].y
= mp_copy(y
);
234 return (s
->t
- s
->i
);
237 /* --- @share_combine@ --- *
239 * Arguments: @share *s@ = pointer to share context
241 * Returns: The secret, as a multiprecision integer.
243 * Use: Reconstructs a secret, given enough shares.
246 mp
*share_combine(share
*s
)
255 /* --- Sanity checking --- */
257 assert(((void)"Not enough shares yet", s
->i
== s
->t
));
259 /* --- Initialization --- */
261 mpbarrett_create(&mb
, s
->p
);
262 mp_build(&ii
, &iiw
, &iiw
+ 1);
263 mp_build(&jj
, &jjw
, &jjw
+ 1);
265 /* --- Grind through the shares --- */
267 for (i
= 0; i
< s
->t
; i
++) {
271 for (j
= 0; j
< s
->t
; j
++) {
275 if (s
->v
[j
].x
>= s
->v
[i
].x
)
276 m
= mp_sub(m
, &jj
, &ii
);
278 m
= mp_sub(m
, &ii
, &jj
);
279 m
= mp_sub(m
, s
->p
, m
);
281 mp_gcd(0, 0, &m
, s
->p
, m
);
282 c
= mp_mul(c
, c
, &jj
);
283 c
= mpbarrett_reduce(&mb
, c
, c
);
285 c
= mpbarrett_reduce(&mb
, c
, c
);
287 c
= mp_mul(c
, c
, s
->v
[i
].y
);
288 c
= mpbarrett_reduce(&mb
, c
, c
);
293 a
= mpbarrett_reduce(&mb
, a
, a
);
297 mpbarrett_destroy(&mb
);
301 /*----- Test rig ----------------------------------------------------------*/
307 static int verify(grand
*r
)
309 unsigned n
= r
->ops
->range(r
, 16) + 8;
310 unsigned t
= r
->ops
->range(r
, n
- 1) + 1;
311 unsigned len
= r
->ops
->range(r
, 160);
313 mp
**v
= xmalloc(t
* sizeof(mp
*));
314 unsigned *p
= xmalloc(n
* sizeof(unsigned));
315 mp
*sec
= mprand(MP_NEW
, len
, r
, 0);
323 for (i
= 0; i
< n
; i
++)
325 for (i
= 0; i
< t
; i
++) {
326 unsigned long j
= r
->ops
->range(r
, n
- i
);
334 share_mkshares(&s
, r
);
335 for (i
= 0; i
< t
; i
++)
336 v
[i
] = share_get(&s
, MP_NEW
, p
[i
]);
339 assert(mparena_count(MPARENA_GLOBAL
) + mparena_count(MPARENA_SECURE
) == t
+ 2);
343 for (i
= 0; i
< t
; i
++)
344 share_add(&s
, p
[i
], v
[i
]);
345 ss
= share_combine(&s
);
348 if (!MP_EQ(sec
, ss
)) {
350 fprintf(stderr
, "\nbad recombination of shares\n");
356 for (i
= 0; i
< t
; i
++)
362 assert(mparena_count(MPARENA_GLOBAL
) + mparena_count(MPARENA_SECURE
) == 0);
368 grand
*r
= fibrand_create(0);
372 fputs("share: ", stdout
);
373 for (i
= 0; i
< 40; i
++) {
381 fputs(" ok\n", stdout
);
383 fputs(" failed\n", stdout
);
384 return (ok ? EXIT_SUCCESS
: EXIT_FAILURE
);
389 /*----- That's all, folks -------------------------------------------------*/