#include #include #define eb emplace_back #define all(a) (a).begin(),(a).end() #define RD(T,...) T __VA_ARGS__;li(__VA_ARGS__) #define LL(...) RD(ll,__VA_ARGS__) #define VE(a,...) au a=vec(__VA_ARGS__) #define JO(a,b) a##b #define jo(a,b) JO(a,b) #define uq_symbol(a) jo(a,__LINE__) #define FE(v,e,...) for(au&&__VA_OPT__([)e __VA_OPT__(,__VA_ARGS__]):v) #define FO(n) for(ll uq_symbol(i)=n;uq_symbol(i)-->0;) #define FOR(i,...) for(au[i,i##O,i##E]=rng(0,__VA_ARGS__);i=i##O;i-=i##E) #define I template #define J class #define au auto #define rr return using std::cin,std::cout,std::begin,std::end,std::rbegin; using std::swap,std::move,std::abs,std::prev,std::next; using std::tuple,std::array,std::bitset,std::minmax,std::get; using vo=void;using bo=bool; vo solve(); using is=std::istream;using os=std::ostream; using dd=long double;using ll=long long;using ull=unsigned long long; using lll=__int128_t;using ulll=__uint128_t;using str=std::string; using namespace atcoder;using ml=modint;au&operator<<(os&o,const ml&x){rr o<bo at2(T x,ll i){rr x>>i&1;} IT zz(T x){rr x<0?-x:x;} Ill len(const T&a){rr a.size();} Ibo amax(T&a,const T&b){rr astruct cp{ A a={};B b={}; cp(){} cp(A a,B b):a(a),b(b){} cp(std::pairp):a(p.first),b(p.second){} bo operator==(const cp&c)const{rr a==c.a&&b==c.b;} bo operator<(const cp&c)const{rr a!=c.a?a(const cp&c)const{rr a!=c.a?a>c.a:b>c.b;} friend is&operator>>(is&i,cp&c){rr i>>c.a>>c.b;} friend os&operator<<(os&o,const cp&c){rr o<; Istruct tr{ A a={};B b={};C c={}; tr(){} tr(A a,B b,C c):a(a),b(b),c(c){} bo operator==(const tr&t)const{rr a==t.a&&b==t.b&&c==t.c;} au operator<(const tr&t)const{rr a!=t.a?a>(is&i,tr&t){rr i>>t.a>>t.b>>t.c;} friend os&operator<<(os&o,const tr&t){rr o<; Idecltype(au)rb(T&a){rr*rbegin(a);} Iau pot(T&a){au r=a.top();a.pop();rr r;} Iusing map=std::map; Ios&operator<<(os&o,const map&m){fe(m,e)o<using umap=std::unordered_map; Ios&operator<<(os&o,const umap&m){fe(m,e)o<struct ve; Iconstexpr bo isv=0; Iconstexpr bo isv> =1; Iconstexpr bo isv> =1; Iau rawv(V){if constexpr(isv)rr rawv(V(1)[0]);else rr V();} Istruct ve:std::vector{ using std::vector::vector; using T=decltype(rawv(V())); Ive(const std::vector&v={}){fe(v,e)this->eb(e);} ve&operator+=(const ve&u){au&v=*this;fo(i,len(v))v[i]+=u[i];rr v;} ve&operator-=(const ve&u){au&v=*this;fo(i,len(v))v[i]-=u[i];rr v;} ve&operator^=(const ve&u){fe(u,e)this->eb(e);rr*this;} ve&operator+=(const T&x){au&v=*this;fe(v,e)e+=x;rr v;} ve&operator-=(const T&x){au&v=*this;fe(v,e)e-=x;rr v;} ve&operator*=(const T&x){au&v=*this;fe(v,e)e*=x;rr v;} Iau&operator+=(const bitset&a){fo(i,n)(*this)[i]+=a[i];rr*this;} Iau&operator-=(const bitset&a){fo(i,n)(*this)[i]-=a[i];rr*this;} ve operator+(const ve&u)const{rr ve(*this)+=u;} ve operator-(const ve&u)const{rr ve(*this)-=u;} ve operator^(const ve&u)const{rr ve(*this)^=u;} ve operator+(const T&x)const{rr ve(*this)+=x;} ve operator-(const T&x)const{rr ve(*this)-=x;} ve operator*(const T&x)const{rr ve(*this)*=x;} ve&operator++(){rr*this+=1;} ve&operator--(){rr*this-=1;} ve operator-()const{rr ve(*this)*=-1;} au lower_bound(const V&x)const{rr std::lower_bound(all(*this),x);} au upper_bound(const V&x)const{rr std::upper_bound(all(*this),x);} }; Iau vec(const ll(&s)[n],T x={}){ if constexpr(n==i+1)rr ve(s[i],x); else{au X=vec(s,x);rr ve(s[i],X);} } Ivo setsize(const ll(&l)[n],A&...a){((a=vec(l,rawv(a))),...);} Iusing vve=ve>;using vl=ve;using vvl=vve; using gl=ve>;using ge=ve>; struct io{io(){cin.tie(0)->sync_with_stdio(0); cout<os&operator<<(os&o,const bitset&b){fo(i,n)o<os&operator<<(os&o,const tuple&t){ apply([&](const au&...a){ll i=0;(((o<os&operator<<(os&o,const std::vector&v){fe(v,e)o<?nl:sp);rr o;} Ivo pp(const A&...a){ll i=0;((cout<is&operator>>(is&i,std::vector&v){fe(v,e)i>>e;rr i;} Ivo li(A&...a){(cin>>...>>a);} Ivo vi(A&...a){fo(i,(len(a)&...))(cin>>...>>a[i]);} Istruct rec:F{rec(F&&f):F(std::forward(f)){}Idecltype(au)operator()(A&&...a)const{rr F::operator()(*this,std::forward(a)...);}}; int main(){ll T=1;fo(T)solve();} vo solve(){ LL(N,M); ve
a(N),b(N),p(N);vi(a,b,p); VE(dp,{pw2(N),M+1},dd(-1)); VE(used,{pw2(N),M+1}); au f=rec([&](au&f,ll S,ll k)->dd{ if(dp[S][k]!=-1)rr dp[S][k]; dd res=0; fo(i,N){ if(!at2(S,i)){ amax(res,f(S|pw2(i),k)+1/a[i]); if(k)amax(res,1/p[i]*(f(S|pw2(i),k)+1/b[i])+(1-1/p[i])*f(S,k-1)); } } rr dp[S][k]=res; }); pp(f(0,M)); }