The Preliminary Contest for ICPC Asia Shanghai 2019 BDL

传送门

B:离散化+差分

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
struct node{
	int x,id;
}a[3010];
int cmp(node a,node b)
{
	return a.id<b.id;
}
int main()
{
	int t;
	scanf("%d",&t);
	int tot=0;
	while(t--)
	{
		int n,m;
		scanf("%d%d",&n,&m);
		int t=0;
		for(int i=0;i<m;i++)
		{
			int x,y;
			scanf("%d%d",&x,&y);
			a[++t].id=x;a[t].x=1;
			a[++t].id=y+1;a[t].x=-1;
		}
		sort(a+1,a+1+t,cmp);
		int sum=0;
		int ans=0; 
		for(int i=1;i<t;i++)
		{
			sum+=a[i].x;
			if(sum&1)
			{
				ans+=a[i+1].id-a[i].id;
			}
		}
		printf("Case #%d: %d\n",++tot,ans);
	}
	return 0;
}

D:搜索+剪枝

首先,我们可以打表发现必然有一组是1....2.n n最大也只有3000 所以我们可以直接DFS

多重集合排列公式:\frac{n!}{k_1!k_2!..k_t!},其实k_i表示每个数重复的次数,n是数字的总数。

dfs枚举2~3000每个数字在整个序列中出现的次数(2*3000中最多有13个2,所以每个数字出现的次数很小),总出现次数记为num,乘积为mul,和为sum,那么满足题目要求的序列中1出现的次数num1=mul-sum,且只有当num1≥0&&num1+num≤3000时才能确定一个满足条件的序列, 该序列的n=num1+num,再利用多重集合排列公式求的排列数。 来自 博主

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 3005;
const ll MOD = 1e9 + 7;
const ll mod = 1e9 + 7;
ll F[N], Finv[N], inv[N];//F是阶乘,Finv是逆元的阶乘 
ll ans[N];
void init(){
    inv[1] = 1;
    for(int i = 2; i < N; i ++){
        inv[i] = (MOD - MOD / i) * 1ll * inv[MOD % i] % MOD;
    }
    F[0] = Finv[0] = 1;
    for(int i = 1; i < N; i ++){
        F[i] = F[i-1] * 1ll * i % MOD;
        Finv[i] = Finv[i-1] * 1ll * inv[i] % MOD;
    }
}
ll power(ll a,ll b)
{
	ll ans=1;
	a=a%mod;
	while(b)
	{
		if(b&1)
		ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
} 
bool dfs(ll x,ll mul,ll sum,ll num,ll INV)// num 非1的位数 
{
	if(x==1)
	{
		ll num1=mul-sum;
		if(num1>=0&&num1+num<=3000)
		{
			ans[num+num1]=(ans[num+num1]+((F[num+num1]*INV)%mod*Finv[num1])%mod)%mod;
			return 1;
		}
		else
		return 0;
	}
	dfs(x-1,mul,sum,num,INV);
	for(int i=1;mul*power(x,i)<=2*3000;i++)
	{
		int flag=dfs(x-1,mul*power(x,i),sum+x*i,num+i,INV*Finv[i]%mod);
		if(!flag)
		break;
	}
	return 1;
}
int main()
{
	init();
	int t;
	cin>>t;
	dfs(3000,1,0,0,1);
	while(t--)
	{
		int n;
		cin>>n;
		cout<<ans[n]<<endl;
	}
	return 0;
}

L:先把他们预处理出来即可

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int dp[1000020][15];
int get(int a,int b)
{
	int ans=0;
	while(a)
	{
		ans+=a%b;
		a/=b;
	}
	return ans;
}
int main()
{
	
	for(int i=2;i<=10;i++)
	{
		for(int j=1;j<=1000010;j++)
		{
			int temp=get(j,i);
			dp[j][i]=dp[j-1][i]+temp;
		}
	}
	int t;
	cin>>t;
	int tot=0;
	while(t--)
	{
		int a,b;
		scanf("%d%d",&a,&b);
		printf("Case #%d: %d\n",++tot,dp[a][b]);
	}
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值