码迷,mamicode.com
首页 > 其他好文 > 详细

绵阳东辰国际test.10.7eve

时间:2019-10-08 12:20:17      阅读:101      评论:0      收藏:0      [点我收藏+]

标签:cpp   爸爸   namespace   har   void   pushd   amp   sizeof   格式   

咕咕咕(gugugu)

【题目描述】

有一个长度为N的数轴,数轴的范围从 1N 还有 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(mlog?n).

然后我搞忘区间乘法怎么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;
}

绵阳东辰国际test.10.7eve

标签:cpp   爸爸   namespace   har   void   pushd   amp   sizeof   格式   

原文地址:https://www.cnblogs.com/wzxbeliever/p/11634493.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!