solotzg

表示第一次接触该定理

polya计数，费马小定理

hdu 5072 Coprime 第一次

hdu 4910 Problem about GCD

hdu 4779 Tower Defense

hdu 2481 Toy

数论

反正是

想用c++来写大数，看来还是不太现实，借鉴了java版的 import java.io.*; import java.math.*; import java.util.*; /* * (2x-1)*2^(ky-1) <= n * (2x-1)*2^(ky) <= n*2 */ public class Main { public static void main(String a

hdu 4651 Partition

矩阵乘法的快速幂和乘法的快速幂 题目连接。。。 其中做了点优化，用该矩阵可由 f0,f1 --> f2,f3 减少计算次数，如果数据量再大点应该看得出来效果 最关键的是在矩阵乘法中的取模运算要选 MOD-1，在这被坑的很惨 #include #include typedef __int64 ll; const ll MOD = 1000000007; struct matrix

由题意可以得出将U和I都转化为I的方法

。。。。

找规律，模拟，计算压缩; 最大公因数，最小公倍数 #include #include #include #include #include using namespace std; typedef long long ll; /** ∑abs(i%A - i%B) 0<=i<N */ ll getGCD ( ll a, ll b ) { for ( int t; b;

数学推导，快速求素数，模拟 #include #include #include #include #include using namespace std; typedef long long ll; const int MAXN = 4000010; /** n = a0^b0 * a1*b1 * a2^b2 * a3^b3 ...... a均为素数 因子数

找规律，分析的水题 #include typedef long long ll; const ll MOD = 1000000007; /* 1 2 5 12 28 64 n = n-k f(n) = 2*f(n-1) + 2^(n-2) n>=2 f(n) = 1 n == 0 f(n) = 2 n == 1 f(n) = 2^2*f(

矩阵快速幂、幂相加是重点 还设计到一些ac自动机的dp #include #include #include using namespace std; typedef unsigned long long ULL; const int MAXN = 60; const int szGP = MAXN*MAXN*(sizeof (ULL)); struct tire { int

数学推导，矩阵快速幂 题目连接。。。 题目涉及n次方，自然想到用矩阵的快速幂来简化操作，把根号的部分和整数部分分开考虑，不难得出矩阵 [a,b; 1,a]，初始的数为[a,1] 但是这样一来最后的结果中就要涉及开根号，由于在矩阵幂运算中要mod m ，所以开根号后肯定会出错。 幸好网上有大神提供了一种猜想，即整数部分和带根号的部分向上取整的值相同，证明了下，还真是这样。。。 [(a+

看了网上各大神的解题，才发现不进行数学推导怎么都会超时。。。 选了个最简单粗暴的算法，对每个四边形的四个点(x,y)，若max(x, y) #include #include using namespace std; typedef __int64 ll; const int MAX = 50005; struct node { ll x, y; int maxnum;