aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--gp-scripts/m+1stPrice35
1 files changed, 17 insertions, 18 deletions
diff --git a/gp-scripts/m+1stPrice b/gp-scripts/m+1stPrice
index 8407e96..a6ac17c 100644
--- a/gp-scripts/m+1stPrice
+++ b/gp-scripts/m+1stPrice
@@ -5,6 +5,8 @@
5\\ Adapt the following values to your needs 5\\ Adapt the following values to your needs
6\\\\\\\\\\\\ 6\\\\\\\\\\\\
7 7
8\\ auction parameter
9M = 1
8\\ amount of bidders 10\\ amount of bidders
9n = 2^3 11n = 2^3
10\\ amount of possible prices 12\\ amount of possible prices
@@ -16,20 +18,18 @@ bid = vector(n,i,random(k)+1)
16\\bid = vector(n,i,(i+1)%2) \\ second bidder wins (with ties) 18\\bid = vector(n,i,(i+1)%2) \\ second bidder wins (with ties)
17 19
18\\ prime finite field setup (result may be ambiguous if your prime is too small, 4*n*k seems to work fine) 20\\ prime finite field setup (result may be ambiguous if your prime is too small, 4*n*k seems to work fine)
19q = prime(4*n*k) 21\\p = 263
20 22\\q = (p-1)/2
21\\\\\\\\\\\\ 23\\ use save prime p
22\\ SETUP 24q = prime(2^12)
23\\\\\\\\\\\\ 25p = 2*q + 1
24 26
25\\ p not needed? wat? 27\\ get generator / primitive element for G_q
26\\p = 47 28\\var = 'x \\ copy pasta from internet
27 29\\pe=ffgen(minpoly(ffprimroot(ffgen(ffinit(p,1))),var),var) \\ get primitive element
28\\ get generator / primitive element for Z_q 30\\1/(fforder(pe) == p-1) \\ error out, if ord(pe) is wrong
29var = 'x \\ copy pasta from internet 31\\g = Mod(eval(Str(pe))^2, p) \\ dirty hack to convert t_FFELEM to t_INT
30pe=ffgen(minpoly(ffprimroot(ffgen(ffinit(q,1))),var),var) \\ get primitive element 32g = Mod(4, p)
311/(fforder(pe) == q-1) \\ error out, if ord(pe) is wrong
32g = Mod(eval(Str(pe)), q) \\ dirty hack to convert t_FFELEM to t_INT
33 33
34\\\\\\\\\\\\ 34\\\\\\\\\\\\
35\\ PROLOG 35\\ PROLOG
@@ -66,8 +66,8 @@ beta = matrix(n,k,i,j, g^r[i,j])
66\\ multiplicative shares 66\\ multiplicative shares
67\\ first index level = owning agent id (multiplicative share) 67\\ first index level = owning agent id (multiplicative share)
68\\ second index level = agent id, price id 68\\ second index level = agent id, price id
69Gamma = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j+1,k,alpha[h,d])) * prod(d=1,j-1,alpha[i,d]) * prod(h=1,i-1,alpha[h,j]) )^m[a][i,j] )) 69Gamma = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j,k,alpha[h,d])*prod(d=j+1,k,alpha[h,d])) * prod(d=1,j,alpha[i,d])^(2*M+2) / g^(2*M+1) )^(m[a][i,j]) ))
70Delta = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j+1,k, beta[h,d])) * prod(d=1,j-1, beta[i,d]) * prod(h=1,i-1, beta[h,j]) )^m[a][i,j] )) 70Delta = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j,k, beta[h,d])*prod(d=j+1,k, beta[h,d])) * prod(d=1,j, beta[i,d])^(2*M+2) )^(m[a][i,j]) ))
71 71
72\\\\\\\\\\\\ 72\\\\\\\\\\\\
73\\ ROUND3 73\\ ROUND3
@@ -84,8 +84,7 @@ Phi = vector(n,a,matrix(n,k,i,j, prod(h=1,n,Delta[h][i,j])^x[a] ))
84 84
85\\ winner matrix 85\\ winner matrix
86v = matrix(n,k,a,j, prod(i=1,n,Gamma[i][a,j]) / prod(i=1,n,Phi[i][a,j]) ) 86v = matrix(n,k,a,j, prod(i=1,n,Gamma[i][a,j]) / prod(i=1,n,Phi[i][a,j]) )
87
88vi = lift(v) 87vi = lift(v)
89 88
90print("bids are: ", bid) 89print("bids are: ", bid)
91for(X=1,n, if(vecmin(vi[X,])==1, print("And the winner is ", X) )) 90for(X=1,n, if(vecmin(vi[X,],&i)==1, print("And the winner is ", X, " with the price ", i) ))