luogu P5289 [十二省联考2019]皮配

传送门

首先考虑一个正常的dp,设\(f_{i,j,k}\)为前\(i\)个学校,\(j\)人在\(\color{#0000FF}{蓝阵营}\),\(k\)人在\(\color{#654321}{吔}\)派系的方案,转移枚举选哪个导师就好了,注意一个城市要选同一阵营,所以可以多开一维\(0/1\)表示当前城市在哪个阵营

\(k=0\)的情况,可以发现选\(\color{#654321}{吔}\)派系的\(Yazid\)和\(Zayid\)都会增加\(\color{#654321}{吔}\)派系人数,另一个派系亦然,所以选阵营和选派系是相对独立的,所以可以dp出\(g_{i,j}\)为前\(i\)个城市\(j\)人在\(\color{#0000FF}{蓝阵营}\)的方案,以及\(h_{i,k}\)为前\(i\)个学校,\(k\)人在\(\color{#654321}{吔}\)派系的方案,然后两者之和的乘积就是答案

\(k\le 30\)的情况,一个显然的想法是把所有有限制的城市单独拿出来做\(f\)数组的dp,然后其他的做\(g\)和\(h\)的dp,然后方案拼起来.至于拼起来,就是枚举\(f\)的\(j\)和\(k\),然后其他城市的维度就会限制在一个范围内,就是区间和乘起来,可以前缀和求.

只不过可能这些城市的学校贼多,,,不过那些城市的没限制学校还是可以乱选阵营的,所以可以把选阵营的贡献算到\(h\)里.具体说就是把\(f_{i,j,k}\)改为前\(i\)个限制的城市的学校,\(j\)人在\(\color{#0000FF}{蓝阵营}\),限制学校中\(k\)人在\(\color{#654321}{吔}\)派系的方案,然后转移是碰到无限制学校就转移\(h\)

// luogu-judger-enable-o2
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<cmath>
#include<ctime>
#include<queue>
#include<map>
#include<set>
#define LL long long
#define db double

using namespace std;
const int N=2500+10,M=300+10,mod=998244353;
int rd()
{
    int x=0,w=1;char ch=0;
    while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*w;
}
void ad(int &x,int y){x+=y,x-=x>=mod?mod:0;}
int f[2][2][N][M],g[2][N],h[2][N];
int n,c,c0,c1,d0,d1,kk,m;
bool ht[N];
struct sch
{
    int b,s,m;
    bool operator < (const sch &bb) const {return ht[b]!=ht[bb.b]?ht[b]<ht[bb.b]:b<bb.b;}
}a[N];

int main()
{
    int T=rd();
    while(T--)
    {
        memset(f,0,sizeof(f));
        memset(g,0,sizeof(g));
        memset(h,0,sizeof(h));
        n=rd(),c=rd();
        for(int i=1;i<=c;++i) ht[i]=0;
        c0=rd(),c1=rd(),d0=rd(),d1=rd();
        for(int i=1;i<=n;++i) a[i].b=rd(),a[i].s=rd(),a[i].m=-1;
        kk=rd();
        for(int i=1;i<=kk;++i)
        {
            int x=rd(),y=rd();
            ht[a[x].b]=1,a[x].m=y;
        }
        sort(a+1,a+n+1);
        int m=n;
        while(m&&ht[a[m].b]) --m;
        int nw=1,la=0;
        g[la][0]=1;
        for(int i=1,sm=0,sb=0;i<=m;++i)
        {
            sm+=a[i].s,sb+=a[i].s;
            if(i<m&&a[i].b==a[i+1].b) continue;
            for(int j=max(sm-sb-c1,0);j<=sm-sb&&j<=c0;++j)
            {
                if(!g[la][j]) continue;
                if(j+sb<=c0) ad(g[nw][j+sb],g[la][j]);
                if(sm-j<=c1) ad(g[nw][j],g[la][j]);
            }
            memset(g[la],0,sizeof(int)*(c0+3));
            sb=0;
            nw^=1,la^=1;
        }
        int lg=la;
        nw=1,la=0;
        h[la][0]=1;
        for(int i=1,sm=0;i<=m;++i)
        {
            sm+=a[i].s;
            for(int j=max(sm-a[i].s-d1,0);j<=sm-a[i].s&&j<=d0;++j)
            {
                if(!h[la][j]) continue;
                if(j+a[i].s<=d0) ad(h[nw][j+a[i].s],h[la][j]);
                if(sm-j<=d1) ad(h[nw][j],h[la][j]);
            }
            memset(h[la],0,sizeof(int)*(d0+3));
            nw^=1,la^=1;
        }
        int lh=la,nh=nw,sk=0,sum=0,ans=0;
        for(int i=1;i<=m;++i) sum+=a[i].s;
        nw=1,la=0;
        f[la][0][0][0]=1;
        for(int i=m+1,sm=0;i<=n;++i)
        {
            bool fg=0;
            sm+=a[i].s;
            if(i==m+1||a[i].b!=a[i-1].b)
            {
                for(int j=0;j<=c0;++j)
                    for(int k=0;k<=sk;++k)
                        f[la][0][j][k]=f[la][1][j][k]=(f[la][0][j][k]+f[la][1][j][k])%mod;
            }
            if(~a[i].m)
            {
                sk+=a[i].s;
                for(int j=max(sm-a[i].s-c1,0);j<=sm-a[i].s&&j<=c0;++j)
                    for(int k=0;k<=sk;++k)
                    {
                        if(f[la][0][j][k])
                        {
                            if(a[i].m!=0&&j+a[i].s<=c0) ad(f[nw][0][j+a[i].s][k+a[i].s],f[la][0][j][k]);
                            if(a[i].m!=1&&j+a[i].s<=c0) ad(f[nw][0][j+a[i].s][k],f[la][0][j][k]);
                        }
                        if(f[la][1][j][k])
                        {
                            if(a[i].m!=2&&sm-j<=c1) ad(f[nw][1][j][k+a[i].s],f[la][1][j][k]);
                            if(a[i].m!=3&&sm-j<=c1) ad(f[nw][1][j][k],f[la][1][j][k]);
                        }
                    }
            }
            else
            {
                for(int j=max(sm-a[i].s-c1,0);j<=sm-a[i].s&&j<=c0;++j)
                    for(int k=0;k<=sk;++k)
                    {
                        if(f[la][0][j][k]&&j+a[i].s<=c0) ad(f[nw][0][j+a[i].s][k],f[la][0][j][k]);
                        if(f[la][1][j][k]&&sm-j<=c1) ad(f[nw][1][j][k],f[la][1][j][k]);
                    }
                for(int j=max(sum+sm-sk-a[i].s-d1,0);j<=sum+sm-sk&&j<=d0;++j)
                {
                    if(!h[lh][j]) continue;
                    if(j+a[i].s<=d0) ad(h[nh][j+a[i].s],h[lh][j]);
                    if(sum+sm-sk-j<=d1) ad(h[nh][j],h[lh][j]);
                }
                fg=1;
            }
            for(int j=0;j<=c0;++j)
                for(int k=0;k<=sk;++k)
                    f[la][0][j][k]=f[la][1][j][k]=0;
            nw^=1,la^=1;
            if(fg)
            {
                memset(h[lh],0,sizeof(int)*(d0+3));
                nh^=1,lh^=1;
            }
        }
        for(int j=1;j<=c0;++j) ad(g[lg][j],g[lg][j-1]);
        for(int j=1;j<=d0;++j) ad(h[lh][j],h[lh][j-1]);
        for(int i=m+1;i<=n;++i) sum+=a[i].s;
        for(int j=0;j<=c0;++j)
            for(int k=0;k<=sk;++k)
            {
                if((f[la][0][j][k]+f[la][1][j][k])%mod==0) continue;
                int l1=max(sum-j-c1,0),r1=min(sum,c0)-j,l2=max(sum-k-d1,0),r2=min(sum,d0)-k;
                if(l1<=r1&&l2<=r2)
                {
                    int gg=(g[lg][r1]-(l1?g[lg][l1-1]:0)+mod)%mod,hh=(h[lh][r2]-(l2?h[lh][l2-1]:0)+mod)%mod;
                    ans=(ans+1ll*(f[la][0][j][k]+f[la][1][j][k])%mod*gg%mod*hh%mod)%mod;
                }
            }
        printf("%d\n",ans);
    }
    return 0;
}

原文地址：https://www.cnblogs.com/smyjr/p/10713772.html

时间： 2024-11-09 07:50:18

luogu P5289 [十二省联考2019]皮配的相关文章

[十二省联考2019]皮配

题目链接 选一个派系和一个阵营可以唯一确定一名导师 因为每一个阵营里的导师都分别来自不同派系,所以k=0时,对阵营的选择是不影响对派系的选择的 唯一的限制就是同城市的要在同一个阵营 所以以每个城市为物品,物品大小为该城市的人数,阵营人数为背包容量,做背包dp 再以每个学校为物品,物品大小为该学校的人数,派系人数为背包容量,做背包dp 只用一维记录背包大小即可,因为总人数-背包里的人数=在另一个阵营或派系的人数 然后合并答案即可 方案数是可以相互乘起来的,k很小,所以我们可以暴力做k!=0的情况,

luogu P5284 [十二省联考2019]字符串问题

传送门 如果某个串\(b\)是\(a_j\)的前缀,并且\(a_i\)支配\(b\),那么\(a_i\)后面就可以放\(a_j\),所以如果把对应的图建出来,问题就是求最长链,如果有环就是无限长 说到前缀,我们可以把所有\(a\)串建一棵\(Trie\),然后某个\(a_i\)支配的\(b\)串节点对应的子树内的\(a_j\)串节点都可以连边\((i,j)\).不过每次建\(Trie\)太慢了,我们可以利用\(SAM\)建后缀树,然后每次相当于一个点给某个点子树内的点连边,所以后缀树上父亲向儿子

luogu P5285 [十二省联考2019]骗分过样例

传送门 你管这玩意叫传统题? 1_998244353通过样例可以发现就是\(19^a\mod\ 998244353\),然后如果读入的数过大,根据费马小定理,那个读进来的数对\(998244352\)取模就行了 1?和1?+就是\(19^a\mod\ ?\),前者可以暴力求出,后者比较牛皮,可以找到数据中的两个数\(x,y(x<y)\),并且\(y-x\)较小,然后可以知道\(19^x,19^{y-x},19^y\),因为\(19^x*19^{y-x}\equiv 19^y\mod\ ?\),然

【题解】Luogu P5283 [十二省联考2019]异或粽子

原题传送门 看见一段的异或和不难想到要做异或前缀和\(s\) 我们便将问题转化成:给定\(n\)个数,求异或值最靠前的\(k\)对之和 我们珂以建一个可持久化01trie,这样我们就珂以求出每个值\(s[a]\)与之前所有的值异或值最大的值\(b\)是多少,把这些所有\((b,a)\)塞进一个堆中 每次从堆顶取元素,设这个元素为\((b,a)\),要将\(b\)加入答案,并且在版本\(a\)的01trie中减去\(s[a]\)^\(b\),再取出\(s[a]\)与01trie中的数异或最大值(原

[十二省联考2019]异或粽子（可持久化tire，堆）

[十二省联考2019]异或粽子(luogu) Description 题目描述 小粽是一个喜欢吃粽子的好孩子.今天她在家里自己做起了粽子. 小粽面前有 nn 种互不相同的粽子馅儿,小粽将它们摆放为了一排,并从左至右编号为 11 到 nn.第 ii 种馅儿具有一个非负整数的属性值 a_iai?.每种馅儿的数量都足够多,即小粽不会因为缺少原料而做不出想要的粽子.小粽准备用这些馅儿来做出 kk 个粽子. 小粽的做法是:选两个整数数 ll, rr,满足 1 \leqslant l \leqslant r

「十二省联考 2019」字符串问题

「十二省联考 2019」字符串问题 解题思路 傻逼题.. 考虑问题转化为一个A串向其支配的所有B串的后缀A串连边,如果有环答案 \(-1\) 否则是这个 \(\text{DAG}\) 上最长路径,直接建图是 \(n^2\) 的,考虑优化建图即可. 由于 \(A,B\) 都是原串的一个子串,那么对原串的反串建 SAM,一个子串的后缀就是其所在节点上比它长的串以及,其子树里的所有串. 首先将所有 \(A,B\) 串在 SAM上用倍增定位并新建节点,把SAM上每个节点拆成入点和出点,对于SAM每一个节

「十二省联考 2019」字符串问题 解题报告

「十二省联考 2019」字符串问题 当场就去世了,我这菜人改了一下午 考虑一个A,B之间的连边实际表示了两个A之间的有向边,然后把A的连边处理好,就转成了拓扑排序找环+最长链 但是边数很多,考虑优化连边 A,B之间的连边显然没法优化的,考虑一个B可以表示所有它的后缀A 把串反向建出SAM,然后一个B的后缀就是par树的子树 可以拿倍增定位 好了这题就没了 注意到一个事情,定位的点可能重复,于是对SAM拆点,每个点挂一个vector表示一个A或者B的点在SAM的这个位置 然后考虑如何连边 一个B所

[十二省联考2019] 异或粽子 解题报告 (可持久化Trie+堆)

interlinkage: https://www.luogu.org/problemnew/show/P5283 description: solution: 显然有$O(n^2)$的做法,前缀和优化一下即可 正解做法是先确定一个右端点$r$,找到最优的$l$使得该区间的异或和最大,这个可以用可持久化$Trie$实现.不懂的话可以在我的博客里搜索 对每个点取出来后把答案放进一个堆里,显然当前的堆顶一定会对答案产生贡献 然后我们考虑每次取出的右端点,它依旧可能产生贡献.即上一次取的最优的$l$把

P5283 [十二省联考2019]异或粽子

传送门 超级钢琴+可持久化$Trie$ 同样设三元组 $(o,l,r)$ 表示左端点为 $o$,右端点 $\in [l,r]$ 的区间的最大异或值,这个东西可以用可持久化 $Trie$ 来维护 一开始把所有 $(i,i,n)$ 扔到堆里,然后每次取出计算贡献,设取得最大异或值的位置为 $t$,然后再把 $(o,l,t-1)$ 和 $(o,t+1,r)$ 扔到堆里 具体还是看代码,很容易理解 注意可能爆 $int$,所以要开 $unsigned\ int$,要注意代码常数,我代码 $luogu$