一句话代码

for(int i=1,j;i<=n;i=j+1) //求Σn/i
    j=n/(n/i),sum+=(j-i+1)*(n/i);

int inv(int a,int p){  ///求a关于p的逆元
    return a==1?1:(1LL*p*(a-inv(p%a,a))+1)/a%p;
}        

for(int i=1;a[i]--==0;a[i++]=9); //高精度a--

bitset<N>a[N]
rep(j,1,n)rep(i,1,n)
    if(a[i][j])a[i]|=a[j];//bitset求传递闭包