#include using namespace std; typedef long long ll; int n,tot; ll p[200],f[2000],ans; int main(){ scanf("%d",&n); for (int i=2;i=p[i];--j) for(int k=p... ...
分类：
其他好文 时间：
2019-04-22 14:02:52
阅读次数：
141
Divisibility 给定n个数，问是否能从中选出恰好k个数，使得这些数两两之差可以被m整除。 【输入格式】 第一行输入三个正数n，k，m。 接下来一行n个数。 【输出格式】 若不能选出k个数，则输出"No"（不包含引号）。 若可以，第一行输出"Yes"（不包含引号），第二行输出k正整数，用空格 ...
分类：
其他好文 时间：
2019-04-18 00:46:44
阅读次数：
155
毒瘤出题人，坑了我们好久，从基本的素数筛选,到埃氏筛法，到随机数快速素数判定，到费马小定理，好好的水题做成了数论题。 结果答案是 2*n=n+3*n，特判1，2. 以下为毒瘤题目： 题目大意： 输入一个数n, 输出两个合数（即非素数）a,b 实现 n+a=b 解题思路 3n=n+2n; 特判1、2 ...
有些事情要明着写出来才会去干。 这里是一个不断更新的Todo List，大致按照重要度和列出时间排序。 主要着眼短期计划，其中的大部分事务应该在一周内解决~~，争取不做一只鸽子~~。 填好 "模板库" 的坑。 学习树状数组、数论（矩阵快速幂、博弈论）。 写链表、hash的总结。 填上[栈&队列]中单 ...
分类：
其他好文 时间：
2019-04-14 20:40:04
阅读次数：
212
要点 998244353的原根g = 3，意味着对于任意$$1 include include include include using namespace std; typedef long long ll; const int p = 998244353, g = 3; int K, b[10 ...
分类：
其他好文 时间：
2019-04-13 19:16:45
阅读次数：
150
1 /** 2 problem: http://www.fjutacm.com/Problem.jsp?pid=3700 3 按二进制将k个待查点分类分别跑dijkstra 4 **/ 5 #include 6 #include 7 #include 8 using namespace std; 9... ...
分类：
其他好文 时间：
2019-04-12 23:06:06
阅读次数：
199
数论卷积： 对于两个数论函数f(x),g(x) f(n)g(n)=∑ f(d)g(n/d) d|n 莫比乌斯函数： 设一个数n=(p1^k1)*(p2^k2)*(p3^k3)*..........*(pn^kn) / 1,n=1 则μ(n)= 0,max(k1,...kn)=1 \ -1^r ,ma ...
分类：
其他好文 时间：
2019-04-12 20:46:04
阅读次数：
239
扩展欧几里得算法 顾名思义，扩欧就是扩展欧几里得算法，那么我们先来简单地回顾一下这个经典数论算法。 对于形如$ax+by=c$的不定方程，扩展欧几里得算法可以在$O(5log_{10}\min\{a,b\})$的时间内找到该方程的一组特解，或辅助$gcd$判断该方程无解。 对于扩欧的详细讲解，可见 ...
分类：
其他好文 时间：
2019-04-12 16:16:00
阅读次数：
163
题解 一道数论欧拉函数和欧拉定理的入门好题。 虽然我提交的时候POJ炸掉了，但是在hdu里面A掉了，应该是一样的吧。 首先我们需要求的这个数一定可以表示成$\frac{(10^x 1)}{9}\times 8$。 那么可以列出一个下面的方程 $$\frac{(10^x 1)}{9}\times 8= ...
分类：
其他好文 时间：
2019-04-12 15:04:09
阅读次数：
126
同余 同余是数论中一个重要的概念，若整数$a$与整数$b$除以正整数$m$的余数相等，则称$a$，$b$再模$m$意义下同余，记为$a\equiv b(mod\ m)$或$m|(a b)$。 同余基础性质 $1.$$a≡a (mod\ m)$，自反性 $2.$若$a≡b (mod\ m)$，则$b≡ ...
分类：
其他好文 时间：
2019-04-11 22:12:58
阅读次数：
195
