【题目记录】——第 46 届 ICPC 国际大学生程序设计竞赛亚洲区域赛(上海)

文章目录


题目集地址 第 46 届 ICPC 国际大学生程序设计竞赛亚洲区域赛(上海)
放假第一次训练,效果还不错,做了DE两个签到题
本次补题:
DE(签到)
G(树上DP)
I(背包DP)
J(位运算)
K(构造)

D Strange_Fractions 思维+数学

题目地址D Strange_Fractions
题目大意:给出一个分数 p q \frac{p}{q} qp​ ,现在找到两个正整数a , b a,ba,b使得等式 p q = a b + b a \frac{p}{q}=\frac{a}{b}+\frac{b}{a} qp​=ba​+ab​ ,如果无解就输出两个0
思路:将原问题转换,设 x = a b , t = p q x=\frac{a}{b},t=\frac{p}{q} x=ba​,t=qp​ ,那么原等式就可以变成 x + 1 x = t x+\frac{1}{x}=t x+x1​=t,等价于求解方程 x 2 − t x + 1 = 0 x^2-tx+1=0 x2−tx+1=0,可得解为 x = t ± ( t 2 − 4 ) 2 x=\frac{t±\sqrt{(t^2-4)}}{2} x=2t±(t2−4) ​​,可知题目条件是有理数解,那么 ( t 2 − 4 ) \sqrt{(t^2-4)} (t2−4) ​就必须是有理数,得出来的解分别对应 a b , b a \frac{a}{b},\frac{b}{a} ba​,ab​,因此直接对给定的q,p带入方程求出解,判断是否是有理数解即可,如果是,则进一步将求得的解化为分数形式获得a,b。

AC代码:

#include <bits/stdc++.h>
#include <iostream>
#include <cstring>
#include <stdlib.h>
#include <stdio.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;
const int maxn = 1e5+10;

const int mod = 1e9+7;

ll gcd(ll a,ll b)
{
    if(a<b)
        swap(a,b);
    ll yu;
    yu=a%b;
    while(true)
    {
        a=b;
        b=yu;
        if(yu==0)
            break;
        yu=a%b;
    }
    return a;
}

void solve()
{
    bool bol=true;
    ll p,q;
	scanf("%lld%lld",&p,&q);
	ll k=p*p-4*q*q;
	ll kk=sqrt(k);
	if(kk*kk==k)
        bol=true;
    else
        bol=false;
    ll gcdd=gcd(p-kk,q*2);
    if(!bol)
        printf("0 0\n");
    else
    {
        printf("%lld %lld\n",(p-kk)/gcdd,q*2/gcdd);
    }
}

int main()
{
    int t = 1;
	scanf("%d",&t);
	while(t--)
	{
		solve();
	}
    return 0;
}

E Strange_Integers 思维

题目地址E Strange_Integers
题目大意:给定一个长为n的序列,从序列中选择m个数字,满足这m个数字中任意两个数的差的绝对值大于等于k值,问m的最大值
思路:将原序列排序,先选出一个最小值t,再从序列中找大于等于t+k的最小值,直到找遍整个序列,找到的数字的个数就是最大的m值
AC代码:

#include <bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;
const int N = 1e5+5;
int a[N];
void solve()
{
    int n,k;
    scanf("%d%d",&n,&k);
    for(int i = 1;i <= n;i++)
    {
        scanf("%d",a+i);
    }
    sort(a+1,a+n+1);
    int ans = 0;
    int t = a[1];
    int pos = 1;
    while(pos!=n+1)
    {
        t=a[pos];
        t+=k;
        ans++;
        pos = lower_bound(a+1,a+1+n,t)-a;
    }
    printf("%d\n",ans);
}

int main()
{
	int t = 1;
	while(t--)
	{
		solve();
	}
    return 0;
}

G Edge Groups 思维+数学

题目地址G Edge Groups
题目大意:给出一个n点n-1边的无向连通图,n必为奇数,现在要将n-1条边分成(n-1)/2组,每组满足的条件如下:一组只有两条边, 两条边有一个共同点,输出所有满足条件的分组方案的方案数对998244353取模后结果

思路:n (奇数)个点的树,n−1 条边分成 n − 1 2 \frac{n−1}{2} 2n−1​ 组,每组两条边并且这两条边要有一个公共点,询问分组的方案数。子树 x 中若有偶数个点,那么奇数条边必然无法分组,需要 x 连向其父亲的边。如果有奇数个点,那么偶数条边可以分组。

设 d[x] 为分组 x 子树中的边的方案数,设 x 的孩子 y 中,有 a 个 y 需要 (x,y) 这条边与 y 子树中的边配对,有 b 个 y 是 (x,y) 在 x 这里配对的。那么当 b 是奇数时还需要 x 与其父亲连接的边。也就是说,x 有 a 个子树点数是偶数,b 个子树点数是奇数。现在只需要考虑把 b 个子树两两分组即可,如果 b 是奇数,那么还需要 (x,fa) 加入其中。
d [ x ] = { f [ b ] ∗ ∏ d [ y ] f [ b + 1 ] ∗ ∏ d [ y ] d[x]=\left\{ \begin{matrix} f[b]*\prod{d[y]} \\ f[b+1]*\prod{d[y]} \end{matrix} \right. d[x]={f[b]∗∏d[y]f[b+1]∗∏d[y]​
n 个元素,每组两个分成 n 2 \frac{n}{2} 2n​ 组的方案数是 f[n],有递推式 f[n]=f[n−2]∗(n−1) 。
AC代码:

#include <bits/stdc++.h>

using namespace std;
#define rep(i,j,k) for(int i=int(j);i<=int(k);i++)
#define per(i,j,k) for(int i=int(j);i>=int(k);i--)
typedef long long ll;
const int N = 100010, mod = 998244353;
int n, sz[N];
vector<int> g[N];
ll d[N], f[N];
void dfs(int x, int fa)
{
    sz[x] = 1;
    d[x] = 1;
    int cnt = 0;
    for(auto &y : g[x])
    {
        if(y == fa) continue;
        dfs(y, x);
        sz[x] += sz[y];
        d[x] = d[x] * d[y] % mod;
        if(sz[y] & 1) cnt ++;
    }
    if(cnt & 1) cnt ++;
    d[x] = d[x] * f[cnt] % mod;
}
int main()
{
    scanf("%d", &n);
    for(int i=1; i<n; i++)
    {
        int x, y;
        scanf("%d%d", &x, &y);
        g[x].push_back(y);
        g[y].push_back(x);
    }
    f[0] = 1;
    for(int i=2; i<=n; i+=2)
    {
        f[i] = f[i-2] * (i-1) % mod;
    }
    dfs(3, 0);//这里其实可以以任意一个结点为根结点,其实以其他节点为根节点来遍历也是可以的
    printf("%lld\n", d[3]);
    return 0;
}

I Steadily Growing Steam DP

题目地址I Steadily Growing Steam
题目大意: 给出n张牌,每张牌有一个数字标签和一个价值,可以至多选择s张牌翻倍数字标签,现在从加倍之后的牌中选择两组牌,使得牌的数字标签之和相等且价值和最大
思路:可看作01背包问题,使用动态规划。
dp[i][j][k],从前i张牌中,至多使用j次技能,两组牌的ti之和的差+2600为k,k之所以定义为差+2600,是因为差最小为-2600,防止数组下标为负数,加一个偏移量,根据不同的情况:放入S不加倍,放入S加倍,放入T不加倍,放入T加倍直接进行状态转移即可

AC代码:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n,s,dp[101][101][5201],v[101],t[101];
const int inf=0x3f3f3f3f;
int main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin >>n>>s;
    for(int i=1; i<=n; i++)
        cin >>v[i]>>t[i];
    for(int i=0; i<=s; i++)//初始化
        for(int j=0; j<=5200; j++)
            dp[0][i][j]=-inf*(j!=2600);
    //这个地方不能动2600,因为2600为解,一开始解为0
    for(int i=1; i<=n; i++)
        for(int j=0; j<=s; j++)
            for(int k=0; k<=5200; k++) {
                dp[i][j][k]=dp[i-1][j][k];
                if(k>=t[i])dp[i][j][k]=max(dp[i][j][k],dp[i-1][j][k-t[i]]+v[i]);
                //装入T
                if(k+t[i]<=5200)dp[i][j][k]=max(dp[i][j][k],dp[i-1][j][k+t[i]]+v[i]);
                //装入S
                if(j&&k>=2*t[i])dp[i][j][k]=max(dp[i][j][k],dp[i-1][j-1][k-2*t[i]]+v[i]);
                //加倍后装入S
                if(j&&k+2*t[i]<=5200)dp[i][j][k]=max(dp[i][j][k],dp[i-1][j-1][k+2*t[i]]+v[i]);
                //加倍后装入T
            }
    cout <<dp[n][s][2600]<<endl;
    return 0;
}
上一篇:日常记录(46)寄存器模型、工厂、回调


下一篇:algs4 1.3.46栈可生成性问题中禁止出现的排列