[反演] LOJ #509. 「LibreOJ NOI Round #1」动态几何问题
2017-07-23 13:22
751 查看
μ2(n)=∑d2|nμ(d)
然后就是xjb推
反正退役了 我也就弃坑了 95分代码
复杂度分析及优化详见官方题解
然后就是xjb推
反正退役了 我也就弃坑了 95分代码
复杂度分析及优化详见官方题解
#include<cstdio> #include<cstdlib> #include<cmath> #include<algorithm> using namespace std; typedef long long ll; ll n,m; const int maxn=122474490; const int N=maxn+5; int prime[7000000],num; int mu ,mu2 ; const int P=10000007; inline void Pre(int n){ mu[1]=1; int *vst=mu2; for (int i=2;i<=n;i++){ if (!vst[i]) prime[++num]=i,mu[i]=-1; ll t; for (int j=1;j<=num && (t=(ll)i*prime[j])<=n;j++){ vst[t]=1; if (i%prime[j]==0){ mu[t]=0; break; }else mu[t]=-mu[i]; } } for (int i=1;i<=n;i++) mu2[i]=((bool)mu[i])+mu2[i-1],mu[i]=mu[i-1]+mu[i]; } inline ll S(ll n){ if (n<=maxn) return mu2 ; int x=sqrt(n); ll ret=0; for (int i=1,j;i<=x;i=j+1){ ll t=n/i/i; j=sqrt(n/t); ret+=(mu[j]-mu[i-1])*t; } return ret; } int main(){ Pre(maxn); scanf("%lld%lld",&n,&m); if (n>m) swap(n,m); ll ans=0; ll last=0,cur; for (ll i=1,j;i<=n;i=j+1){ ll t1=sqrt(n/i),t2=sqrt(m/i); j=min(n/(t1*t1),m/(t2*t2)); cur=S(j); ans+=(cur-last)*t1*t2; last=cur; } printf("%lld\n",ans); return 0; }
相关文章推荐
- 动态处理几何问题一则
- #510. 「LibreOJ NOI Round #1」动态几何问题
- DP动态规划问题 1300
- 动态规划之 0/1背包 问题
- 通过动态代理的方式解决中文乱码问题
- 数组上计数问题的动态规划-TopCoder-SRM584-DIV2-1000pt
- 常见的动态规划问题分析与求解
- android - ListView的动态刷新问题——用notifyDataSetChanged没作用
- 简单的动态规划——装箱问题
- 有效解决vue动态绑定多个class的官方实例语法无效的问题
- hibernate动态更新问题
- LinearLayout中动态添加Fragment不能填充整个Activity的问题
- 新学的插件jquery.tablesorter很好用,说说遇到的表格动态添加的数据不能排序的问题
- 动态规划3——最大公共子序列问题
- 有关动态多维数组的一个问题
- 意外解决js动态生成的html代码中一些事件失效问题
- 动态规划之投资问题
- 查libevent不能编译出动态库问题
- 动态规划0—1背包问题
- 在论坛中出现的比较难的sql问题:36(动态行转列 解析json格式字符串)