首页 > 代码库 > BZOJ2820: YY的GCD
BZOJ2820: YY的GCD
题解:
面对这么神的题,蒟蒻只好搬运题解。。。
跪烂iwtwiioi!!!
http://www.cnblogs.com/iwtwiioi/p/4132095.html
代码:
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #include<iostream> 7 #include<vector> 8 #include<map> 9 #include<set>10 #include<queue>11 #include<string>12 #define inf 100000000013 #define maxn 10000000+514 #define maxm 100000+515 #define eps 1e-1016 #define ll long long17 #define pa pair<int,int>18 #define for0(i,n) for(int i=0;i<=(n);i++)19 #define for1(i,n) for(int i=1;i<=(n);i++)20 #define for2(i,x,y) for(int i=(x);i<=(y);i++)21 #define for3(i,x,y) for(int i=(x);i>=(y);i--)22 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)23 #define mod 100000000724 using namespace std;25 inline int read()26 {27 int x=0,f=1;char ch=getchar();28 while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}29 while(ch>=‘0‘&&ch<=‘9‘){x=10*x+ch-‘0‘;ch=getchar();}30 return x*f;31 }32 int tot,p[maxn],mu[maxn];33 ll g[maxn];34 bool v[maxn];35 void get()36 {37 mu[1]=1;38 for2(i,2,10000000)39 {40 if(!v[i])p[++tot]=i,mu[i]=-1,g[i]=1;41 for1(j,tot)42 {43 int k=i*p[j];44 if(k>10000000)break;45 v[k]=1;46 if(i%p[j])mu[k]=-mu[i],g[k]=mu[i]-g[i];47 else {mu[k]=0;g[k]=mu[i];break;}48 }49 }50 for1(i,10000000)g[i]+=g[i-1];51 }52 int main()53 {54 freopen("input.txt","r",stdin);55 freopen("output.txt","w",stdout);56 get();57 int T=read();58 while(T--)59 {60 ll n=read(),m=read(),ans=0;61 if(n>m)swap(n,m);62 for(int i=1,pos;i<=n;i=pos+1)63 {64 pos=min(n/(n/i),m/(m/i));65 ans+=(g[pos]-g[i-1])*(n/i)*(m/i);66 }67 printf("%lld\n",ans);68 } 69 return 0;70 }
2820: YY的GCD
Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 782 Solved: 391
[Submit][Status]
Description
神犇YY虐完数论后给傻×kAc出了一题
给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对
kAc这种傻×必然不会了,于是向你来请教……
多组输入
Input
第一行一个整数T 表述数据组数
接下来T行,每行两个正整数,表示N, M
Output
T行,每行一个整数表示第i组数据的结果
Sample Input
2
10 10
100 100
10 10
100 100
Sample Output
30
2791
2791
HINT
T = 10000
N, M <= 10000000
BZOJ2820: YY的GCD
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。