链接
思路:
求1到n中的因数个数和等价于求,设x为因子,就是求x在1到n里出现了几次,求1到n里是x的倍数的数有几个,即n/x。需要用整除分块,n/i的值是分块分部的,右端点是n/(n/i)。
代码:
#include<bits/stdc++.h>
using namespace std;
const int N=4e5+10;
const int inf=0x3f3f3f3f;
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
//#define int long long
//const ll P=2281701377;
const ll P=998244353;
const int mod=1e9+7;void solve(){ll n;cin>>n;ll ans=0;for(int i=1;i*i<=n;i++){ans+=(n/i-i+1)*2-1;}cout<<ans<<endl;
}
signed main(){ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);int t=1;cin>>t;while(t--){solve();}}
链接
思路:
即求,在分块的时候n/i是一个定值可以提取出来,然后剩下的就是i的和,这是一个等差数列求和,直接套公式。
代码:
#include<bits/stdc++.h>
using namespace std;
const int N=4e5+10;
const int inf=0x3f3f3f3f;
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
//#define int long long
//const ll P=2281701377;
const ll P=998244353;
const int mod=1e9+7;void solve(){ll n;cin>>n;//1--n [n/i]ll ans=0;for(ll l=1,r;l<=n;l=r+1){r=min(n,n/(n/l));ans+=(l+r)*(r-l+1)*(n/l)/2;}n=ans,ans=0;for(ll l=1,r;l<=n;l=r+1){r=min(n,n/(n/l));ans+=(l+r)*(r-l+1)*(n/l)/2;}cout<<ans<<endl;
}
signed main(){ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);int t=1;//cin>>t;while(t--){solve();}}