牛客网NOIP赛前集训营-提高组（第六场）B-选择题[背包]

题意

题目链接

分析

直接背包之后可以 $O(n)$ 去除一个物品的影响。
注意特判 $[p==1]$ 的情况。
总时间复杂度为 $O(n^2)$ 。

代码

#include<bits/stdc++.h>
using namespace std;
#define go(u) for(int i=head[u],v=e[i].to;i;i=e[i].last,v=e[i].to)
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define pb push_back
typedef long long LL;
inline int gi(){int x=0,f=1;char ch=getchar();while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-48;ch=getchar();}return x*f;
}
template<typename T>inline bool Max(T &a,T b){return a<b?a=b,1:0;}
template<typename T>inline bool Min(T &a,T b){return b<a?a=b,1:0;}
const int N=2004,mod=998244353;
int n;
LL p[4][N],w[4][4],f[N],tmp[N],ans[N];
LL Pow(LL a,LL b){LL res=1;for(;b;b>>=1,a=1ll*a*a%mod) if(b&1) res=res*a%mod;return res;
}
void solve(int k){memset(f,0,sizeof f);f[0]=1;rep(i,1,n){LL p=::p[k][i];for(int j=n;~j;--j)f[j]=((j-1>=0?f[j-1]:0)*p+f[j]*(1-p))%mod;}rep(i,1,n){LL p=::p[k][i];LL inv=Pow(1-p,mod-2);if(p^1){tmp[0]=f[0]*inv%mod;rep(j,1,n) tmp[j]=(f[j]-.tmp[j-1]*p)%mod*inv%mod;}elserep(j,0,n-1) tmp[j]=f[j+1];LL sum=0;for(int j=n/2+1;j<n;++j) (sum+=tmp[j])%=mod;rep(a,0,3)(ans[i]+=(sum+(k==a)*tmp[n/2])*::p[a][i]%mod*w[k][a]%mod)%=mod;}
}
int main(){n=gi();rep(i,1,n) rep(j,0,3) p[j][i]=gi();rep(i,0,3) rep(j,0,3) w[i][j]=gi();rep(k,0,3) solve(k);rep(i,1,n) printf("%lld\n",(ans[i]+mod)%mod);return 0;
}

转载于:https://www.cnblogs.com/yqgAKIOI/p/9829494.html

牛客网NOIP赛前集训营-提高组（第六场）B-选择题[背包]相关推荐

牛客网NOIP赛前集训营提高组第5场 T2 旅游
[题解] 我们可以发现不在最小生成树上的边一定不能多次经过,因为一条不在最小生成树上的边(u,v)的边权比最小生成树上(u,v)之间的路径更长,选择不在最小生成树上的边一定不划算. 我们还需要确定最小 ...
牛客网NOIP赛前集训营提高组（第七场）
中国式家长 2 链接:https://www.nowcoder.com/acm/contest/179/A 来源:牛客网时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K ...
牛客网NOIP赛前集训营-提高组（第六场）B-选择题
题目描述有一道选择题,有 a,b,c,d 四个选项. 现在有 n 个人来做这题,第 i 个人有 pi,j 的概率选第 j 个选项. 定义$cnt(x)$为选第$ x $个选项的人数. 令\(mx ...
[NowCoder]牛客网NOIP赛前集训营-提高组（第六场）题解
A.最长路题意:给定有向图,每条边有个字符$[0,10^9]$,求每个点最长路字典序最小的方案.$N,M\le 10^6$ 建反图跑拓扑排序,显然入过队的点都有最长路,考虑如何判断字典序大小 ...
牛客网NOIP赛前集训营-提高组（第七场）C 洞穴
洞穴思路: 矩阵乘法由于只需要知道 A ^ l 的第a行第b个元素所以我们每次在做矩阵乘法时只需要算第a行就可以了还要像矩阵快速幂一样预处理A ^ (1<<d) 代码: #pra ...
牛客网NOIP赛前集训营-提高组（第七场）C-洞穴
题目描述有一天,牛牛找到了一个巨大的洞穴.洞穴可以描述成一个有向图,一共有$N$个节点(从$1$到$N$编号)和$M$条长度为$1$的有向边,每条边从某一个节点$u$连向另一 ...
【题解】[牛客网NOIP赛前集训营-提高组（第七场）]C.洞穴倍增优化DP+bitset
题目链接 #include<cstdio> #include<bitset> using namespace std; const int N=110; int n,m,q; ...
牛客网NOIP赛前集训营-普及组（第一场）
前三题略 T4: 题目描述小A有n个长度都是L的字符串.这些字符串只包含前8个小写字符,'a'~'h'.但这些字符串非常的混乱,它们几乎长得互不相同.小A想通过一些规则,让它们长得尽可能相同.小A现 ...
Nowcoder | [题解-N165]牛客网NOIP赛前集训营-普及组（第二场）
啊...表示一大早还没睡醒就开始打比赛(开始前一分钟的我还在桌子上趴着休眠)...表示题目思路清奇(尤其C题)...但是我还是太蒻了...$D$题暴力都没打...题解正式开始之前先$\%$一下 ...

牛客网NOIP赛前集训营-提高组（第六场）B-选择题[背包]

题意

分析

代码

牛客网NOIP赛前集训营-提高组（第六场）B-选择题[背包]相关推荐

最新文章

热门文章