whzzt的博客

本文介绍了一种 $ O(n^{2/3} \log^{-1} n) $ 的积性函数求和方法。

好像很久没写博客了…随便写点东西吧。在很多数论题里都要用到整数分解，大家好像正常使用的都是 Pollard RhoPollard Rho\texttt{Pollard Rho} 算法，但是大家好像都不太会这个算法的正确写法(包括之前的我)，所以在这也算是普及一下吧。Pollard Rho Algorithm首先先考虑一下我们要做的事：我们想要找到一个 nn n 的非平凡因子 1&...

#include"bits/stdc++.h"#define last last_using namespace std;typedef long long LL;const int P=1000000007;struct node{ int n,m; inline LL hash(){return (LL)n<<30|m;} node(int _n=0,int _m=0){

简单数学题

一道小清新（？）的数论题

一道要用到许多数论知识的题

本蒟蒻的题解#include"bits/stdc++.h"using namespace std;typedef long long ll;#define regi register intconst int N=50005;int pri[N],cnt,mu[N],n;bool f[N];void Linear_Shaker(){ f[1]=true;mu[1]=1; for(

4522: [Cqoi2016]密钥破解Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 459  Solved: 223[Submit][Status][Discuss]Description 一种非对称加密算法的密钥生成过程如下：1.任选两个不同的质数p,q2.计算N=pq,r=(p−1)(q−1)3.选取小于r

虽然有点长还是建议看这篇博客，大意就是通过子集和变换构造出函数F使得可以快速计算∑[0,+∞)F（当其收敛），然后通过消无穷，最后得出式子就可以辣。题解 #include"bits/stdc++.h"using namespace std; const double eps=1e-6;const int N=1048576+5;double a[N],ans;int n,m,o[N]

我蒟蒻讲不清 = =#include "stdio.h"#include "algorithm"#include "iostream"#include "string.h"#include "stdlib.h"#include "math.h"#include "vector"#include "map"#include "set" using namespace std;

具体怎么反演的请看POPOQQQ大爷的莫比乌斯反演论文TATATAT然后注意，当用到函数F(i)的时候再将函数筛入树状数组里。。时间复杂度O(n log^2 n+q√n log n)如果不自然溢出会慢很多TAT#include "stdio.h"#include "iostream"#include "algorithm"#include "stdlib.h"#includ

2142: 礼物Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 1018  Solved: 426[Submit][Status][Discuss]Description一年一度的圣诞节快要来到了。每年的圣诞节小E都会收到许多礼物，当然他也会送出许多礼物。不同的人物在小E心目中的重要性不同，在小E心中分量越重的人，收到的礼物会

任务1直接快速幂，时间复杂度 O(q log y)        任务2用exgcd或者快速幂，看看费马小定理就清楚了，快速幂时间复杂度 O(q log P)        任务3比较麻烦。首先我们考虑设x=km+t，其中m为一个自定义的常数（这个好像叫baby steps giant step 也不知道什么鬼名字= =）        然后我们就得到了         根据费马小

显然如果长和宽互素那么就能看到。    所以只要求1~n-1的欧拉函数和，然后*2+1即可#include"iostream"#include"stdio.h"using namespace std;int n,eular[40005],prime[10005],ans,tmp;void Eular(){ eular[1]=1; int i,j; for (i=2;

先链popoqqq大爷的题解    然后看本蒟蒻口胡（雾），毕竟50题了qwq    嗯首先莫比乌斯反演在popoqqq大爷的题解里面已经很详细了。。    然后我就贴代码放筛法了……莫比乌斯反演一下，O(n)预处理，O(T √T)分块查询。。    然后就A了啊 - - （其实我还不是很懂啊qwq所以太弱也讲不清什么qwq#include "stdio.h"#include

[ 题外话 ：          = =看了一眼题目就知道是求最值      然后就不会做了╮(╯▽╰)╭所以，数学渣就去学了一发拉格朗日乘数法 - -]    那么上正文TAT 由于公式太多，我就直接截图哈qwq    附代码：#include "stdio.h"#include "iostream"#define rep(f,a,b) for(f=a

我们筛出N以内所有的素数，同时求出从1-N所有的欧拉函数值，最后枚举素数累加答案即可。#include "iostream"using namespace std;typedef long long ll;const int N=10000005;int n,tmp;int eular[N],prime[N];ll sum,esum[N];void Eular_All(){ ll

根据题目要求，我们要求出 ∑ gcd(i,N) 。初看这题，仿佛很难下手，因为N的范围是到Max_int，答案已经达到longlong型，不可能全部枚举。 我们不妨枚举N的约数，这样枚举的效率是O (log n)。然后我们要求出针对每一个约数k， gcd(i,N)=k的答案数。    那么我们通过如下过程：       ∵gcd(i,N)=k ∴ gcd(i/k,N/k)=1     

k次剩余