【JZOJ2700】【GDKOI2012模拟02.01】数字
2017-07-15 21:37
323 查看
Description
Data Constraint
Solution
比赛时离正解就差一步……我们打个表找找规律后发现:X*D(X)=i *i+9 *i *k,i∈[1,9],那么对于一个询问,我们可以通过枚举i来确定k的最大值计算数量。对于询问[l,r]我们可以用[1,r]-[1,l-1]来解决。
问题来了:一个x*d(x)可能在不同的i中出现(如125*d(125)=1000*d(1000))如何去重?
我们考虑i∗i+9∗i∗k=j∗j+9∗j∗li∗i−j∗j=9∗(j∗l−i∗k)找找规律找一下可能出现x*d(x)相同的i,j,一定满足i+j=9.(i+j)∗(i−j)=9∗(j∗l−i∗k)i−j=j∗l−i∗k当l=k=-1时,这是一个解。那么我们现在要讨论ax+by=c的解的数量。
若满足ax−by=ca(x+bgcd(a,b))−b(y+agcd(a,b))=c所以相邻连个x的姐之间间隔是b/gcd(a,b)。讨论减一下就好。
Code
#include<iostream> #include<cmath> #include<cstring> #include<cstdio> #include<algorithm> #define ll long long using namespace std; const int maxn=1e7+5; ll test,n,m,i,t,j,k,l,x,y,z,ans; ll gcd(ll x,ll y){ ll r=x%y; while (r) x=y,y=r,r=x%y; return y; } int main(){ // freopen("data.in","r",stdin); scanf("%lld",&test); for (;test;test--){ scanf("%lld%lld",&m,&n);m--;ans=0; for (i=1;i<=9;i++){ if (i*i>n) break; t=(n-i*i)/(9*i)+1; ans+=t; if (i<=4 && (9-i)*(9-i)<=n){ k=(9-i)/gcd(i,9-i); ans-=t/k; } } for (i=1;i<=9;i++){ if (i*i>m) break; t=(m-i*i)/(9*i)+1; ans-=t; if (i<=4 && (9-i)*(9-i)<=m){ k=(9-i)/gcd(i,9-i); ans+=t/k; } } printf("%lld\n",ans); } }
相关文章推荐
- jzoj2700 【GDKOI2012模拟02.01】数字
- 【JZOJ2702】【GDKOI2012模拟02.01】探险
- 【GDKOI2012模拟02.01】数字(找规律&&数论)
- 【JZOJ2701】【GDKOI2012模拟02.01】矩阵
- jzoj2701 【GDKOI2012模拟02.01】矩阵
- 【GDKOI2012模拟02.01】矩阵
- JZOJ 3056. 【NOIP2012模拟10.27】数字
- 【GDKOI2012模拟02.01】探险(最短路&&SPFA)
- {题解}[jzoj3084]【NOIP2012模拟11.3】超级变变变
- jzoj3082. 【NOIP2012模拟11.1】骰子游戏
- JZOJ3051. 【NOIP2012模拟10.25】单元格
- {题解}[jzoj3083]【NOIP2012模拟11.1】塔
- jzoj3083. 【NOIP2012模拟11.1】塔(加强)
- [JZOJ 3052]【NOIP2012模拟10.25】剪草
- jzoj 2570. 【NOIP2011模拟9.17】数字生成游戏
- JZOJ3052. 【NOIP2012模拟10.25】剪草
- 高中OJ3056【NOIP2012模拟10.27】数字
- JZOJ 5379. 【NOIP2017提高A组模拟9.21】Victor爱数字
- 中国剩余定理模板(JZOJ 3093. 【NOIP2012模拟11.7】合唱队形)
- JZOJ5379. 【NOIP2017提高A组模拟9.21】Victor爱数字 数位DP