摘要：
#include #include #include #include #include #include using namespace std; int n,m,jdg[10000001],prime[10000001],cnt; int main() { cin>>n>>m; jdg[1]=1; for(int i=2;in) bre... 阅读全文
摘要：
#include #include #include #include #include #include using namespace std; int jdg[10000001],n,m; int main() { cin>>n>>m; jdg[1]=1; for(int i=2;i>a; if(jdg[a]==0) cout... 阅读全文
摘要：
//三分 //劣于二分 但可以解决某些二分无法解决的问题 //求函数顶点可用 #include #include #include #include #include #include #define eps 1e-7 using namespace std; int n; double l,r,a[14]; double get(double x) { double num=0; ... 阅读全文
摘要：
//卢卡斯定理 可求在模p意义下的组合数 //公式：C(x,y)=C(x/p,y/p)*C(x%p,y%p) (mod p) #include #include #include #include #include #include #define ll long long #define ull unsigned long long #define INF 0x3f3f3f3f ll T,... 阅读全文
摘要：
//KMP //求匹配子串 #include #include #include #include #include #include using namespace std; int lena,lenb,next[1000001];//next[a]存放在a长度的子字符串中前后缀相同最长的长度 char a[1000001],b[1000001]; void getnxt() { ... 阅读全文
摘要：
#include #include #include #include #include #include using namespace std; long long a,b,MOD,ans=1; int main() { scanf("%lld%lld%lld",&a,&b,&MOD); while(b) { if(b&1) a... 阅读全文
摘要：
//矩阵快速幂 //满足c[x][y]=∑(i=1,i #include #include #include #include #include #define ll long long #define MOD 1000000007 using namespace std; ll n,k; struct uio{ ll squ[101][101],n,m; }mat; uio mult... 阅读全文
摘要：
//Manacher //可求最长回文子串长度 #include #include #include #include #include #include #define ma 11000005 using namespace std; int n,ans,length[mamxr) { mxr=length[i]+i; ... 阅读全文
摘要：
//高斯消元 //可解n元一次方程 //消成左上到右下的对角线 精度相对于上三角损失更大 #include #include #include #include #include #include #define eps 1e-7 using namespace std; int n; double coe[101][101],num[101]; int main() { scan... 阅读全文
摘要：
#include #include #include #include #include #include #include #include using namespace std; const int MAX_L=2005; //最大长度，可以修改 class bign { public: int len, s[MAX_L];//数的长度，记录数组 //构造函数 big... 阅读全文