随笔分类 -  数论

摘要：问区间内有多少素数，确定了起始位置一直往后挪再标记就行，不知道之前为什么老是TLE，纪念一下吧。 阅读全文

摘要：看了两个多小时，大概意思我觉得莫比乌斯是一个利用一个数的能被整除的一些数，通过容斥来计算出这个数的一些信息，很高级。这个题就当做是模板吧。然后把两种形式贴出来： 第一种： 第二种： 然后是容斥的函数： 然后有一些性质： 再耍两道题巩固一下。 BZOJ 2301 MDZZ......只能说有毒。。。 阅读全文

摘要：#define MAX_INPUT_SIZE 10000 #define MAX_RESULT_SIZE (MAX_INPUT_SIZE * 2 - 1) char a[MAX_INPUT_SIZE],b[MAX_INPUT_SIZE],res[MAX_RESULT_SIZE]; void mult(char result[], char a[], char b[])//res一定要清空 { ... 阅读全文

摘要：nlogn的欧拉函数 阅读全文

摘要：看网上说是套路题，所以就安利一下咯。应该是按照论文积分做的，并不会。。。 这份代码最快。 http://blog.csdn.net/xuanandting/article/details/52577406 阅读全文

摘要：计算1-(2^n-1)*...*(2^n-(k-1))/(2^n(k-1)). 先算公约数，可以看出公约数只能是2的n次幂，求每个分子和2的n次幂的最大公约数，因为gcd(a,b)=gcd(b-a,a)，所以直接求2^n和(k-1)!的公约数，然后同除就行。这里有勒让德定理：在正数n!的素因子标准分 阅读全文

摘要：第一道FFT，研究了kuangbin神的版子，发现想不太通，不过基本的意思都懂了，就是用一个数组的权和另一个数组的权求和，计算个数。题目中有不少坑点。 阅读全文

摘要：题意是，一堆操作，I是添加一个数，D是去掉一个数，Q是询问这个数如果通过乘除质数，x次能到集合中任何一个数，问最少次数。 NUM是每个数的质因子个数，通过筛法筛出来；这里规定乘一个质数称之为一次，Ci,j是指i这个数，j次能到达集合中的数的数量。ans中是一个二进制数，表示这个可以几次到（相应位置是 阅读全文

摘要：// 质因子个数int num[MAXN];inline void init(){ num[1] = 0; for (int i = 2; i < MAXN; ++i) { if (!num[i]) { num[i] = 1; for (int j = i + i; j < MAXN; j += i 阅读全文

摘要：首先是一个公式： gcd(a^m-b^m,a^n-b^n)=a^(gcd(m,n))-b^(gcd(m,n)) (a>b) 由此可得，本题要求x^gcd(a,b)-1; 再考虑每一个gcd值(d)可能对应多个数对，而只要求出对应多少个数对(sd)，就能直接由乘法得出结论。 那么怎么求N以内最大公约数 阅读全文

摘要：存模板 乘法逆元 ON 阅读全文

摘要：盗一波素数筛 const int mx = 1000000 + 1; ///在(1,mx)的范围内寻找素数 const int sqrt_mx = (int)sqrt((double)mx); bool vis[mx]; int prime[mx / 10]; ///在mx>65000时建议写成 i 阅读全文

摘要：typedef long long LL;const int MOD=1e9+7;LL quick_mod(LL a,LL b){ LL ans=1; a%=MOD; while(b) { if(b&1) { ans=ans*a%MOD; b--; } b>>=1; a=a*a%MOD; } ret 阅读全文

摘要：#include <iostream>#include <functional>#include <algorithm>#include <complex>#include <cstdlib>#include <cstring>#include <fstream>#include <iomanip> 阅读全文

摘要：给你一个地图10*10，从1到100，问掷骰子的次数的期望，中间会有传送门从a到b。 高斯消元基础题，学了一发板子，x存放多出来的，a存放系数； 要么是dp[i]=(dp[i+1]+...+dp[i+6]+6)/6; 要么是dp[i]=dp[go[i]]; 值得注意的是哪怕是没有dp[i+6]也要加 阅读全文

摘要：#include <stdio.h> #include <string.h> int main( ) { int m,n; int i,j,k,t; while(scanf("%d",&m)==1) { if(m<0) break; n=0; int val[m],num[m];//neng fan 阅读全文