poj 3252 Round Numbers 数位DP or 组合数学

POJ 3252

题意 : 给你两个数  a  < b ,a到b的二级制中 0 的个数大于 1 的个数 的数有多少个


数位DP :刚开始这样搞得,dp[i][j],i代表第几个数  j代表可以弄几个0 ,被雷到了,就是搞不出,各种烦
,后来改成dp[i][j][k],i代表几位,j代表前面有j个1,k代表前面有k个0


组合数学 :首先对一个二进制有num位的数,首先处理出有num位的有多少个,然后处理出小于num位的数


组合数学code

#include <stdio.h>
#include <string.h>
#define max(a,b) (a>b?a:b)
#define LL __int64
LL c[88][88] ;
int a[66];
LL solve(LL x)
{
	if(x == 0)
		return 0;
	int cnt[2] , i, j , num = 0;
	LL ans = 0;
	while(x)
	{
		a[++num] = x%2;
		x /= 2;
	}
	cnt[0] = cnt[1] = 0;
	for(i = 1;i <= num; i++)
		cnt[a[i]]++;
	if(cnt[0]>=cnt[1])
		ans++;
	for(i = 1;i< num; i++)
	{
		if(a[i] == 0)
			continue;
		cnt[0] = cnt[1] = 0;
		for(j = num;j > i; j--)
			cnt[a[j]]++;
		cnt[0]++;
		int m = (num+1)/2-cnt[0];
		m = max(m,0);
		for(j = m;j < i;j ++)
		{
			if(c[i-1][j]==0)
				printf("%d %d\n",i,j);
			ans += c[i-1][j];
		}
	}
	for(i = num-1;i >= 1;i --)
	{
		for(j = (i+1)/2;j < i; j++)
			ans += c[i-1][j];
	}
	return ans;
}	
int main()
{
	int i,j;
	LL x,y;
	memset(c,0,sizeof(c));
	c[0][0] = 1;
	for(i = 1;i <= 66; i++)
		for(j = 0;j <= 66; j++)
		{
			c[i][j] = c[i-1][j] + c[i-1][j-1];
		}
	scanf("%I64d%I64d", &x, &y) ;
	LL ans1 = solve(x-1);
	LL ans2 = solve(y);
	printf("%I64d\n",ans2-ans1);
	return 0;
}

数位DP code

#include <stdio.h>
#include <string.h>
int dp[66][66][66],a[66] ;
void init()
{
	memset(dp,-1,sizeof(dp));
}
int dfs(int len,int num1,int num0,int doing,int first)
{
	if(len==0)
	{
		if(first)
			return 1;
		if(num0>=num1)
			return 1;
		return 0;
	}
	if(!doing && !first && dp[len][num1][num0]!=-1)
		return dp[len][num1][num0];
	int end = doing ? a[len] : 1;
	int ans = 0;
	for(int i = 0;i <= end; i++)
	{
		if(first)
		{
			if(i==0)
			{
				ans += dfs(len-1,0,0,doing && end==i,1);
			}
			else
				ans += dfs(len-1,num1+1,num0,doing && end==i,0);
		}
		else
		{
			if(i==0)
				ans += dfs(len-1,num1,num0+1,doing && i==end,0);
			else
				ans += dfs(len-1,num1+1,num0,doing && i==end,0);
		}
	}
	if(!doing && !first)
		dp[len][num1][num0] = ans;
	return ans;
}
int solve(int x)
{
	int num = 1;
	while(x)
	{
		a[num++] = x%2;
		x /= 2;
	}
	num--;
//	puts("fuck");
	int ans = dfs(num,0,0,1,1);
	return ans;
}
int main()
{
	init();
	int x,y;
	scanf("%d%d", &x, &y);
	x--;
	int ans1 = solve(x);
	int ans2 = solve(y);
//	printf("%d %d\n" ,ans1,ans2);
	printf("%d\n", ans2-ans1);
	return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 给出一个$n\times m$的矩阵,每个位置上有一个非负整数,代表这个位置的海拔高度。一开始时,有一个人站在其中一个位置上。这个人可以向上、下、左、右四个方向移动,但是只能移动到海拔高度比当前位置低或者相等的位置上。一次移动只能移动一个单位长度。定义一个位置为“山顶”,当且仅当从这个位置开始移动,可以一直走到海拔高度比它低的位置上。请问,这个矩阵中最多有多少个“山顶”? 输入格式 第一行两个整数,分别表示$n$和$m$。 接下来$n$行,每行$m$个整数,表示整个矩阵。 输出格式 输出一个整数,表示最多有多少个“山顶”。 样例输入 4 4 3 2 1 4 2 3 4 3 5 6 7 8 4 5 6 7 样例输出 5 算法1 (递归dp) $O(nm)$ 对于这道题,我们可以使用递归DP来解决,用$f(i,j)$表示以$(i,j)$为起点的路径最大长度,那么最后的答案就是所有$f(i,j)$中的最大值。 状态转移方程如下: $$ f(i,j)=\max f(x,y)+1(x,y)是(i,j)的下一个满足条件的位置 $$ 注意:这里的状态转移方程中的$x,y$是在枚举四个方向时得到的下一个位置,即: - 向上:$(i-1,j)$ - 向下:$(i+1,j)$ - 向左:$(i,j-1)$ - 向右:$(i,j+1)$ 实现过程中需要注意以下几点: - 每个点都需要搜一遍,因此需要用双重for循环来枚举每个起点; - 对于已经搜索过的点,需要用一个数组$vis$来记录,防止重复搜索; - 在进行状态转移时,需要判断移动后的点是否满足条件。 时间复杂度 状态数为$O(nm)$,每个状态转移的时间复杂度为$O(1)$,因此总时间复杂度为$O(nm)$。 参考文献 C++ 代码 算法2 (动态规划) $O(nm)$ 动态规划的思路与递归DP类似,只不过转移方程和实现方式有所不同。 状态转移方程如下: $$ f(i,j)=\max f(x,y)+1(x,y)是(i,j)的下一个满足条件的位置 $$ 注意:这里的状态转移方程中的$x,y$是在枚举四个方向时得到的下一个位置,即: - 向上:$(i-1,j)$ - 向下:$(i+1,j)$ - 向左:$(i,j-1)$ - 向右:$(i,j+1)$ 实现过程中需要注意以下几点: - 每个点都需要搜一遍,因此需要用双重for循环来枚举每个起点; - 对于已经搜索过的点,需要用一个数组$vis$来记录,防止重复搜索; - 在进行状态转移时,需要判断移动后的点是否满足条件。 时间复杂度 状态数为$O(nm)$,每个状态转移的时间复杂度为$O(1)$,因此总时间复杂度为$O(nm)$。 参考文献 C++ 代码
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值