做网站后期要收维护费吗,北京网站建设 性价比,做策划的工资高吗,网站 模板 下载题目 思路来源
登录 - Luogu Spilopelia
题解
参考了两篇洛谷题解#xff0c;第一篇能得出这个式子#xff0c;第二篇有比较严格的复杂度分析
结合去年蓝桥杯洛谷P9238#xff0c;基本就能得出这题的正确做法 代码
#includebits/stdc.h
#includeiostream第一篇能得出这个式子第二篇有比较严格的复杂度分析
结合去年蓝桥杯洛谷P9238基本就能得出这题的正确做法 代码
#includebits/stdc.h
#includeiostream
#includecstdio
#includemap
#includeunordered_map
#includeset
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int N1e710;
bool ok[N];
int pr[N/10],mu[N],ans[N],cnt,up;
unordered_mapint,intsmu;
unordered_mapll,llsmu2;
ll n,m;
void sieve(ll n){mu[1]1;ans[1]1;for(ll i2;iN;i){if(!ok[i]){pr[cnt]i;mu[i]-1;}for(int j0;jcnt;j){ll ki*pr[j];if(kN)break;ok[k]1;if(i%pr[j]0){mu[k]0;break; }mu[k]-mu[i];}ans[i]ans[i-1](mu[i]!0);mu[i]mu[i-1];}
}
int djsmu(int n){if(nN)return mu[n];if(smu.count(n))return smu[n];int ans1;for(int l2,r;ln;lr1){rn/(n/l);ansans-(r-l1)*djsmu(n/l);}return smu[n]ans;
}
ll cal(ll n){if(nN)return ans[n];if(smu2.count(n))return smu2[n];ll res0;for(ll l1,r,v;l*ln;lr1){vn/l/l;rsqrt(n/v);resv*(djsmu(r)-djsmu(l-1));}return smu2[n]res;
}
int main(){cinnm;sieve(n);ull ans0;for(ll l1,r,x,y;lmin(n,m);lr1){xsqrt(n/l),ysqrt(m/l);rmin(n/(x*x),m/(y*y));ans1ull*(cal(r)-cal(l-1))*x*y;}coutansendl;return 0;
}
