パパ活算が面白いというのは
月極払いの存在が計算を複雑にするからだな。

phimoseパパ活問題のRでの計算

d30=30
m80=69
f3=3
h1=1
t20=20
p2=-2


gr1=expand.grid(0:30,0:30,0:30)
gr2=cbind(gr1,apply(gr1,1,\(x) 30-sum(x)))
gr3=gr2[gr2[,4]>=0,]
colnames(gr3)=NULL
ans=gr3[apply(gr3,1,\(x) f3*x[1]+h1*x[2]+t20*as.numeric(x[3]>0))==m80,]
rownames(ans)=NULL
ans