一只懒猫的博客

原题链接poj 1273问题图论中有一个重要问题，假设有一个水管网图，每根水管都规定了一个可以流的方向，从v号结点到u号结点最大允许通过的水量为c(v,u)，有一个源点s存储了无限多的水，有一个汇点可以存储无限多的水。问：从源点流向汇点的水流的最大速率是多少？思路采用经典的EK算法，时间复杂度为O(n*m^2)n为结点数，m为边数，虽然上界很高，但是一般情况很能到上界，一般的情况要好的多。四个概念，一个定理：网络流的概念和定理非常多，为了不混淆，先只讲两个会用到的定理。1.反向边：这个概

import numpy as npfrom scipy import integrateimport sympy#sol:通解 ics:初始条件 x:自变量 known_params:已知参数def apply_ics(sol, ics, x, known_params): #得到未知参数 free_params = sol.free_symbols - set(known_params) #将初始条件代入零阶偏导和一阶偏导 eqs = [(sol.lhs.d..

拉格朗日插值法：可以用于少样本点的情况，一定可以返回较好的插值函数。当样本点较多时，函数复杂度急剧上升。样条插值法：低阶多项式进行拟合，误差比拉格朗日稍大，但是避免了龙格现象。...

#求解一个函数fun的极值，注意是局部最优from numpy.core.fromnumeric import argsortfrom scipy.optimize import minimizeimport numpy as npdef fun(args): #返回一个args/x+x的函数 a=args v = lambda x:a/x[0]+x[0] return vif __name__ == "__main__": args = (1) ..

#求解线性规划问题代码from scipy import optimizeimport numpy as np#规划函数的系数矩阵c = np.array([2,3,-5])#不等式约束的系数矩阵，规定是Ax<=C的形式，C是常数A = np.array([[-2,5,-1],[1,3,1]])#不等式约束的右边部分B = np.array([-10,12])#等式约束的系数矩阵Aeq = np.array([[1,1,1]])#等式约束的右边部分Beq = np.arra.