- 2024-11-188
8.4点击查看代码#学号17importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivpdefsystem(t,state):x,y=statedxdt=-x-ydydt=x-yreturn[dxdt,dydt]t_span=(0,30)y0=[1,0.5]sol=so
- 2024-11-178.4 8.5
8.4点击查看代码importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivpdefsystem(t,state):x,y=statedxdt=-x*3-ydydt=x-y*3return[dxdt,dydt]t_span=(0,30)y0=[1,0.5]sol=solve_ivp(
- 2024-11-17The sol to print
Thesoltoprinthttps://oier.team/problems/93思路用两个优先队列。一个用于存储没有打印任务的打印机,一个存储有任务的打印机。如果有打印机没有打印任务直接选择里面最小的。否则,找到等待时间最小的那一个。Code#include<bits/stdc++.h>usingnamespacestd;consti
- 2024-11-17第八章
8.4importnumpyasnpfromscipy.integrateimportodeintimportmatplotlib.pyplotaspltdefsystem(z,t):x,y=zdxdt=-x3-ydydt=x-y3return[dxdt,dydt]z0=[1,0.5]t=np.linspace(0,30,1000)sol=odeint(system,z0,t)x_sol=sol[:,0]y_sol=sol[:,1]plt.sub
- 2024-11-17数学建模习题8.5
`importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivp定义微分方程模型defmodel(t,y):f,df_dm,d2f_dm2,T,dT_dm=yd3f_dm3=-3*f*d2f_dm2+2*(df_dm)**2-Td2T_dm2=-2*f*dT_dmreturn[df_dm,d2f_dm2,d3f_dm
- 2024-11-17习题8.5
1.代码实现点击查看代码importnumpyasnpfromscipy.integrateimportodeintimportmatplotlib.pyplotasplt#定义微分方程组defsystem(state,eta):f,df,d2f,T,dT=stated3f=-3*d2f+2*df**2-Td2T=-2.1*df*dTreturn[df,
- 2024-11-17习题8.4
1.代码实现点击查看代码importnumpyasnpfromscipy.integrateimportodeintimportmatplotlib.pyplotasplt#定义微分方程组defsystem(state,t):x,y=statedxdt=-x**3-ydydt=x-y**3return[dxdt,dydt]#初始条件x0=1y
- 2024-11-17数学建模习题8.4
`importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivp定义微分方程系统defsystem(t,state):x,y=statedxdt=-x-ydydt=x-y**3return[dxdt,dydt]定义时间跨度和初始条件t_span=(0,30)y0=[1,0.5]解微分方程so
- 2024-11-17The sol to Bismuth / Linear Sieve
ThesoltoBismuth/LinearSievehttps://www.luogu.com.cn/problem/P11169思路因为懒惰,所以直接转载了当时参考的这篇博客https://www.luogu.com.cn/article/ys9ualoj首先观察样例发现第一个样例一定是\(\lceil\frac{n}{2}\rceil\)。手推一下加入数字的过程:当处理到
- 2024-11-178.5
importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivpplt.rcParams['text.usetex']=Falsedefmodel(t,y):f,df_dm,d2f_dm2,T,dT_dm=yd3f_dm3=-3*f*d2f_dm2+2*(df_dm)**2-Td2T_dm2=-2.1*f*dT
- 2024-11-178.4
importnumpyasnpfromscipy.integrateimportsolve_ivpimportmatplotlib.pyplotaspltdefdifferential_equations(t,z):x,y=zdx_dt=-x**3-ydy_dt=x-y**3return[dx_dt,dy_dt]z0=[1,0.5]t_span=(0,30)sol=solve_ivp(differential_eq
- 2024-11-16第八章课后习题
习题8.4点击查看代码fromscipy.integrateimportodeintimportnumpyasnpimportmatplotlib.pyplotasplt#设置Matplotlib不使用LaTeXplt.rc('font',size=15)plt.rc('text',usetex=False)#定义微分方程系统dz=lambdaz,t:[-z[0]**3-z[1],z[0]-
- 2024-11-13第八章
8.4求微分方程组的数值解x'=-x3-y,x(0)=1,y'=x-y3,y(0)=0.5,0<=t<=30,要求画出x(t)和y(t)的解曲线图形,再相平面上画出轨线点击查看代码importmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivpdefsystem(t,state):x,y=statedxdt=-x**3-y
- 2024-11-13The sol to pairing
Thesoltopairinghttps://www.luogu.com.cn/problem/P11187思路把答案序列中相邻而相等的两个数,我们称之为“块”。那么可以发现,对于以某块为结尾的一个答案序列,其一定是由一个结尾不为该块的序列转移而来。因而,本题具有最优子结构性质,可以使用动态规划求解。\(1.\)对于
- 2024-11-13The sol to expr
Thesoltoexprhttps://www.luogu.com.cn/problem/P11186思路递归加上一点点优化……重点:根据区间范围预处理出所有情况。具体的看以下题解https://www.luogu.com.cn/article/i0oqvl99e……水完了Code#include<bits/stdc++.h>usingnamespacestd;constintN=1e5+5;
- 2024-11-128.4 求微分方程组的数值解 x'=-x^3-y,x(0)=1,y'=x-y^3,y(0)=0.5,0<=t<=30,要求画出x(t)和y(t)的解曲线图形,再相平面上
importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.integrateimportsolve_ivpdefsystem(t,state):x,y=statedxdt=-x3-ydydt=x-y3return[dxdt,dydt]t_span=(0,30)y0=[1,0.5]sol=solve_ivp(system,t_span,y0,t_eval=np.linsp
- 2024-11-10Sol - P2900 [USACO08MAR] Land Acquisition G
完整准确地理解FlushHu的题解。0x00初步分析我们发现对于矩形\(i,j\)满足\(h_i\leqh_j,w_i\leqw_j\),那么选\(j\)的时候一定可以并购\(i\),因此将\(i\)删去。将剩下的矩形按照\(h\)从大到小排序,此时\(w\)从小到大。因为如果合并的不是一段连续区间,那么中间未被
- 2024-11-08讲座の题解
讲座配套题单的题解喵每题的文字解释会逐渐补充,如果有疑问直接私信喵目录讲座配套题单的题解喵目录A-看看你会不会用电脑B-求求你不要用内置函数C-GPAD-minE-for循环大神F-居然有人说这个是线性代数G-高三同学秒了H-无穷级数I-不要用内置函数
- 2024-10-28薄壁铣削中 chatter 振动的建模及其对工件表面粗糙度的影响研究(Matlab代码实现)
- 2024-10-24The sol of coin(SA)
Thesolofcoin(SA)https://www.luogu.com.cn/problem/P3878这题是模拟退火的板子。虽然一般的模拟退火不能支持多测的题目,因为时间可能卡不好,但在没有办法的办法下(其实就是乱搞),可以自己酌情设置SA的参数和SA的次数。AC_code#include<random>#include<algorithm>#defineY
- 2024-10-24The sol of coin(搜索减脂版)
Thesolofcoin(搜索减脂版)https://www.luogu.com.cn/problem/P3878这题是模拟退火的板子,但这里先讲搜索(刚好练练搜索)搜索减脂\(1.\)按价值从大到小排序,你一不小心取的价值太大会被剪枝\(2.\)最多取n/2个金币,你取得太多是要被剪枝的codewith注解#include<bits/stdc++.
- 2024-10-18arc板刷
RT记录一下板刷的过程a.aarc058a题意买家想买一个价格为N的物品,但他又讨厌k个数字,分别为D_1,D_2,……,D_K。问他最少出多少钱,才能在保证买下这个物品的同时使自己出的钱不包括自己讨厌的数字。sol枚举即可submissionb问题陈述我们有一个行数为\(H\)列数为\(W\)的大正
- 2024-10-04At_pakencamp_2023_day1_p sol
题面给你两个序列\(A,B\),\(\forallu,v(u\not=v)\)之间边的权值为\(a_ua_v+b_ub_v\)。求最小生成树的边权和。原题目editorial朴素的想法考虑类似题目的做法,考虑每一次寻找最小的然后加入。发现这种思想和Boruvka比较相似。于是我们考虑Boruvka的方式来做。对现有的连
- 2024-10-04openzeppelin/contracts/utils/Counters.sol" not found
运行以下//SPDX-License-Identifier:MITpragmasolidity0.7.5;import{Counters}from"@openzeppelin/contracts/utils/Counters.sol";contractMyTokenisERC721,Pausable,Ownable{usingCountersforCounters.Counter;Counters.Cou
- 2024-09-24CF1270 sol
题目大意cflink给定一个长度为\(n\)的序列\(a\),保证任何时间序列\(a\)两两不同,\(i\)和\(j\)有边当且仅当\(a_i<a_j\)。询问连通块的个数,带单点修。做法step1观察性质结论:若\(i,j\)连通,则\(\forallk\in(i,j),i,j\)和\(k\)连通。$\mathcal{proof}:$分讨:\(a_i