数论
文章平均质量分 65
HDACMer
I am from Soochow University
展开
-
反素数
反素数问题描述:对于任何正整数x,起约数的个数记做g(x).例如g(1)=1,g(6)=4.定义:如果某个正整数x满足:对于任意i(0现在给一个N,求出不超过N的最大的反素数.比如:输入1000 输出 840思维过程:求[1..N]中最大的反素数–>求约数最多的数如果求约数的个数 756=2^2*3^3*7^1(2+1)*(3+1)*(1+1)=24基于上述原创 2013-07-07 20:52:57 · 718 阅读 · 0 评论 -
POJ 1061 扩展欧几里得
#include#includetypedef long long ll;void gcd(ll a,ll b,ll& d,ll& x,ll& y){ if(b==0){ d=a; x=1; y=0; return ; } gcd(b,a%b,d,y,x); y-=(a/b)*x;}int main(){原创 2013-09-03 09:42:53 · 584 阅读 · 0 评论 -
HDU 1724 自适应辛普森法
//很裸的积分题,直接上模板#include#includeint aa, bb;//函数double F(double x){ return sqrt((1-x*x/aa/aa)*bb*bb);}//三点simpson法double simpson(double a,double b){ double c=a+(原创 2013-09-04 23:00:44 · 924 阅读 · 0 评论 -
SGU 107
#include#include//#define suf 987654321int main(){ //测试使用,暴力解出n=9时的8个解,平方后缀是否为suf只与其后9位有关 /*for(long long i=sqrt(suf*1.0);i if((i*i)%(1000000000)==suf){ pri原创 2013-09-04 23:01:29 · 696 阅读 · 0 评论 -
Codeforces 300C
题目链接: http://codeforces.com/contest/300/problem/C本来是道不难的题目,还是自己的数学功底不扎实。从该题又一次巩固了关于乘法逆的概念,在剩余系中,如果要做除法,如 a / b%n , 此时应该计算 a * (b在n下的逆), 而不是直接计算 a / b另外一点值得注意的是,如果 n 为一素数,那么 b 的逆就是 pow_mod(a,原创 2014-01-07 19:24:42 · 1246 阅读 · 0 评论