2017多校-8

本文解析了HybridCrystals和KillerNames两道算法题目的解决方案。HybridCrystals通过分析给定数列的特点,利用数学方法确定特定数值是否可由数列中的元素组合而成。KillerNames则探讨了N个球放入M个盒子的问题,并提供了高效的求解方法。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Hybrid Crystals
http://acm.hdu.edu.cn/showproblem.php?pid=6140

贴一个别人的题解
我们考虑到,n的范围达到了1e3,如果要从这n个数中选择一些数直接凑的话,近乎不可能完成。
题目中又写到,为降低题目难度,增加了两个约束条件:
①.a_1一定是1,它的下标一定是N。
②.【重点】给出后n-1个数的范围:a_i ≤ (j从1——i-1的)(所有下标为N的a_j之和)+ (a_i下标为L ? a_j下标为L的数值)+ (a_i下标为D ? a_j下标为D的数值)。然而这个条件有什么用呢?
它间接表明第一个下标为L的数,只能取0或1;
(在不出现N时)第二个下标为L的数只能是0或1或2;
(在不出现N时)第三个下标为L的数取值范围为0-4;
(在不出现N时)第四个下标为L的数取值范围为0-8;
……
而这范围有什么用呢?
它保证了新增加下标为L的数时,整个序列可表示的数的范围的上界随之扩大,即可以用所给的数表示0-上界中的任意一个数。
(如在不出现N时前四个下标为L的a_i分别是1,2,4,8,则我可以从这些数中选出0-4个数来表示0-15中的任意一个数)
同理,把上面的L换为D,可得到整个序列能表示的数的范围的下界。
因此,在碰到下标为L的数就更新序列所能表示的数的范围的上界,碰到下标为D时就更新下界。因为下标为N的数可正可负,所以当出现N时,更新L边界的过程中就把N看作L,更新D边界过程中就把N看作D(即同时更新上下界)。
最后看k是否在能表示的范围内就可以了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>

using namespace std;

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n,k;
        scanf("%d%d",&n,&k);
        int a[1010];
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
        }
        getchar();
        char ch;
        int sum1=0,sum2=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%c",&ch);
            if(ch=='N')
            {
                sum1+=a[i];
                sum2-=a[i];
            }
            else if(ch=='L')
            {
                sum1+=a[i];
            }
            else if(ch=='D')
            {
                sum2-=a[i];
            }
            getchar();
        }
        if(sum1>=k&&sum2<=k)
        {
            printf("yes\n");
        }
        else
        {
            printf("no\n");
        }
    }
    return 0;
}

Killer Names
http://acm.hdu.edu.cn/showproblem.php?pid=6143

N个球放M个盒子问题:https://bbs.qzzn.com/thread-14856448-1-1.html

标程

#include <assert.h>
#include <string.h>
#include <stdlib.h>
#include <stdio.h>
#include <algorithm>
#include <iostream>
#include <map>
#include <queue>
#include <string>
#include <vector>
using namespace std;
#define debug printf("%s %d\n", __FUNCTION__, __LINE__)

const int mod = 1e9 + 7;
const int maxn = 2000 + 10;

int power(int x, int times) {
    int rt = 1;
    int base = x;
    while (times) {
        if (times & 1) rt = 1LL * rt * base % mod;
        base = 1LL * base * base % mod;
        times >>= 1;
    }
    return rt;
}

int dp[maxn][maxn];
void calc_dp(int n, int m) {
    dp[1][1] = m;
    for (int i = 2; i <= n; i++) {
        for (int j = 1; j <= min(m, i); j++) {
            dp[i][j] = (1LL * dp[i - 1][j] * j + 1LL * dp[i - 1][j - 1] * (m - j + 1)) % mod;
        }
    }
}

int c[maxn][maxn];
int f[maxn];
int main(int argc, char **argv) {
    for (int i = 0; i < maxn; i++) {
        c[i][0] = c[i][i] = 1;
        for (int j = 1; j < i; j++) 
            c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
    }
    int cases;
    int n, m;
    scanf("%d", &cases);
    while (cases--) {
        scanf("%d%d", &n, &m);
        assert(1 <= n && n <= 2000);
        assert(1 <= n && n <= 2000);
        int ans = 0;
        for (int i = 1; i < m; i++) {
            if (i > n) break;
            f[i] = power(i, n);
            for (int j = 1; j < i; j++) {
                f[i] = (f[i] - 1LL * c[i][j] * f[j]) % mod;
            }
            ans = (ans + 1LL * f[i] * c[m][i] % mod * power(m - i, n)) % mod;
        }

        calc_dp(n, m);
        int dp_ans = 0;
        for (int i = 1; i < m; i++) {
            int c = dp[n][i];
            for (int j = 0; j < n; j++) 
                c = 1LL * c * (m - i) % mod;

            dp_ans = (dp_ans + c) % mod;
        }
        ans = (ans + mod) % mod;
        printf("%d\n", ans);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值