【数论】正约数集合

前言

一开始做了50个数据，然后内存$5GB$，不得不放弃了45个测试点，改成了5个，后面两个比较大

大意

求1到$n$的所有自然数的所有正约数集合

思路

两种做法，一种$O(n\sqrt{n})$，一种$O(nlogn)$

nlogn代码

#include<cstdio>
#include<vector>
using namespace std;int n;
vector<int>factor[500001];
void write(int x){if(x>9)write(x/10);putchar(x%10+48);return;}
void writes(int x){write(x);putchar(32);return;}//输出流加速
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n/i;j++)
      factor[i*j].push_back(i);//初始化
    for(int i=1;i<=n;putchar(10),i++)
     for(int j=0;j<factor[i].size();j++)
      writes(factor[i][j]);//输出
}

n根号n代码

#include<cstdio>
#include<cmath>
using namespace std;
int n,a[500001],len;
void write(int x){if(x>9)write(x/10);putchar(x%10+48);return;}
void writes(int x){write(x);putchar(32);return;}
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        len=0;
        for (int j=1;j*j<=i;j++)
            if (i%j==0)
             a[++len]=j,writes(j);//存下左边一半的约数
        for (int j=len;j>0;j--)//再把剩下的约数利用数学方法求出来
         if(i/a[j]>a[j])
          writes(i/a[j]);//这里大于是为了避免完全平方数
        putchar(10);
    }
}