分割范围Codeforces Round #181 (Div. 2)

每日一贴,今天的内容关键字为分割范围

    标题链接

    A ,B虽然是水题,但要是没想清晰就易容掉进某个坑里。。

    C:预处理出有所可能的和,然后个每个每和分给a b两个数

    D:竞赛的时候标题没看懂,题意就是说在一个n*n的正方形格网中分割k次有几种分法,每次分完后剩下的四个子矩形必须为正方形。斟酌到最后可以分割log级别的次数,而且每分割一次,又现出了四个空白的子正方形,天然以可就想到dp,很明显的子结构啊。我们不妨从小到大结构,每多一层,正方形的巨细就扩展两倍,dp[i][j]示表前i层,经已画了k次的总案方数,那么转移的时候我们只需要枚举以后层的四个角落各自画了多少次,然后酿成一个两倍巨细的正放形(多画一步)。天然不能四重循环去枚举,间中要套一个的d【】组数来多进行一次小DP,算是DP套着DP吧。

    E:求满意p!%(a1!*a2!*a3!*...*ak!)=0;的最小的p。k的范围是100w,a[i]的范围是1000w.

    没什么难想的,就是考你有没有把握什么速快解分质因子的法方。。。。

    学了一种欧拉法筛,,可以在线性时间内筛选出素数,还可以失掉个每数的最小质因子lp[i]。

    http://wenku.baidu.com/view/542961fdba0d4a7302763ad5.html

    O(n)明证的关键是个每数只会被筛一次

    然后接下来的问题就是怎么利用lp[i]来解分这100w个数,其实这一步的算计量跟数的个数关无,跟数的巨细有关。

    我们可以用一个组数cnt[i]示表i这个素数现出的次数,那么接下来这段代码就是精华了。

for(int i = 0; i < n; i++) scanf("%d",&a[i]),cnt[a[i]]++,sum+=a[i];
	for(int i = M - 10; i >= 2; i--) cnt[i] += cnt[i+1];
	for(int i = M - 10; i >= 2; i--)
	{
		if(lp[i]!=i) cnt[lp[i]] += cnt[i];
		cnt[i/lp[i]] += cnt[i];
	}

    自己研讨研讨还是挺乐快的。。。。。

    

    每日一道理
漫漫人生路,谁都难免会遭遇各种失意或厄运。在凄风苦雨 惨雾愁云的考验面前,一个强者,是不会向命运低头的。风再冷,不会永远不息;雾再浓,不会经久不散。风息雾散,仍是阳光灿烂。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long lld;
const int maxn = 1000010;
const int M   =  10000010;
int a[maxn];
int prime[M];
int lp[M]; 
lld cnt[M];
int tot;
void init() // precalc in O(n) and get lp[i](the minimum prime factor of i)
{
	for(int i = 2; i <= M-10; i++)
	{
		if(!lp[i]){
			lp[i] = i;
			prime[tot++] = i;
		}
		for(int j = 0; j < tot && i*prime[j]<=10000000; j++) {
			lp[i*prime[j]] = prime[j];
			if(i%prime[j]==0) break;
		}
	}
}
bool check(lld mid)
{
	for(int i = 0;i < tot; i++)
	{
		lld tmp = mid,sum=0;
		while(tmp) {
			tmp /= prime[i];
			sum += tmp;
		}
		if(sum < cnt[prime[i]]) return false; 
	}
	return true;
}
int main()
{
	init();
	int n;
	scanf("%d",&n);
	lld sum = 0;
	for(int i = 0; i < n; i++) scanf("%d",&a[i]),cnt[a[i]]++,sum+=a[i];
	for(int i = M - 10; i >= 2; i--) cnt[i] += cnt[i+1];
	for(int i = M - 10; i >= 2; i--)
	{
		if(lp[i]!=i) cnt[lp[i]] += cnt[i];
		cnt[i/lp[i]] += cnt[i];
	}
	lld l = 1 , r = sum , best = -1;
	while(l <= r) {
		lld mid = l + r >> 1;
		if(check(mid)) {
			best = mid;
			r = mid - 1;
		} else {
			l = mid + 1;
		}
	}
	printf("%I64d\n",best);
	return 0;
}

文章结束给大家分享下程序员的一些笑话语录: 一个程序员对自己的未来很迷茫,于是去问上帝。
"万能的上帝呀,请你告诉我,我的未来会怎样?"
上帝说"我的孩子,你去问Lippman,他现在领导的程序员的队伍可能是地球上最大的"
于是他去问Lippman。
Lippman说"程序员的未来就是驾驭程序员"
这个程序员对这个未来不满意,于是他又去问上帝。
"万能的上帝呀,请你告诉我,我的未来会怎样?"
上帝说"我的孩子,你去问Gates,他现在所拥有的财产可能是地球上最多的"
于是他去问Gates。
Gates说"程序员的未来就是榨取程序员"
这个程序员对这个未来不满意,于是他又去问上帝。
"万能的上帝呀,请你告诉我,我的未来会怎样?"
上帝说"我的孩子,你去问侯捷,他写的计算机书的读者可能是地球上最多的"
于是他去问侯捷。
侯捷说"程序员的未来就是诱惑程序员"
这个程序员对这个未来不满意,于是他又去问上帝。
"万能的上帝呀,请你告诉我,我的未来会怎样?"
上帝摇摇头"唉,我的孩子,你还是别当程序员了")

posted @ 2013-04-30 17:40  xinyuyuanm  阅读(283)  评论(0编辑  收藏  举报