LTang

摘要： 解题思路:对A质因子分解A=p1^n1*p2^n2*...*pm^nm,则A^B=p1^(n1*B)*p2^(n2*B)*...*pm^(nm*B);A^B的因子之和为(1+p1^1+...+p1^n1*B)*(1+p2^1+...+p2^n2*B)*...*(1+pm^1+pm^2+...+pm^nm*B);1+pi^1+...+pi^m=(1+pi^m-1)/(pi-1),因此因子之和可以化成p/q的形式，且q,M互质时:(p/q)mod(M)=p*q^(-1)mod(M)q^(-1)为q的乘法逆元，M为质数时q^(-1)=q^(M-2)，所以(p/q)mod(M) = (p * q^(p 阅读全文

摘要： 解题思路:1 首先判断prime：H-numbers, 下标i表示4*i+1,j表示4*j+1,则(4*i+1)*(4*j+1)=4(4*i*j+i+j)+1,即下表i,j的乘积为下标4*i*j+i+j，判断prime方法类似判断普通prime方法，Line 10~122 判断Semi-prime，prime以及prime与no-prime的乘积肯定为非Semi-prime，Line 15~223 统计代码 阅读全文

摘要： 解题思路:　　枚举10^6内的所有素数　　诡异的是，采用10^10进制的时间比采用10^3进制时间大,不明真相　　程序中采用的是10^12进制代码 阅读全文

摘要： 这还是以前在spoj上A的一道题，不过直接在poj提交时Time Limit Exceeded，本打算有时间把算法重新优化一下，昨天无聊把输入改为scanf试试，没想到竟然过了，意料之外。简单的讲一下思路:我们先考虑平面求面积的情况:图1图2考虑与x轴平行的线段，沿y轴从下往上1）出现重叠就删除重叠部分2）出现相接部分就连接起来例 图1首先出现线段(0,0)-(10,0),长度l=10，往上遇到(0,3)-（3,3）、(7,3)-(10,3)，相对高度h=3，此时面积增加Δs=l*h=30;　　按照原理1)此时线段变为(3,3)-(7,3), l=4往上遇到(0,7)-（3,7）、(7,7)- 阅读全文

摘要： 解题思路:组合数学1-递推:C(m,n)=c(m-1,n)+C(m,n-1),时间复杂度O(m*n)2-组合，时间复杂度O(min(m,n))代码 阅读全文