[问题] 一个算法的问题

楼主: bravooooo (bravooooo)   2022-05-16 09:44:42
开发平台(Platform): (Ex: Win10, Linux, ...)
ubuntu 12.x
编译器(Ex: GCC, clang, VC++...)+目标环境(跟开发平台不同的话需列出)
g++
额外使用到的函数库(Library Used): (Ex: OpenGL, ...)
问题(Question):
大家好 最近遇到一个算法的问题
题目是要求解 a1/(1+a1*b1*s) + a2/(1+a2*b2*s)... + an/(1+an*bn*s)
加总后的各个系数
以2项为例就是:
a1/(1+a1*b1*s) + a2/(1+a2*b2*s) =
分母通分 → [a1(1+a2*b2) + a2(1+a1*b1)] *s / (1+a1*b1*s)(1+a2*b2*s) =
[a1(1+a2*b2) + a2(1+a1*b1)] *s / [1 + (a1*b1+a2*b2)s + (a1*b1*a2*b2)s^2]
上式中黄色部份即为所求
题目的多项式用暴力展开可以得到分母部份可以用组合公式把每一项求出
但若 n的值太大组合数会超多 计算量也跟着爆炸
不知道像这种题目有没有什么更好的解法呢
述叙的有点乱请大家见谅 也谢谢大家帮忙 Q_Q
喂入的资料(Input):
预期的正确结果(Expected Output):
错误结果(Wrong Output):
程式码(Code):(请善用置底文网页, 记得排版,禁止使用图档)
补充说明(Supplement):
作者: j0958322080 (Tidus)   2022-05-16 10:15:00
可以去数学版问问
作者: closer76 (克楼瑟)   2022-05-16 11:50:00
你的 (a1, a2,...,an) 和 (b1,b2,...bn) 是有规则的吗?
楼主: bravooooo (bravooooo)   2022-05-16 13:02:00
没有规则的哦 QQ
作者: hongsiangfu   2022-05-16 14:39:00
s domain model? 只能找通分后的系数规则了
作者: jack7775kimo (阿庞)   2022-05-16 14:54:00
Iteration的作法应该就是你说的暴力法? ((1,2)处理后跟3) Divide and Conquer有用上了吗? 会变O(logN)
作者: TheOneisNEO (Thomas Anderson)   2022-05-16 15:17:00
bn是不是没有单独用到...?可以先把an*bn算一下?
作者: jack7775kimo (阿庞)   2022-05-16 16:04:00
更:D&C应无助益,因为s的最高次方在不同层的时候不同
作者: TheOneisNEO (Thomas Anderson)   2022-05-16 16:44:00
应该是没特别办法 没有规则又每一个都要算出来大概就是把算过的先存起来 或考虑一些恒等式
作者: hellophoenix (Rainey)   2022-05-16 17:39:00
能不能看看n=3的长相,也许很规律
作者: mikemike1021 (mike)   2022-05-16 18:18:00
直接算会太多吗?an, bn 如果都整数用阵列存s多项式分母部分直接乘或d&c来做大数乘法 假设结果为P ,分子只需要在利用大数除法算出 sum ai * P/(1+aibis) 这样?
作者: oToToT (屁孩)   2022-05-16 19:54:00
推荐个Prob_Solve板然后分治应该做得到T(n)=2T(n/2)+O(多项式乘法)的复杂度吧如果多项式乘法naive套个fft的话应该可以O(n log^2 n)?(分母部分)

Links booklink

Contact Us: admin [ a t ] ucptt.com