Let [5] be the set of the first 5 positive integers. We let α_=(αA)A≠∅,A⊆[5] consist of a vector with 31 real entries, where each αA is associated with a subset A⊆[5].
Define OBJ(α_)=∑A⊆[5],A≠∅αAlog(|A|), v(α_)=∑A⊆[5],A≠∅αA, and E(α_)=∑A,B:A∩B≠∅αAαB, where the sum for E(α_) is taken over all unordered pairs of disjoint nonempty sets A and B.
Also define FEAS(1/4) to be the set of all such vectors α_ with nonnegative real entries such that v(α_)=1 and E(α_)≥1/4.
I want to learn how to program the following optimization problem: OPT(1/4):=max
I was told that I can do this in SageMath. I have some basic knowledge of how to use Sage. How could I create the set FEAS(1/4)? I think that from there I may be able to figure out how to maximize OBJ(\underline{\alpha}) over this set.