蒟蒻 lxw的博客

题目链接 题目大意： n种硬币，每种无限多，问在膜m的情况下能拼出的多少面额。设可以拼出的面额为x，每种硬币的面额为s1,s2,s3,s4….sn 则x=a1*s1+a2*s2+a3*s3+a4*s4······an*sn+am*m 方程当且仅当x%gcd(s1,s2,s3,s4…sn,m)==0时有解 由于am<=0,所以一定存在一组非负解#include<bi...

题目链接：https://nanti.jisuanke.com/t/38226sum可以分块求，F(T)和T*T*T放到一起计算前缀和就可以了。#include <bits/stdc++.h>#define rep(i, a, b) for(int i = (a); i <= (b); i++)#define ll long longusing namespa...

参考题解地址：https://blog.csdn.net/firstlucker/article/details/49336427上面倒数第二行写错了点东西，应该是 k*n/k = n.最后应该是因子的平方和-n*因子的个数。#include <bits/stdc++.h>#define rep(i, a, b) for(int i =(a); i <= ...

题目顺序与难度无关！一：2019南昌邀请赛网络赛：https://nanti.jisuanke.com/t/38226二：2015ACM/ICPC亚洲区长春站：https://vjudge.net/problem/HDU-5528 题解地址：https://blog.csdn.net/firstlucker/article/details/49336427...

基础教学博客：https://www.cnblogs.com/zhouzhendong/p/8627380.html可能退役前都学不完的VFK(跳骚国王)blog :http://vfleaking.blog.uoj.ac/slide/87#/N以内质数的个数：N/logNN的因子个数：sqrt(N)N的质因子个数：log10(N) （误）整除分块for(int l=1...

模板题，多项式乘法题目链接 0 1 2 a[0] - a[2] X a b c b[0] - b[2]________________ 0c 1c 2c 0b 1b 2...

题目链接：https://ac.nowcoder.com/acm/contest/549/J参考题解：https://www.cnblogs.com/henry-1202/p/10699252.html#_label7反演后迪利克雷卷积，然后线性筛卷积，再分块处理。预处理复杂度O(N) ，每次询问O(sqrt(N))。并没有使用杜教筛继续优化。#include<b...

题目链接：http://acm.sdut.edu.cn/onlinejudge2/index.php/Home/Index/problemdetail/pid/4219.html读错题目，一开始以为是任意两个元素互不相等。写了好长实际后才发现是x1​≠x2​,x2​≠x3​,x3​≠x4​,x4​≠x1。思路：用全集减去非法集合的并。非法集合有四个A：x1=x2 B:x2=...

题目链接：https://www.cometoj.com/contest/34/problem/A以后再也不能写rep(i, 1, sqrt(x)) 或者rep(i, 0, s.length())题目的核心点在于求一个数的因子个数以及因子和二者都是积性函数，所以用唯一分解定理处理每一个质因子对答案的贡献就可以了。当然也可以直接暴力求，但是一定不能在for循环里写sqr...

只能解决常系数线性递推式要求数域中每个的非0数存在乘法逆元#include&lt;bits/stdc++.h&gt;using namespace std;#define rep(i,a,n) for (int i=a;i&lt;n;i++)#define pb push_back#define mp make_pair#define SZ(x) ((int)(x).size(...

bzoj 2876骑行川藏套拉格朗日乘子法二分λ后二分速度贴一个很难看的代码#include&lt;bits/stdc++.h&gt;using namespace std;const int N = 1e4+100;double s[N],k[N],u[N],v[N],eu,t;double eps = 1e-14;int n;double pow2(double ...

https://www.cnblogs.com/flipped/p/5716603.html

自适应辛普森公式时间复杂度O(能过)double f(double x){ return x*x+sqrt(x); //积分函数}double simpson(double a,double b){ double c=(a+b)/2.0; return (f(a)+f(b)+4.0*f(c))*(b-a)/6.0;}double ars(double a,d...

线性筛int prime[N]; //素数表下标从1开始bool mark[N]; //mark[i]==1表示i不是素数int tot; void get_list(int n){ mark[1] = 1; for(int i = 2;i &amp;amp;amp;lt;= n;i ++){ if(!mark[i]) prime[++tot] = i; f...

题目链接：https://www.luogu.org/problem/P3389#include <bits/stdc++.h>#define rep(i, a, b) for(int i = (a); i <= (b); i++)#define per(i, a, b) for(int i = (a); i >= (b); i--)#define pb pus...