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

DP Good Sequences

时间:2014-08-17 18:38:52      阅读:197      评论:0      收藏:0      [点我收藏+]

标签:blog   io   for   ar   2014   amp   log   ef   

给N个升序的数字,要求找出一个子串,每相邻两个数字不互质,求最长串的长度

提示

1)dp[i]表示到第i个数字的最长串

2)dp[i]用前i-1项中与第i项不互质的最大项更新

3)寻找与第i项不互质,即找与第i项有公公因数,所以建立数组dig[i]表示该因数出现的次数

#include <stdio.h>
#include <string.h>
#define maxn 100005
int dp[maxn];
int dig[maxn];
int mark[maxn];
int gcd(int a,int b)
{
	int k;
	while(b!=0)
	{
		k=a%b;
		a=b;
		b=k;
	}
	return a;
}
int max(int a,int b)
{
	return a<b?b:a;
}
int find(int a)
{
	int maxx=0;
	if(a==1||a==2||a==3)
	{
		dig[a]++;
		return 1;
	}
	int i=0,j=2;
	int flag=1;
	while(j*j<=a)
	{
		if(a%j==0)
		{
			dig[j]++;
			flag=0;
			if(j*j!=a)
				dig[a/j]++;
			maxx=max(max(maxx,dig[j]),dig[a/j]);
		}
		j++;
	}
	j=2;
	while(j*j<=a)             //将所有第i项的因数出现次数更新为最大值!!!
	{
		if(a%j==0)
		{
			dig[j]=maxx;
			if(j*j!=a)
				dig[a/j]=maxx;
		}
		j++;
	}
	if(flag==1)             //素数的dig[a]直接加1
	{
		dig[a]++;
		return 1;
	}
	return maxx;
}

void init()
{
	int i;
	for(i=0;i<maxn;i++)
	{
		mark[i]=0;
		dp[i]=0;
		dig[i]=0;
	}
}
int main()
{
	int n;
	int i,j,k;
	int num[maxn];
	int m;
	int tot;
	while(scanf("%d",&n)!=-1)
	{
		init();
		m=0;
		for(i=1;i<=n;i++)
		{
			scanf("%d",&num[i]);
			mark[num[i]]=i;
		}
		for(i=1;i<=n;i++)
		{
			tot=find(num[i]);
			dp[i]=tot;
			m=max(m,dp[i]);
		}
		printf("%d\n",m);
	}
	return 0;
}


DP Good Sequences,布布扣,bubuko.com

DP Good Sequences

标签:blog   io   for   ar   2014   amp   log   ef   

原文地址:http://blog.csdn.net/trq1995/article/details/38640253

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