BestCoder Round #79 (div.2) 1001 1004 HDU 5663 Hillan and the girl

链接：戳这里

jrMz and angles

 Time Limit: 2000/1000 MS (Java/Others)   Memory Limit: 65536/65536 K (Java/Others)

问题描述
jrMz有两种角，第一种角都是正nn边形的内角，第二种角都是正mm边形的内角。jrMz想选出其中一些，某种角可以选多个或一个都不选，使得选出的所有角的度数之和恰好为360度。jrMz想知道这是否可能实现。
输入描述
有多组测试数据，第一行一个整数T(1≤T≤10)，表示测试数据的组数。
对于每组测试数据，仅一行，两个整数n,m(3≤n,m≤100)，之间有一个空格隔开。
输出描述
对于每组测试数据，仅一行，一个字符串，若可能实现则为Yes，若不可能实现则为No。
输入样例
3
4 8
3 10
5 8
输出样例
Yes
Yes
No
Hint
第一组数据中，jrMz可以选择1个第一种角和2个第二种角，因为90+135+135=360。
第二组数据中，jrMz可以选择6个第一种角，因为6*60=360。
第三组数据中，jrMz无法选出一些度数之和为360度的角。

思路：分别先求出两个正多边形的角度想x,y,  然后就是x*i+y*j==360  n*n枚举i,j就可以了,打的时候hack了几发  数据是 5 10

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
#include <ctime>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<iomanip>
#include<cmath>
#define mst(ss,b) memset((ss),(b),sizeof(ss))
#define maxn 0x3f3f3f3f
#define MAX 1000100
///#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef unsigned long long ull;
#define INF (1ll<<60)-1
using namespace std;
int T,n,m;
const double eps=1e-8;
int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        int x=(n-2)*180/n;
        int y=(m-2)*180/m;
        ///printf("%.2f %.2f\n",x,y);
        int flag=0;
        for(int i=0;i<=1000;i++){
            for(int j=0;j<=1000;j++){
                if(i*x+j*y==360){
                    flag=1;
                    break;
                }
            }
            if(flag) break;
        }
        if(flag) cout<<"Yes"<<endl;
        else cout<<"No"<<endl;
    }
    return 0;
}<strong>
</strong>

Hillan and the girl

 Time Limit: 12000/6000 MS (Java/Others)   Memory Limit: 524288/524288 K (Java/Others)

问题描述
“WTF！泥萌要么有妹子要么有男人，而我什么都没有！”
机房里，Hillan看着一个个人赢，突然爆发了，于是他下定决心表白。
他找到了他喜欢的妹子，妹子对他说：“好啊，我问你一个问题，如果你回答正确，我就和你在一起。请求出
的值，其中当gcd(i,j)是完全平方数时f(i,j)=0，否则f(i,j)=1   （gcd(x,y)表示x和y的最大公约数）。”
但是Hillan的智商很低，不可能回答正确，于是他找到了你。如果你回答错误的话，你会被万念俱灰的Hillan嘿嘿嘿。
输入描述
有多组测试数据，第一行一个整数T(1≤T≤10,000)，表示测试数据的组数。
对于每组测试数据，仅一行，两个整数n,m(1≤n,m≤10,000,000)，之间有一个空格隔开。
输出描述
对于每组测试数据，仅一行，一个整数，即问题的答案。
输入样例
2
1 2333333
10 10
输出样例
0
33
Hint
第一组数据中，f(i,j)显然始终为0，因为ii始终为1，gcd(i,j)始终为完全平方数（始终为1）。

思路：前段时间做了一下莫比乌斯反演的题,这道题和bzoj 2820很像 (它是gcd(i,j)=prime素数)  ,我的博客也写了这道题  可以去看看

这道题是要求gcd(i,j)=完全平方数

bzoj2820的具体分析：戳这里

我感觉还是在理解吧,直接把求素数倍数的前缀和改成完全平方数的倍数的前缀和就可以了

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
#include <ctime>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<iomanip>
#include<cmath>
#define mst(ss,b) memset((ss),(b),sizeof(ss))
#define MAX 10000001+5
#define maxn 1e6+5
///#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef unsigned long long ull;
#define INF (1ll<<60)-1
using namespace std;
int mu[MAX],prime[MAX/10+5],vis[MAX],d[MAX/10+5],num;
ll g[MAX];
int cnt;
void Moblus(){
    cnt=0;
    num=0;
    mst(vis,0);
    mu[1]=1;
    for(int i=2;i<=MAX;i++){
        if(!vis[i]){
            prime[++cnt]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=cnt;j++){
            if(i*prime[j]>MAX) break;
            vis[i*prime[j]]=1;
            if(i%prime[j]==0){
                mu[i*prime[j]]=0;
                break;
            } else mu[i*prime[j]]=-mu[i];
        }
    }
    for(int i=1;i<=MAX;i++){
        int t=(int)sqrt(i*1.0);
        if(t*t==i) d[++num]=i;
    }
    for(int i=1;i<=num;i++){
        int p=d[i];
        for(int j=1;j*p<=MAX;j++){
            g[j*p]+=mu[j];
        }
    }
    for(int i=1;i<=MAX;i++) g[i]+=g[i-1];
}
int main(){
    Moblus();
    int T,n,m;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        if(n>m) swap(n,m);
        ll ans=0;
        for(int i=1,last=0;i<=n;i=last+1){
            last=min((n/(n/i)),(m/(m/i)));
            ans+=(ll)(g[last]-g[i-1])*(n/i)*(m/i);
        }
        printf("%I64d\n",(ll)n*m-ans);
    }
    return 0;
}