// URL: https://yukicoder.me/problems/no/766 import std.algorithm, std.array, std.container, std.math, std.range, std.typecons, std.string; version(unittest) {} else void main() { int N, M, P; io.getV(N, M, P); int[] V; io.getA(N, V); V.sort!"a>b"; auto w = new mint[](N); foreach (i; 0..N) w[i] = V[i]; auto wc = w.cumulativeSum; auto pp = new mint[](M+1); pp[0] = 1; foreach (i; 1..M+1) pp[i] = pp[i-1] * P / 100; auto p1p = new mint[](N+1); p1p[0] = 1; foreach (i; 1..N+1) p1p[i] = p1p[i-1] * (100-P) / 100; auto fact = Fact!mint(N+M); auto r = mint(0); foreach (i; 0..M) r += wc[0..N] * fact.combi(N-1+i, i) * p1p[N] * pp[i]; foreach (i; 0..N) r += wc[0..i] * fact.combi(M-1+i, i) * p1p[i] * pp[M]; io.put(r); } class CumulativeSum(T) { const size_t n; T[] s; this(T[] a) { n = a.length; s = new T[](n+1); s[0] = 0; foreach (i; 0..n) s[i+1] = s[i]+a[i]; } pure T opSlice(size_t l, size_t r) { return s[r]-s[l]; } pure size_t opDollar() { return n; } } CumulativeSum!T cumulativeSum(T)(T[] a) { return new CumulativeSum!T(a); } struct Fact(T, bool inv = true) { int n; T[] table; alias table this; static if (inv) T[] invTable; this(int n) { this.n = n; table = new T[](n+1); table[0] = 1; foreach (i; 1..n+1) table[i] = table[i-1]*i; static if (inv) { invTable = new T[](n+1); invTable[n] = T(1)/table[n]; foreach_reverse (i; 1..n+1) invTable[i-1] = invTable[i]*i; } } static if (inv) { pure T perm(size_t a, size_t b) in { assert(n >= a && a >= b); } do { return table[a]*invTable[a-b]; } pure T combi(size_t a, size_t b) in { assert(n >= a && a >= b); } do { return table[a]*invTable[b]*invTable[a-b]; } pure T homo(size_t a, size_t b) in { assert(n >= a+b-1); } do { return combi(a+b-1, b); } } } const mod = 10^^9+7; alias mint = ModInt!mod; auto lowerBoundBy(alias conv = "a", alias comp = "a tuple(e, convFun(e))).assumeSorted!((a, b) => compFun(a[1], b[1])); } pure T isqrt(T)(T n) { static if (is(T == int)) auto max = 46341; else static if (is(T == long)) auto max = 3037000500L; return iota(T(0), max).lowerBoundBy!("a^^2", "a<=b")(n).back; } pure T icbrt(T)(T n) { static if (is(T == int)) auto max = 1291; else static if (is(T == long)) auto max = 2097152L; return iota(T(0), max).lowerBoundBy!("a^^3", "a<=b")(n).back; } pure T powr(alias pred = "a*b", T, U)(T a, U n, T one) { import std.functional; alias predFun = binaryFun!pred; if (n == 0) return one; auto r = one; for (; n > 0; n >>= 1) { if (n&1) r = predFun(r, a); a = predFun(a, a); } return r; } pure T powr(alias pred = "a*b", T, U)(T a, U n) { return powr!(pred, T, U)(a, n, T(1)); } pure T extGcd(T)(T a, T b, out T x, out T y) { auto g = a; x = 1; y = 0; if (b) { g = extGcd(b, a%b, y, x); y -= a/b*x; } return g; } struct ModInt(int m, bool pos = false) { version(BigEndian) union { long l; struct { int i2; int i; } } else union { long l; int i; } alias M = ModInt!(m, pos), Op = string; @property static init() { return M(0); } @property int value() { return i; } @property void value(int v) { i = nm(v); } alias value this; this(T)(T v = 0) { i = nm(v); } ref M opAssign(T)(T v) { i = nm(v); return this; } static if (!pos) pure M opUnary(Op op: "-")() { return M(-i); } static if (m < int.max / 2) { pure M opBinary(Op op)(int r) if(op=="+"||!pos&&op=="-") { return M(mixin("i"~op~"r")); } ref M opOpAssign(Op op)(int r) if(op=="+"||!pos&&op=="-") { i=nm(mixin("i"~op~"r")); return this; } } else { pure M opBinary(Op op)(int r) if(op=="+"||!pos&&op=="-") { return M(mixin("l"~op~"r")); } ref M opOpAssign(Op op)(int r) if(op=="+"||!pos&&op=="-") { i=nm(mixin("l"~op~"r")); return this; } } pure M opBinary(Op op: "*")(int r) { return M(l*r); } ref M opOpAssign(Op op: "*")(int r) { i=nm(l*r); return this; } pure M opBinary(Op op)(M r) if(op=="+"||!pos&&op=="-"||op=="*") { return opBinary!op(r.i); } ref M opOpAssign(Op op)(M r) if(op=="+"||!pos&&op=="-"||op=="*") { return opOpAssign!op(r.i); } pure M opBinary(Op op: "^^", T)(T n) { return this.powr(n); } static if (!pos) { pure M opBinary(Op op: "/")(M r) { return M(l*r.inv.i); } ref M opOpAssign(Op op: "/")(M r) { i=nm(l*r.inv.i); return this; } pure M opBinary(Op op: "/")(int r) { return opBinary!op(M(r)); } ref M opOpAssign(Op op: "/")(int r) { return opOpAssign!op(M(r)); } pure M inv() { int x = i, a, b; extGcd(x, m, a, b); return M(a); } } private { pure int nm(int v) { static if (pos) return v%m; else return (v%m+m)%m; } pure int nm(long v) { static if (pos) return cast(int)(v%m); else return cast(int)((v%m+m)%m); } } } auto io = IO!()(); import std.stdio; struct IO(alias IN = stdin, alias OUT = stdout) { import std.conv, std.format, std.meta, std.traits, core.stdc.stdlib; auto getV(T...)(ref T v) { foreach (ref w; v) get(w); } auto getA(T)(size_t n, ref T v) if (hasAssignableElements!T) { v = new T(n); foreach (ref w; v) get(w); } auto getC(T...)(size_t n, ref T v) if (allSatisfy!(hasAssignableElements, T)) { foreach (ref w; v) w = new typeof(w)(n); foreach (i; 0..n) foreach (ref w; v) get(w[i]); } auto getM(T)(size_t r, size_t c, ref T v) if (hasAssignableElements!T && hasAssignableElements!(ElementType!T)) { v = new T(r); foreach (ref w; v) getA(c, w); } template getS(E...) { auto getS(T)(size_t n, ref T v) { v = new T(n); foreach (ref w; v) foreach (e; E) mixin("get(w."~e~");"); } } const struct PutConf { bool newline = true, flush, exit; string floatFormat = "%.10f", delimiter = " "; } auto put(alias conf = "{}", T...)(T v) { mixin("const PutConf c = "~conf~"; putMain!c(v);"); } auto putB(alias conf = "{}", S, T)(bool c, S t, T f) { if (c) put(t); else put(f); } auto putRaw(T...)(T v) { OUT.write(v); OUT.writeln; } private { dchar[] buf; auto sp = (new dchar[](0)).splitter; void nextLine() { IN.readln(buf); sp = buf.splitter; } auto get(T)(ref T v) { if (sp.empty) nextLine(); v = sp.front.to!T; sp.popFront(); } auto putMain(PutConf c, T...)(T v) { foreach (i, w; v) { putOne!c(w); if (i < v.length-1) OUT.write(c.delimiter); } static if (c.newline) OUT.writeln; static if (c.flush) OUT.flush(); static if (c.exit) exit(0); } auto putOne(PutConf c, T)(T v) { static if (isInputRange!T && !isSomeString!T) putRange!c(v); else if (isFloatingPoint!T) OUT.write(format(c.floatFormat, v)); else OUT.write(v); } auto putRange(PutConf c, T)(T v) { auto w = v; while (!w.empty) { putOne!c(w.front); w.popFront(); if (!w.empty) OUT.write(c.delimiter); } } } }