咕咕咕(gugugu)
【题目描述】
有一个长度为N的数轴,数轴的范围从 1到 N 还有 M个区间,其中第i个区间是 [Li,Ri]并且有一个权值 Vi。定义一个选择区间的方案是合法的,当且仅当数轴上每一个整点都被至少一个区间覆盖了,同时该方案的权值是所有选择区间的权值乘积。
请你求出所有合法方案权值的和。由于答案可能很大,请输出其对1e9+7取模的结果。
【输入格式】
第一行两个正整数 N,M 。
接下来 M行,第 i 行先是两个正整数 Li,Ri描述区间,然后一个非负整数 Vi描述区间的权值。
【输出格式】
输出仅一行,表示答案对 1e9+7取模的结果。
【数据范围】
1<=N,M<=2e5+5,Vi<1e9+7
蒟蒻分析:首先来个排序就nlogn,然后在On扫一遍得到答案,结果发现扫不来呀!
蒟蒻总结:思维就陷进了出题人的圈套里,总是想怎么样组成方案,组成方案后答案怎么统计,结果正解根本不是这么想的,其实当时就明白了这样想是死路一条,但就是偏偏往里钻
solution by jklover %%
将所有区间按照左端点 Li从小到大排序后,
设 f(i,j) 表示考虑了前 i 个区间,所有覆盖了 [1,j] 的方案的权值和.
转移时,若不选第 i 个区间,则
f(i,j)+=f(i−1,j)
若选第 i个区间,则只能从 f(i−1,j)转移过来,否则
f(i,max(j,Ri))+=vi⋅f(i−1,j),j≥Li−1
利用线段树优化转移,时间复杂度 O(mlogn).
然后我搞忘区间乘法怎么pushdown了(先说下面这pushdown不是我写的,是不是觉得超级丑?)
void pushdown(int root, int l, int r){ int m=(l+r)/2; //根据我们规定的优先度,儿子的值=此刻儿子的值*爸爸的乘法lazytag+儿子的区间长度*爸爸的加法lazytag st[root*2].v=(st[root*2].v*st[root].mul+st[root].add*(m-l+1))%p; st[root*2+1].v=(st[root*2+1].v*st[root].mul+st[root].add*(r-m))%p; //很好维护的lazytag st[root*2].mul=(st[root*2].mul*st[root].mul)%p; st[root*2+1].mul=(st[root*2+1].mul*st[root].mul)%p; st[root*2].add=(st[root*2].add*st[root].mul+st[root].add)%p; st[root*2+1].add=(st[root*2+1].add*st[root].mul+st[root].add)%p; //把父节点的值初始化 st[root].mul=1; st[root].add=0; return ; }
code by theshaow
#include<bits/stdc++.h> #define del(a,i) memset(a,i,sizeof(a)) #define ll long long #define inl inline #define il inl void #define it inl int #define ill inl ll #define re register #define ri re int #define rl re ll #define INF 0x3f3f3f3f #define lowbit(x) (x&(-x)) #define mid ((l+r)>>1) using namespace std; namespace io { const int SIZE = (1 << 21) + 1; char ibuf[SIZE], *iS, *iT, obuf[SIZE], *oS = obuf, *oT = oS + SIZE - 1, c, qu[55]; int f, qr; #define gc() (iS == iT ? (iT = (iS = ibuf) + fread (ibuf, 1, SIZE, stdin), (iS == iT ? EOF : *iS ++)) : *iS ++) inline void flush () { fwrite (obuf, 1, oS - obuf, stdout); oS = obuf; } inline void putc (char x) { *oS ++ = x; if (oS == oT) flush (); } template <class I> inline void read (I &x) { for (f = 1, c = gc(); c < '0' || c > '9'; c = gc()) if (c == '-') f = -1; for (x = 0; c <= '9' && c >= '0'; c = gc()) x = x * 10 + (c & 15); x *= f; } template <class I> inline void print (I x) { if (!x) putc ('0'); if (x < 0) putc ('-'), x = -x; while (x) qu[++ qr] = x % 10 + '0', x /= 10; while (qr) putc (qu[qr --]); } struct Flusher_ {~Flusher_(){flush();}}io_flusher_; } using io :: read; using io :: putc; using io :: print; const int MAXN = 2e5+5,mod = 1e9+7; int n,m,dp[MAXN],a[MAXN]; struct Node{ int l,r,val; bool operator <(const Node &t) const{ return l==t.l?r<t.r:l<t.l; } }node[MAXN]; inl bool cmp(Node x,Node y){ return x.l==y.l?x.r<y.r:x.l<y.l; } it add(int x,int y){return x+y>=mod?x+y-mod:x+y;} it mul(int x,int y){return 1ll*x*y%mod;} #define lc (cur<<1) #define rc (cur<<1|1) struct Seg_Tree{ int sum,tag_add,tag_mul; }T[MAXN<<2]; il build(int cur,int l,int r){ T[cur].tag_mul=1; if(l==r) return ; build(lc,l,mid),build(rc,mid+1,r); } il pushup(int cur){T[cur].sum=add(T[lc].sum,T[rc].sum);} il pushdown(int cur,int l,int r){ if(!T[cur].tag_add&&T[cur].tag_mul==1) return ; if(T[cur].tag_mul!=1){ T[lc].sum=mul(T[lc].sum,T[cur].tag_mul); T[lc].tag_mul=mul(T[lc].tag_mul,T[cur].tag_mul); T[lc].tag_add=mul(T[lc].tag_add,T[cur].tag_mul); T[rc].sum=mul(T[rc].sum,T[cur].tag_mul); T[rc].tag_mul=mul(T[rc].tag_mul,T[cur].tag_mul); T[rc].tag_add=mul(T[rc].tag_add,T[cur].tag_mul); T[cur].tag_mul=1; } if(T[cur].tag_add){ T[lc].sum=add(T[lc].sum,mul(T[cur].tag_add,mid-l+1)); T[lc].tag_add=add(T[lc].tag_add,T[cur].tag_add); T[rc].sum=add(T[rc].sum,mul(T[cur].tag_add,r-mid)); T[rc].tag_add=add(T[rc].tag_add,T[cur].tag_add); T[cur].tag_add=0; } } il updata_add(int cur,int l,int r,int pos,int k){ if(l==r) return T[cur].tag_add=add(T[cur].tag_add,k),T[cur].sum=add(T[cur].sum,k),void(); pushdown(cur,l,r); if(mid>=pos) updata_add(lc,l,mid,pos,k); else updata_add(rc,mid+1,r,pos,k); pushup(cur); } il updata_mul(int cur,int l,int r,int L,int R,int k){ if(l>=L&&r<=R) return T[cur].sum=mul(T[cur].sum,k),T[cur].tag_add=mul(T[cur].tag_add,k),T[cur].tag_mul=mul(T[cur].tag_mul,k),void(); pushdown(cur,l,r); if(mid>=L) updata_mul(lc,l,mid,L,R,k); if(R>mid) updata_mul(rc,mid+1,r,L,R,k); pushup(cur); } it query(int cur,int l,int r,int L,int R){ if(l>=L&&r<=R) return T[cur].sum; pushdown(cur,l,r);ri res=0; if(mid>=L) res=add(res,query(lc,l,mid,L,R)); if(R>mid) res=add(res,query(rc,mid+1,r,L,R)); return res; } int main(){ freopen("gugugu.in","r",stdin); freopen("gugugu.out","w",stdout); read(n),read(m); for(ri i=1;i<=m;++i) read(node[i].l),read(node[i].r),read(node[i].val); sort(node+1,node+1+m,cmp); if(node[1].l!=1) return puts("0"),0; build(1,0,n),updata_add(1,0,n,0,1); for(ri i=1;i<=m;++i){ ri res=query(1,0,n,node[i].l-1,node[i].r); updata_add(1,0,n,node[i].r,mul(res,node[i].val)); updata_mul(1,0,n,node[i].r+1,n,node[i].val+1); } print(query(1,0,n,n,n)); return 0; }
第二题gugu太简单了,
送分题,就是出题人比较傻逼,没说条件能否成立,样例又没这种情况
相信大家只要注意到了,这题都会A掉,不然就只有40分哦

直接solution by jklover%%%%
gugu

code by theshadow
#include<bits/stdc++.h> #define del(a,i) memset(a,i,sizeof(a)) #define ll long long #define inl inline #define il inl void #define it inl int #define ill inl ll #define re register #define ri re int #define rl re ll #define INF 0x3f3f3f3f #define lowbit(x) (x&(-x)) #define mid ((l+r)>>1) using namespace std; namespace io { const int SIZE = (1 << 21) + 1; char ibuf[SIZE], *iS, *iT, obuf[SIZE], *oS = obuf, *oT = oS + SIZE - 1, c, qu[55]; int f, qr; #define gc() (iS == iT ? (iT = (iS = ibuf) + fread (ibuf, 1, SIZE, stdin), (iS == iT ? EOF : *iS ++)) : *iS ++) inline void flush () { fwrite (obuf, 1, oS - obuf, stdout); oS = obuf; } inline void putc (char x) { *oS ++ = x; if (oS == oT) flush (); } template <class I> inline void read (I &x) { for (f = 1, c = gc(); c < '0' || c > '9'; c = gc()) if (c == '-') f = -1; for (x = 0; c <= '9' && c >= '0'; c = gc()) x = x * 10 + (c & 15); x *= f; } template <class I> inline void print (I x) { if (!x) putc ('0'); if (x < 0) putc ('-'), x = -x; while (x) qu[++ qr] = x % 10 + '0', x /= 10; while (qr) putc (qu[qr --]); } struct Flusher_ {~Flusher_(){flush();}}io_flusher_; } using io :: read; using io :: putc; using io :: print; const int MAXN = 3e3+5,mod = 1e9+7; int n,m,k,u,v,a[MAXN],b[MAXN],s[MAXN*2],sz,f1[MAXN][MAXN],f2[MAXN],ans; char lt[MAXN<<1][MAXN<<1]; it add(int x,int y){return x+y>=mod?x+y-mod:x+y;} int main(){ freopen("gu.in","r",stdin); freopen("gu.out","w",stdout); read(n),read(m),read(k); for(ri i=1;i<=n;++i) read(a[i]),s[++sz]=a[i]; for(ri i=1;i<=m;++i) read(b[i]),s[++sz]=b[i]; sort(s+1,s+1+sz),sz=unique(s+1,s+1+sz)-s-1; for(ri i=1;i<=n;++i) a[i]=lower_bound(s+1,s+1+sz,a[i])-s; for(ri i=1;i<=n;++i) b[i]=lower_bound(s+1,s+1+sz,b[i])-s; for(ri i=1;i<=k;++i){ read(u),read(v); ri x=lower_bound(s+1,s+1+sz,u)-s; ri y=lower_bound(s+1,s+1+sz,v)-s; if(s[x]==u&&s[y]==v) lt[x][y]=1; } for(ri i=1;i<=n;++i){ ri sum=0; for(ri j=1;j<=m;++j){ if(a[i]==b[j]) f1[i][j]=add(sum,1),ans=add(ans,f1[i][j]); if(lt[b[j]][a[i]]) sum=add(sum,f2[j]); f2[j]=add(f2[j],f1[i][j]); } } print(ans); return 0; }