zsben - 博客园

uva11754 中国剩余定理+暴力搜索

摘要：是当y的组合数较小时，暴力枚举所有组合，然后用中国剩余定理求每种组合的解，对解进行排序即可注意初始解可能是负数，所以如果凑不够S个，就对所有解加上M,2M。。。。当y的组合数较大时，选择一个k/x最小的序列，枚举N=x*t+y，外层枚举t，内层枚举y，然后验证N是否是可行解为什么要选k/x最小阅读全文

posted @ 2019-02-26 17:28 zsben 阅读(241) 评论(0) 推荐(0) 编辑

中国剩余定理

摘要： https://blog.csdn.net/niiick/article/details/80229217 拓展中国剩余定理注意拓展欧几里得只能解出a,b,x,y,gcd(a,b)的一组特解模板1 模板2：拓展定理阅读全文

posted @ 2019-02-25 23:14 zsben 阅读(170) 评论(0) 推荐(0) 编辑

摘要： /* 给定n，对于所有的对（i,j）,if[x*t]+=x*phi[t] 筛出phi表 */ #include using namespace std; #define maxn 4000005 //#define ll long long bool check[maxn+10]; int phi[maxn+10],prime[maxn+10],tot; void init(){ m... 阅读全文

posted @ 2019-02-25 13:58 zsben 阅读(149) 评论(0) 推荐(0) 编辑

lightoj1197 素数双筛，可以参考poj的那题双筛

摘要： /* 判断一个数是否是素数，只要判断这个数有没有在[2,sqrt(n)]区间的因子同样，对于大数短区间的筛选，同样可以用这种判断方式，先筛出sqrt(n)范围内的素数，然后用这些素数去筛出区间内的因子 */ #include using namespace std; #define maxn 1000000 #define ll long long int v[maxn],prime[m... 阅读全文

posted @ 2019-02-25 07:26 zsben 阅读(207) 评论(0) 推荐(0) 编辑

lightoj1214 大数取模模板

摘要： #include using namespace std; #define maxn 300 #define ll long long ll a,b; char s[400]; inline ll calc(int i){ return s[i]-'0'; } int main(){ int T; cin>>T; for(int tt=1;tt>s>>b; ... 阅读全文

posted @ 2019-02-20 19:59 zsben 阅读(168) 评论(0) 推荐(0) 编辑

lightoj 1220 唯一分解定理

摘要： #include using namespace std; #define maxn 1000005 #define ll long long int v[maxn],prime[maxn],m; void init(int n){ memset(v,0,sizeof v); memset(prime,0,sizeof prime); m=0; for(int ... 阅读全文

posted @ 2019-02-20 19:21 zsben 阅读(173) 评论(0) 推荐(0) 编辑

lightoj1234 打表技巧：分块打表

摘要： /* 打不了那么大的表，所以只记录分块的信息即可 */ #include using namespace std; #define maxn 1000005 double arr[maxn]; void init(){ arr[0]=0; double tmp=0; for(int i=1;i>T; while(T--){ cin>>n; ... 阅读全文

posted @ 2019-02-20 13:12 zsben 阅读(258) 评论(0) 推荐(0) 编辑

light1236 素数打表，质因数分解

摘要：不知道为什么会错这是网上的ac代码阅读全文

posted @ 2019-02-19 22:20 zsben 阅读(278) 评论(0) 推荐(0) 编辑

lightoj1259 线性筛的另一种写法 v变成bool标记数组

摘要：也是用线性筛，但是v用int会爆，所以这个线性筛用的是另外一种写法阅读全文

posted @ 2019-02-19 08:40 zsben 阅读(209) 评论(0) 推荐(0) 编辑

lightoj 1282 取对数的操作

摘要： /* 前三位 len=log10n^k(乘积的长度) len=klog10n n^k=x*10^(len-1) x=n^k/10^(len-1) log10x = k*log10n - (len-1) x=pow(10,k*log10n - (len-1)) 后三位快速幂解决 */ #include using namespace std; #define ll long long l... 阅读全文

posted @ 2019-02-16 11:12 zsben 阅读(181) 评论(0) 推荐(0) 编辑