码迷,mamicode.com
首页 > 其他好文 > 详细

HH的军训(容斥)

时间:2015-11-22 18:53:53      阅读:138      评论:0      收藏:0      [点我收藏+]

标签:

1248: HH的军训

时间限制: 1 Sec  内存限制: 128 MB 提交: 95  解决: 11 [提交][状态][讨论版]

题目描述

大学里,最难忘的事情莫过于军训了,白白的HH童鞋就被无情的晒黑了,话说军训期间,最怕的就是被教官逗了,俗话说:珍爱生命,远离教官。教官总爱让同学们站方阵,由于同学很多,难免会有前面的同学挡着后面的同学,为了远离教官,大家都希望自己站的位置教官看不到自己,现在假如同学们站在一个有坐标点的矩阵上,同学们都站在每一组整数点上,而教官站在坐标(0,0)上,请问此时有多少个位置教官可以直接看到!

输入

第一行输入一个T,代表T组测试数据。

接下来每行输入两个数字M,N,对于任意一个同学的坐标位置(x,y),都有1<=x<=M<=10000,1<=y<=N<=10000.

输出

输出有多少位置教官能够看到。

样例输入

2
1 4
2 3

样例输出

4
5
题解:wa。。。无语,错了一下午。。。还是wa,先贴着;
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;
#define mem(x,y) memset(x,y,sizeof(x))
typedef long long LL;
vector<int>p;
void getp(int x){
	for(int i=2;i*i<x;i++){
		if(x%i==0){
			p.push_back(i);
			while(x%i==0)x/=i;
		}
	}
	if(x>1)p.push_back(x);
}
int getn(int x){
	int sum=0;
	/*for(int i=0;i<p.size();i++){
		cur=p[i];
		sum+=x/cur;
		for(int j=i+1;j<p.size();j++){
			cur*=-p[j];
			sum+=x/cur;
		}
	}*///这种方法不行如果是2,3,4,则2,4选不了。。。。
	for(int i=1;i<(1<<p.size());++i){
		int num=0,mult=1;
		for(int k=0;k<(int)p.size();k++){
			if(i&(1<<k)){
				mult*=p[k];
				num++;
			}
		}
		if(num&1)sum+=x/mult;//少了等号。。 
		else sum-=x/mult;
	} 
	//printf("%d\n",sum);
	return sum;
}
int main(){
	int T,M,N;
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&M,&N);
		int ans=M;
		for(int i=2;i<=N;i++){
			p.clear();
			getp(i);
			ans+=M-getn(M);
		}
		printf("%d\n",ans);
	}
	return 0;
}

  

HH的军训(容斥)

标签:

原文地址:http://www.cnblogs.com/handsomecui/p/4986334.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!