2693: jzptab

来源：互联网发布：网络对大学生的好处编辑：程序博客网时间：2024/06/06 07:32

吐血！！！

其实理论上我本能A掉这题，但由于中途少算了个x^2，又烧了好几个小时，最后还是在连击哥的帮助下、看了答案才查出错来。

表示已经失去自己解决问题的能力了。。。

讲道理以前考虑过报数竞，但由于怕算错，果断报了信竞，现在这些数学题让我情何以堪？

不管怎样，还是感谢Xrzl的耐心教导。

恶心的代码：

#pragma GCC optimize("O3")#include<cstdio>#include<iostream>#include<algorithm>#include<cmath>using namespace std;#define rep(i,j,k) for(i=j;i<=k;++i)#define per(i,j,k) for(i=j;i>=k;--i)#define ll long long#define pli pair<ll,int>#define mkp make_pair#define X first#define Y second#define N 10000005#define MOD 100000009ll _,n,m,k,ans;inline ll add(ll a,ll b){return (a+b)%MOD;}inline ll sub(ll a,ll b){return (a-b)%MOD;}inline ll mul(ll a,ll b){return a*b%MOD;}inline ll cal(ll a){return a*(a+1)/2%MOD;}ll pri[N],cnt,mu[N],Sf[N],f[N];bool vis[N];void shai(){ll i,j,v,rg=N-5;mu[1]=f[1]=Sf[1]=1;rep(i,2,rg){if(!vis[i])mu[pri[++cnt]=i]=-1,f[i]=1-i;for(j=1;j<=cnt&&(v=pri[j]*i)<=rg;j++){vis[v]=1;if(i%pri[j]==0){f[v]=f[i];break;}mu[v]=-mu[i];f[v]=mul(f[i],1-pri[j]);}Sf[i]=add(f[i]=mul(f[i],i),Sf[i-1]);}}int main(){shai();scanf("%lld",&_);ll i,j,rg;while(_--){scanf("%lld%lld",&n,&m);rg=min(n,m);ans=0;for(i=1;i<=rg;i=j+1){j=min(n/(n/i),m/(m/i));ans=add(ans,mul(mul(cal(n/i),cal(m/i)),sub(Sf[j],Sf[i-1])));}printf("%lld\n",add(ans,MOD));}return 0;}

0 0