- 博客(7)
- 资源 (3)
- 收藏
- 关注
RSS订阅原创 LA 3485 Bridge
自适应辛普森公式模板。 #include #include #include #include #include #include #include #include #include #define LL long long #define CLR(a, b) memset(a, b, sizeof(a)) using namespace std; double a; double F(d
2013-12-04 18:30:27
1042
原创 poj 2480 Longge's problem
欧拉函数 或 积性函数应用。显然，我们设L = n / gcd(i, n) ，t = i /gcd(i, n)。则gcd(L, t) = 1；于是，gcd(i, n) 的个数为phi(L)。于是题目就可以用欧拉函数求解了。 #include #include #include #include #include #include #include #include #define LL lon
2013-12-03 10:36:50
1220
原创 poj 2462 Period of an Infinite Binary Expansion
欧拉定理。根据分数转换成2进制的过程，分子每次都乘2。对于循环节x，当2^x = 1(mod b)时肯定是循环节。显然当分母不能整除2的时候，即分母和2互质的话，就可以利用欧拉定理，使得2^(Euler(b)) = 1(mod b)。然后对于Euler(b)，枚举其因子，找到最小循环节就可以了。
2013-12-02 18:40:06
1009
原创 hdu 2462 The Luckiest number
欧拉定理。这道数论题好凶残啊！ 首先，由题意可以得出，(10^x - 1)/ 9 * 8 = L * p(p是一个未知数，但必定是整数)。 然后对上式进行移项处理，得：(10^x - 1) = 9 * L * p / 8。 设m = 9 * L / gcd(L, 8)，则有(10^x - 1) = m * p'。p’是必然存在的
2013-12-02 17:44:05
1481
原创 hdu 3292 No more tricks, Mr Nanguo
佩尔方程+矩阵快速幂。先暴力求出佩尔方程的特解（最小解），然后利用佩尔方程的矩阵递推式，推出第k项。 #include #include #include #include #include #include #include #define LL long long #define REP(i, n) for(int i = 0; i < n; i ++) #defin
2013-12-02 15:09:55
1107
原创 poj 2891 Strange Way to Express Integers
线性同余方程模板。求x满足，x%ai = ri。 #include #include #include #include #include #include #include #define LL long long #define REP(i, n) for(int i = 0; i < n; i ++) #define CLR(a, b) memset(a, b, sizeof(a)) us
2013-12-02 12:02:32
924
原创 poj 2154 Color
Polya + Euler优化。由于n太大，所以需要用Euler优化。设L = n / gcd(n， i)，t = i / gcd(i, n)，则gcd(L, t) = 1.又有i sum(Euler(L) * pow(n, n / L - 1))。 #include #include #include #include #include #include #incl
2013-12-02 11:26:16
1255
数位计数问题解法研究
2013-10-07
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人