Codeforces round 1106
Div 2 536
我还是太菜了.jpg
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的更多相关文章
- 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 ...
- Codeforces Round #354 (Div. 2) ABCD
Codeforces Round #354 (Div. 2) Problems # Name A Nicholas and Permutation standard input/out ...
- Codeforces Round #368 (Div. 2)
直达–>Codeforces Round #368 (Div. 2) A Brain’s Photos 给你一个NxM的矩阵,一个字母代表一种颜色,如果有”C”,”M”,”Y”三种中任意一种就输 ...
- cf之路,1,Codeforces Round #345 (Div. 2)
cf之路,1,Codeforces Round #345 (Div. 2) ps:昨天第一次参加cf比赛,比赛之前为了熟悉下cf比赛题目的难度.所以做了round#345连试试水的深浅..... ...
- Codeforces Round #279 (Div. 2) ABCDE
Codeforces Round #279 (Div. 2) 做得我都变绿了! Problems # Name A Team Olympiad standard input/outpu ...
- 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 ...
- 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 ...
- Codeforces Round #370 - #379 (Div. 2)
题意: 思路: Codeforces Round #370(Solved: 4 out of 5) A - Memory and Crow 题意:有一个序列,然后对每一个进行ai = bi - bi ...
- Codeforces Round #371 (Div. 1)
A: 题目大意: 在一个multiset中要求支持3种操作: 1.增加一个数 2.删去一个数 3.给出一个01序列,问multiset中有多少这样的数,把它的十进制表示中的奇数改成1,偶数改成0后和给 ...
随机推荐
- 走通Django的基本流程
工程目录及文件的说明 manage.py:一个命令行工具,可以使我们用多种方式对Django项目进行交互 __init__.py:一个空文件,它告诉Python这个文件的上级目录应该看做一个pytho ...
- JavaScript 事件处理机制
DOM 事件流 DOM 是一个树型结构,当一个HTML元素产生一个事件时,该事件会在该元素结点与根节点之间按特定的顺序传播,路径所经过的节点都会收到该事件,这个传播过程可称为DOM事件流.而事件传播的 ...
- Intellij idea Cannot start internal HTTP server.
错误提示:Cannot start internal HTTP server. Git integration, JavaScript debugger and LiveEdit may operat ...
- 安装nvm管理不同的node版本
在工作或者学习中,偶尔会遇到需要切换不同node版本的需求,幸好有神器nvm可以帮我们解决问题.下面我们就来讲解如何在window系统上安装nvm!
- PHP编译安装时常见错误解决办法,php编译常见错误
PHP编译安装时常见错误解决办法,php编译常见错误 1.configure: error: xslt-config not found. Please reinstall the libxslt & ...
- 葡萄城报表V11 SP2新版本震撼发布!
葡萄城报表正式发布 v11.2 版本,强势推出国内首创的基于HTML5的在线报表设计器,从此报表设计告别桌面应用程序,随时随地修改报表,真正跨平台操作,从而使任何报表用户更快速的响应报表业务变化! 在 ...
- SurfaceViewVideoList网络获取视频播放
主布局: <?xml version="1.0" encoding="utf-8"?> <RelativeLayout xmlns:andro ...
- 【SPL标准库专题(1)】 SPL简介
什么是SPL SPL是Standard PHP Library(PHP标准库)的缩写. 根据官方定义,它是"a collection of interfaces and classes th ...
- 转:sqlserver 临时表、表变量、CTE的比较
1.临时表 1.1 临时表包括:以#开头的局部临时表,以##开头的全局临时表. 1.2 存储 不管是局部临时表,还是全局临时表,都会放存在tempdb数据库中. 1.3 作用域 局部临时表:对当前连接 ...
- 重置AD用户密码
$cc = import-csv D:\Powershell\Tauba.csv foreach ($c in $cc) { $c.username $pwd = $c.password get-ad ...