#include<iostream>
#include<algorithm>
#include<string>
#include<cmath>
#include<cstdio>
using namespace std;
int n,a[10]={2,3,5,7},b[5000];
int c[20000];
void asd(int k,int s)
{
for(int i=1;i<=9;i++)
{
int q=0;
for(int j=1;j*j<=k*10+i;j++)
asd(k*10+i,s+1);
}
return;
}
int main()
{
//freopen("crypt1.in","r",stdin);
//freopen("crypt1.out","w",stdout);
cin>>n;
for(int i=0;i<4;i++)
asd(a[i],1);
return 0;
}

真好用

吼吼吼...map慢出翔..换了bool数组。。0.281sVS>34s

回复 @Mike is Fool :
这..

表，提供检查用
by Clairs

啊啊啊……
卡我最后一个点！！！！！
话说漆子超的论文中神奇方法为啥会快？

这个保证有解是肯定的，因为我找了一些非常大的数然后随机一堆数取模，所以-1骗不了分。
不是质数的情况就要分解质因数取所有质数的指数的最高项即可（因为$x$ $mod$ $a$ $= c$,$x$ $mod$ $b$ $=$ $c$,则$x$ $mod$ $lcm(a,b) = c$,$lcm$为最小公倍数）而$(2^1,2^5,...... 2^x)$最小公倍数肯定是$2^{max(x)}$，所以我们取新的$P_i$为$2^{max(x)}$即可，然后让使得取得最高项的$A_i$ $mod$ 新的$P_i$作为新的$A_i$,这样的话所有的$P_i$必定互质，构造出新的方程后就按一般互质的情况计算即可
例如下面一组数据：
10
40 39
60 19
14 1
95 39
9 7
85 59
87 55
88 63
96 31
5 4
我们进行转换后得
32 31 //2^5 from 96 31
9 7 //3^2 from 9 7
5 4 //5^1 from 5 4
7 1 //7^1 from 14 1
11 8//11^1 from 88 63(63 mod 11 =8)
17 8//17^1 from 85 59(59 mod 17=8)
19 1//19^1 from 95 39(39 mod 19 =1)
29 26//29^1 from 87 55(55 mod 29=26)

读了N遍题，一直认为是按顺序建造。。。。。怎么都算不出来样例答案

kmp大法吼！~

好不科学，树链剖分竟然比倍增慢

一水就过，测试高精度模板用

粗鄙刷榜@coolkid

我是sb，让merge给挂了，真是智障！

回复 @Janis :
有吗{{莫不是来逗我玩的。。

题解有三种做法，一种是$O(n \log^2 n+q)$,一种是$O(n +q\log^2 n)$,一种是$O(n+q)$,

打表大法好！！！

解方程
$ f(x)=0 $
则可以选一个初始值$x_{0}$不断进行迭代
$ x_{n+1}=x_{n}-\frac{f(x)}{f'(x)}\ $
牛顿迭代法，维基百科上的比较详细

150纪念