Div 2 536

E

傻逼DP直接做

我居然调了1.5h

我真的是太菜了.jpg

堆+扫描线直接维护每个位置的贪心结果

然后要么使用干扰

要么就接受贪心的结果

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <iostream>
#include <bitset>
using namespace std;
#define N 100005
#define M 205
#define ll long long
int n,m,K,p[N];ll f[N][M];
struct node{int s,t,d,w;}a[N];
bool cmp(const node &a,const node &b){return a.s==b.s?a.t<b.t:a.s<b.s;}
priority_queue<pair<pair<int ,int > ,int > >q;
int main()
{
scanf("%d%d%d",&n,&m,&K);memset(f,0x3f,sizeof(f));
for(int i=1;i<=K;i++)scanf("%d%d%d%d",&a[i].s,&a[i].t,&a[i].d,&a[i].w);sort(a+1,a+K+1,cmp);
for(int i=1,j=1;i<=n;i++)
{
while(a[j].s<=i&&j<=K)q.push(make_pair(make_pair(a[j].w,a[j].d),j)),j++;
while(!q.empty()&&a[q.top().second].t<i)q.pop();
if(!q.empty())p[i]=q.top().second;
}
f[1][0]=0;
for(int i=1;i<=n;i++)
for(int j=0;j<=m;j++)
{
int x=p[i];
f[i+1][j+1]=min(f[i][j],f[i+1][j+1]);
if(x)f[a[x].d+1][j]=min(f[i][j]+a[x].w,f[a[x].d+1][j]);
else f[i+1][j]=min(f[i][j],f[i+1][j]);
}
ll ans=1ll<<60;
for(int i=0;i<=m;i++)ans=min(ans,f[n+1][i]);
printf("%lld\n",ans);
}

F

  • 这不是一个特征多项式优化常系数线性齐次递推裸题嘛!

然后我就开始写了...

然后我发现我不会求$K$次剩余...

然后我就GG了...

所以这个题不用会求$K$次剩余...

那么根据原根的性质,我们可以发现,$K$次剩余可以表达为$\frac{q}{p} \mod 998244352$的形式,其中$q$表达为$m = 3^q \mod 998244353$,$p=k$

如果$K$不存在逆元的话,就没有对应的$K$次剩余...

然后,就可以通过BSGS+exgcd求

所以矩阵乘法就能做的题,为啥我要用特征多项式啊!!!!!

附上代码:

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <bitset>
#include <map>
using namespace std;
#define N 205
#define ll long long
#define mod 998244353
#define mmod 998244352
int a[N],b[N],mo[N],tmp[N],n,k,m,ans;
int q_pow(int x,int n){int ret=1;for(;n;n>>=1,x=(ll)x*x%mod)if(n&1)ret=(ll)ret*x%mod;return ret;}
void mul(int *a,int *b,int *ret)
{
memset(tmp,0,sizeof(tmp));
for(int i=0;i<k;i++)
for(int j=0;j<k;j++)
tmp[i+j]=(tmp[i+j]+(ll)a[i]*b[j])%mmod;
for(int i=(k<<1)-2;i>=k;i--)if(tmp[i])for(int j=1;j<=k;j++)
tmp[i-j]=(tmp[i-j]-(ll)tmp[i]*mo[k-j])%mmod;
for(int i=0;i<k;i++)ret[i]=tmp[i];
}
map<int ,int >mp;
int ex_gcd(int a,int b,int &x,int &y)
{
if(!b)return x=1,y=0,a;int ret=ex_gcd(b,a%b,y,x);
y=y-a/b*x;return ret;
}
int BSGS(int x)
{
int s=5000,t=1;
for(int i=0;i<s;i++)mp[int((ll)t*x%mod)]=i,t=t*3ll%mod;
for(int i=1,now=1;;i++)
{
now=(ll)now*t%mod;
if(mp.find(now)!=mp.end())return i*s-mp[now];
}
}
int main()
{
scanf("%d",&k);
for(int i=1;i<=k;i++)scanf("%d",&mo[k-i]),mo[k-i]=mmod-mo[k-i];
scanf("%d%d",&n,&m);b[1]=1;a[0]=1;
if(k==1)b[0]=mmod-mo[0],b[1]=0;
for(n--;n;n>>=1){if(n&1)mul(a,b,a);mul(b,b,b);}
int t=(a[k-1]+mmod)%mmod,x=0,y=0;
int tt=ex_gcd(t,mmod,x,y),kk=BSGS(m);
if(kk%tt)return puts("-1"),0;
x=((ll)x*(kk/tt)%mmod+mmod)%mmod;
printf("%d\n",q_pow(3,x));
}

Codeforces round 1106的更多相关文章

  1. Codeforces Round #366 (Div. 2) ABC

    Codeforces Round #366 (Div. 2) A I hate that I love that I hate it水题 #I hate that I love that I hate ...

  2. Codeforces Round #354 (Div. 2) ABCD

    Codeforces Round #354 (Div. 2) Problems     # Name     A Nicholas and Permutation standard input/out ...

  3. Codeforces Round #368 (Div. 2)

    直达–>Codeforces Round #368 (Div. 2) A Brain’s Photos 给你一个NxM的矩阵,一个字母代表一种颜色,如果有”C”,”M”,”Y”三种中任意一种就输 ...

  4. cf之路,1,Codeforces Round #345 (Div. 2)

     cf之路,1,Codeforces Round #345 (Div. 2) ps:昨天第一次参加cf比赛,比赛之前为了熟悉下cf比赛题目的难度.所以做了round#345连试试水的深浅.....   ...

  5. Codeforces Round #279 (Div. 2) ABCDE

    Codeforces Round #279 (Div. 2) 做得我都变绿了! Problems     # Name     A Team Olympiad standard input/outpu ...

  6. Codeforces Round #262 (Div. 2) 1003

    Codeforces Round #262 (Div. 2) 1003 C. Present time limit per test 2 seconds memory limit per test 2 ...

  7. Codeforces Round #262 (Div. 2) 1004

    Codeforces Round #262 (Div. 2) 1004 D. Little Victor and Set time limit per test 1 second memory lim ...

  8. Codeforces Round #370 - #379 (Div. 2)

    题意: 思路: Codeforces Round #370(Solved: 4 out of 5) A - Memory and Crow 题意:有一个序列,然后对每一个进行ai = bi - bi  ...

  9. Codeforces Round #371 (Div. 1)

    A: 题目大意: 在一个multiset中要求支持3种操作: 1.增加一个数 2.删去一个数 3.给出一个01序列,问multiset中有多少这样的数,把它的十进制表示中的奇数改成1,偶数改成0后和给 ...

随机推荐

  1. element-ui Tag、Dialog组件源码分析整理笔记(五)

    Tag 标签组件 <script> export default { name: 'ElTag', props: { text: String, closable: Boolean, // ...

  2. CentOS7+ApacheServer2.4+MariaDB10.0+PHP7.0+禅道项目管理软件8.0环境部署

    CentOS7+ApacheServer2.4+MariaDB10.0+PHP7.0+禅道项目管理软件8.0环境部署 by:授客 QQ:1033553122 目录 一. 二. 三. 四. 五. 六. ...

  3. Android项目实战(三十):Fresco加载gif图片并播放

    前言: 项目中图文混合使用的太多太多了,但是绝大部分都是静态图片. 然而项目开发中有这么一个需求:显示一个出一个简短的动画(一般都不超过3秒)演示 比如说:一个功能提供很多步骤来教用户做广播体操,那么 ...

  4. Android--WebView 自适应代码

    //WebView自适应代码 private String getHtmlData(String bodyHTML) { String head = "<head>" ...

  5. 《ECMAScript6标准入门》第三版--读书笔记

    2015年6月,ECMAScript 6正式通过,成为国际标准.尽管在目前的工作中还没有使用ES6,但是每项新技术出来总是忍不住想尝尝鲜,想知道ES6能为前端开发带来哪些变化?对自己的工作有哪些方面可 ...

  6. Supervisor 管理进程,Cloud Insight 监控进程,完美!

    Supervisor 是由 Python 语言编写.基于 linux 操作系统的一款服务器管理工具,用于监控服务器的运行,发现问题能立即自动预警及自动重启等. Cloud Insight 是一款次世代 ...

  7. 从 Azure 下载 Windows VHD

    本文介绍如何使用 Azure 门户从 Azure 下载 Windows 虚拟硬盘 (VHD) 文件. Azure 中的虚拟机 (VM) 将磁盘用作存储操作系统.应用程序和数据的位置. 所有 Azure ...

  8. ejb-jar.xml

    所有bean类(无论是会话bean还是实体bean)必须实现的最基本的接口是javax.ejb.EnterpriseBean接口. 所有的会话bean必须实现javax.ejb.SessionBean ...

  9. etc/skel目录介绍

    /etc/skel目录的作用: /etc/skel目录是用来存放新用户配置文件的目录,当我们添加新用户时,这个目录下的所有文件会自动被复制到新添加的用户家目录下,默认情况下,/etc/skel 目录下 ...

  10. 【转】Linux 高级的视角来查看Linux引导过程

    [原文]https://www.toutiao.com/i6594210975480545800/ 1.概述 图 1 是我们在20,000 英尺的高度看到的视图. 当系统首次引导时,或系统被重置时,处 ...