![](https://img-blog.csdnimg.cn/20201014180756757.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
努力变成大白的小白
书读得不多,想得太多
展开
-
整数规划--分支定界法的Matlab实现
Matlab实现:A = [-1 3; 7 1];b = [6; 35];c = [-7; -9];%标准格式是求min,此题为max,需要转换一下lb = [0; 0];%x值的初始范围下界ub=[inf;inf];%x值的初始范围上界optX = [0; 0];%存放最优解的x,初始迭代点(0,0)optVal = 0;%最优解[x, fit, exitF, iter] =...原创 2020-01-03 11:16:32 · 7604 阅读 · 7 评论 -
单纯形法以及对偶单纯形法的Matlab实现
最近在上《运筹与优化》这门课,讲到了单纯形法这部分,老师让我们上机用Matlab实现单纯形法,数学公式实现过程看不懂的我就暴力模拟了单纯形法的整个过程,包含无界解等情况。A=[2,-3,2,1,0; 1/3,1,5,0,1];b=[15,20]';c=[1,2,1,0,0];m=2;n=5;format rat;%结果以分数形式显示[x_opt,fx_opt,iter]=S...原创 2019-12-20 17:05:51 · 6811 阅读 · 7 评论 -
整数规划----割平面法的Matlab实现
在纯整数规划问题中,有一种求解方法叫做割平面法,它是通过不断地求解最优解,然后每次引入一个松弛变量,将其加入约束条件中,最终求得最优解的这样的一个循环过程。根据b矩阵的值选择利用单纯形法或者对偶单纯形法求解松弛问题判断此时解是否满足整数规划问题,若满足则结束循环,否则继续下一步;从解得的b矩阵中选择和对应整数差值最大的值的行位置作为新加入的割平面约束的对应行∑\sum∑选中行中非基矩阵乘...原创 2020-01-03 11:15:40 · 7875 阅读 · 13 评论 -
E - Sigma Function【数学】【规律】
题目链接:https://vjudge.net/contest/333800#problem/E题目大意:根据唯一分解定理,计算[1,n][1,n][1,n]约数是偶数的个数。思路:被数据范围演了,打表找规律即可。我们发现:一个数是平方数的话约数和就是奇数(直接开方向下取整);一个数包含2的奇数次方的约数和也是奇数(除2之后直接开方向下取整),所以用nnn减去即是答案。AC代码:#inc...原创 2019-10-15 21:09:53 · 148 阅读 · 0 评论 -
Codeforces Round #589 (Div. 2)E. Another Filling the Grid【容斥】
题目链接:E. Another Filling the Grid题目大意:在给定n∗nn*nn∗n方格中放数字,数字范围是[1,k][1,k][1,k],问在每一行每一列的最小数字是1(每一行每一列至少有一个1)的方案数有多少种?题目思路:虽然这题数据范围小,但是比赛时想到结束都没想出来要怎么解。(神仙网友的题解,赶紧学习一波)对于给定大小的有限解,考虑容斥!于是就有:可行方案数= 总方...原创 2019-09-30 21:36:00 · 324 阅读 · 0 评论 -
Codeforces Round #338 (Div. 2) D. Multipliers【数学】【欧拉降幂】
题目链接:https://codeforc.es/contest/615/problem/D题目大意:根据唯一分解定理,给定关于nnn的所有素因子,计算关于nnn的所以因数乘积%(1e9+7)\%(1e9+7)%(1e9+7)。思路:记录每个素因子出现的次数cic_ici,预处理素因子的次数+1的前缀积pre[i]pre[i]pre[i]和后缀积suf[i]suf[i]suf[i]。枚举每...原创 2019-08-28 13:51:46 · 113 阅读 · 0 评论 -
HDU 1724【自适应辛普森公式】
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1724AC代码:#include<bits/stdc++.h>using namespace std;int gcd(int x,int y){ return y==0?x:gcd(y,x%y);}double a,b;double ll,rr;double eps=...原创 2019-08-16 20:37:03 · 105 阅读 · 0 评论 -
HDU 4497 GCD and LCM
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4497题目大意:给定两个数L,GL,GL,G,问满足gcd(x,y,z)=G,lcm(x,y,z)=Lgcd(x,y,z)=G, lcm(x,y,z)=Lgcd(x,y,z)=G,lcm(x,y,z)=L的(x,y,z)(x,y,z)(x,y,z)组合有多少种?思路:首先我们可以判断:如果LLL...原创 2019-08-24 21:03:59 · 82 阅读 · 0 评论 -
Codeforces Round #315 (Div. 1) A. Primes or Palindromes?【素数+回文数】
题目链接:http://codeforces.com/contest/568/problem/A题目大意:定义π(n)π(n)π(n)为不大于nnn的素数个数,rub(n)rub(n)rub(n)为不大于nnn的回文数个数,给定两个数p,qp,qp,q,A=pqA=\frac{p}{q}A=qp,找到最大的nnn满足π(n) ≤ A⋅rub(n)π(n) ≤ A·rub(n)π(n) ≤ A⋅...原创 2019-08-19 23:50:33 · 112 阅读 · 0 评论 -
HDU 6639 Faraway【数学】
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6639题目大意:给定n,mn,mn,m以及nnn组xi,yi,ki,tix_i,y_i,k_i,t_ixi,yi,ki,ti,且条件为(∣xi−xe∣+∣yi−ye∣) %ki=ti(|x_i-x_e|+|y_i-y_e|) \ \%k_i=t_i(∣xi−xe∣+∣yi−ye...原创 2019-08-07 23:59:14 · 218 阅读 · 0 评论