const int maxn = 1e6+50;
bool vis[maxn];
int pr[maxn];
int len = 0;
void Seve_prime()
{
memset(vis, true, sizeof(vis));
vis[0] = vis[1] = false;
for(int i = 2;i*i < maxn;i++)
{
if(vis[i])
{
for(int j = i*i;j <maxn;j+=i)vis[j]= false;
}
}
for(int i = 2;i < maxn;i++)
{
if(vis[i])pr[len++] = i;
}
}
ll decom(ll a, set<int>&st)
{
for(int i = 0;a!=1&&pr[i]<= sqrt(a);i++)
{
while(a%pr[i] == 0)
{
a/=pr[i];
st.insert(pr[i]);
}
}
if(a > 1)st.insert(a);
}
数论9——分解质因子
最新推荐文章于 2024-01-19 05:30:00 发布