理论-数值计算方法
数值计算方法
豆沙睡不醒
永远不要停下脚步
展开
-
数值计算方法-曲线拟合问题-最小二乘法
#include<bits/stdc++.h>using namespace std;const int maxn=1100;double x[maxn],y[maxn];double s[maxn],t[maxn];double a[maxn][maxn],b[maxn];double l[maxn][maxn],u[maxn][maxn];void solve(int n){ //Doolitle分解法解线性方程组 for(int i=1;i<=n;i++) l[i]原创 2021-08-30 17:30:50 · 184 阅读 · 0 评论 -
一些实用的工具类网站
保存一波防止换个电脑又忘了1.1.1. https://paste.ubuntu.com/ 点击跳转贴代码的网站,pasteme挂了就用这个了2.2.2. http://oeis.org/ 点击跳转oeis,找规律的网站3.3.3. https://csacademy.com/app/graph_editor/ 点击跳转输入边可以生成图,支持带边权4.4.4. https://www.desmos.com/calculator?lang=zh-CN 点击跳转几何画图网站,有网格线和各种图形原创 2021-08-06 19:50:31 · 413 阅读 · 3 评论 -
转载-高精度模板(综合篇)
原文网址:https://www.cnblogs.com/liuyongliu/p/10315760.html作者:liuyongliuACM-高精度模板(综合篇)在这里,我们约定,能用int表示的数据视为单精度,否则为高精度。所有函数的设计均采用带返回值的形式。本文包含1.高精度加法2.高精度减法3.高精度乘法1)高精度乘高精度的朴素算法2)高精度乘高精度FFT优化算法3)高...转载 2021-08-04 19:39:43 · 215 阅读 · 0 评论 -
2021 年百度之星·程序设计大赛 - 初赛一、二
感觉难度比去年大,也可能是我变菜了第一场1003 鸽子传送门1004 萌新void solve(){ ll a=read,b=read; if(b>a) swap(a,b); ll tmp=a-b; ll res_max=a-b,res_min=0; if(!tmp) res_max=a,res_min=2; rep(i,2,tmp/i){ if(tmp%i==0){ res_min=i;break;原创 2021-08-01 20:40:54 · 2320 阅读 · 2 评论 -
PTA团体程序设计天梯赛-练习集 L2 网红点打卡攻略(模拟)
#include<bits/stdc++.h>using namespace std;const int maxn=1e5+10;int n,m,g[210][210];struct node{ int val,id;}a[210];int vis[210],path[210];int idx=0;bool cmp(node a,node b){ if(a.val==b.val) return a.id<b.id; return a.val<b.val;}i原创 2021-04-20 08:28:50 · 360 阅读 · 3 评论 -
PTA团体程序设计天梯赛-练习集 L2完全二叉树的层序遍历(递归)
#include<bits/stdc++.h>using namespace std;const int maxn=1e5+10;int n,a[35];void dfs(int u){ if(u>n) return ; dfs(2*u); dfs(2*u+1); cin>>a[u];}int main(){ //后序:左子树 右子树 根节点 cin>>n; dfs(1); for(int i=1;i<=n;i++){ cout原创 2021-04-20 08:24:44 · 233 阅读 · 0 评论 -
PTA团体程序设计天梯赛-练习集 L3-020 至多删三个字符 (dp)
#include<bits/stdc++.h>typedef long long ll;using namespace std;const int maxn=1e6+7;ll res,dp[maxn][4];char s[maxn];int main(){ cin>>s+1; int n=strlen(s+1); dp[0][0]=1; ///dp[i][j]表示到第i位时 删除了j位 得到的字符串个数 for(int i=1;i<=n;i++){ f原创 2021-04-20 08:15:04 · 314 阅读 · 1 评论 -
2021年度训练联盟热身训练赛第三场——C,G,I
感谢队友带飞,赛后1min过G,经典操作C——Gerrymandering(枚举)原题链接题意:求每个地区A和B的无效票,以及所有地区的 efficiency gap 之和。代码:#pragma GCC optimize(2)#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef unsigned long long ull;typedef pair<ll, ll>PLL;原创 2021-03-21 19:44:59 · 385 阅读 · 1 评论 -
CF1031B Curiosity Has No Limits(结论+枚举)
原题链接a&b+a∣b=a+b#pragma GCC optimize(2)#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef unsigned long long ull;typedef pair<ll, ll>PLL;typedef pair<int, int>PII;typedef pair<double, double>PDD;#de原创 2021-03-20 22:31:36 · 186 阅读 · 3 评论 -
Codeforces1499——C. Minimum Grid Path(思维+分奇偶+贪心)
原题链接题意:从(0,0)(0,0)(0,0)走到(n,n)(n,n)(n,n),每次只能向上或向右走,并且方向一定要改变。第iii段路的花费为len∗cilen*c_{i}len∗ci,求最小的花费。拐弯次数最多为n−1n-1n−1思路:由于每走一次后方向都要改变一次,所以就想到了分奇偶讨论,然后思路就卡住了。贪心没想出来,只想出来了个O(n3)的dpO(n^{3})的dpO(n3)的dp。下面的思路是KingZhang的,wtcl。首先,ccc数组是不能排序的,因为计算花费的时候是有顺原创 2021-03-19 01:24:59 · 3451 阅读 · 8 评论 -
Codeforces 954D. Fight Against Traffic
原题链接#pragma GCC optimize(2)#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef unsigned long long ull;typedef pair<ll,ll>PLL;typedef pair<int,int>PII;typedef pair<double,double>PDD;#define I_int llinlin原创 2021-03-02 00:34:54 · 2363 阅读 · 9 评论 -
Codeforces987D(逆向思维多源BFS)
原题链接题意:有n(<=1e5)个城市,m(1e5)条边,最多有k(1<=k<=100)个不同类型的物品,你要收集s个不同类型的物品。每个城市会产生一种类型(类型为1~k)的物品,问从每个城市出发收集够s个物品的最小步数。思路:代码:#pragma GCC optimize(2)#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef unsigned long long u原创 2021-02-28 23:14:22 · 563 阅读 · 2 评论 -
LDUOJ-瓦罗兰大陆(素数筛和哥德巴赫猜想)
(原题链接)[http://47.110.142.74/contest/1681/problem/0]题意:给一个数n,序列为1~n。现在将序列分为k组,使得每组内数的和都为质数并且k最小。如果k相等的话,优先输出第一个数小的。思路:哥德巴赫猜想,偶数可以写成两个素数的和,奇数可以写成三个素数的和。这是一定的。题目限制k最小,所以应该先考虑奇数是否能写成两个素数的和。相当于对1+2+……+n进行判断,筛完素数后,枚举其中一个或两个,判断是否合理即可。注意特判几种情况,比如n为2的时候。代码:原创 2021-02-27 23:47:20 · 157 阅读 · 0 评论