BZOJ一句话

 

一句话题解集合。

1061: [Noi2008]志愿者招募

 单纯形，运用对偶原理转化过来，变成标准形然后单纯性裸上即可。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstring>
#include<algorithm>
const double eps=1e-6;
const double oo =1e9;
double a[10001][1001];
double b[10001];
double c[1001];
int n,m;
int read(void)
{
    int ans=0;
    int f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')f=-f;
        ch=getchar();
    }
    while(ch<='9'&&ch>='0')
    {
        ans=ans*10+ch-'0';
        ch=getchar();
    }
    return ans*f;
}
void pivot(double&ans,int l,int e)
{
    b[l]/=a[l][e];
    for(int i=1;i<=n;i++)if(i!=e)a[l][i]/=a[l][e];
    a[l][e]=1.00/a[l][e];
    for(int i=1;i<=m;i++)
    if(fabs(a[i][e])>eps&&i!=l)
    {
        b[i]-=b[l]*a[i][e];
        for(int j=1;j<=n;j++)if(j!=e)
            a[i][j]-=a[i][e]*a[l][j];
        a[i][e]=-a[i][e]*a[l][e];
    }
    ans+=c[e]*b[l];
    for(int i=1;i<=n;i++)if(i!=e)c[i]-=c[e]*a[l][i];
    c[e]=-c[e]*a[l][e];
    return ;
}
double simplex(void)
{
    double ans=0;
    while(true)
    {
        int e,l;
        for(e=1;e<=n;e++)if(c[e]>eps)break;
        if(e>n)return ans;
        double mins=oo;
        for(int i=1;i<=m;i++)
            if(a[i][e]>eps&&mins>b[i]/a[i][e])
                mins=b[i]/a[i][e],l=i;
        if(mins>=oo-eps)return oo;
        pivot(ans,l,e);
    }
    return ans;
}
int main()
{
//    freopen("a.in","r",stdin);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)c[i]=read();
    for(int i=1;i<=m;i++)
    {
        int s,t;
        scanf("%d%d",&s,&t);
        for(int j=s;j<=t;j++)
        {
            if(j<1)j=1;if(j>n)break;
            a[i][j]=1;
        }
        b[i]=read();
    }
    printf("%d\n",(int)(simplex()+0.5));
    return 0;
}

1197: [HNOI2006]花仙子的魔法

神TM递推题，考虑升维的计算仍然是类似的，就有递推式$f_{[i][j]}=f_{[i-1][j-1]}+f_{[i][j-1]}$

#include<cstdio>
typedef long long lnt;
lnt dp[17][1000];
int n,m;
int main()
{
//    freopen("flower.in","r",stdin);
//    freopen("flower.out","w",stdout);
    scanf("%d%d",&m,&n);dp[1][0]=1;
    for(int i=1;i<=m;i++)dp[1][i]=i*2;
    for(int i=2;i<=n;i++)
    {
        dp[i][0]=1;
        for(int j=1;j<=m;j++)dp[i][j]=dp[i][j-1]+dp[i-1][j-1];
    }
    printf("%lld\n",dp[n][m]);
    return 0;
}

2161: 布娃娃

偏不写扫描线，考虑对于同一种p，询问编号越小越先得到答案，先离散化按p排序建立线段树，询问挂在叶节点上，维护区间最小询问。

按照c从大到小的顺序加入查询，区间-1，若最小值变成0，则说明这个区间是某个p的答案，二分删除这个询问。

易知每个询问最多被删除一次，所以时间复杂度仍为$O(nlog_2n)$

#include<map>
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
#define Mod 19921228
#define lll spc<<1
#define rrr spc<<1|1
typedef long long lnt;
struct data{
    lnt add;
    lnt first;
    lnt mod;
    lnt prod;
    void Insert(void)
    {
        scanf("%lld%lld%lld%lld",&add,&first,&mod,&prod);
        return ;
    }
    lnt sta(void)
    {
        return first%mod;
    }
    lnt S(lnt lst,lnt i)
    {
        return (prod*lst+add+i)%mod;
    }
}P,C,L,R;
struct pnt{
    int p;
    int l,r;
    int no;
    int c;
    void sta(void)
    {
        no=1;
        p=P.sta();
        l=L.sta();
        r=R.sta();
        c=C.sta();
        return ;
    }
    void ch(void)
    {
        if(l>r)std::swap(l,r);
        return ;
    }
}p[100001];
struct trnt{
    int minval;
    int to;
    int h;
    int lzt;
}tr[1300000];
int cnt;
int tot;
lnt ans;
int n;
int a[300001];
std::map<int,int>M;
std::vector<int>v[300001];
bool cmp(pnt x,pnt y)
{
    return x.c>y.c;
}
void pushup(int spc)
{
    tr[spc].minval=std::min(tr[lll].minval,tr[rrr].minval);
    return ;
}
void add(int spc,int v)
{
    tr[spc].minval+=v;
    tr[spc].lzt+=v;
    return ;
}
void pushdown(int spc)
{
    if(tr[spc].lzt)
    {
        add(lll,tr[spc].lzt);
        add(rrr,tr[spc].lzt);
        if(!tr[spc].to)tr[spc].lzt=0;
    }
    return ;
}
void build(int l,int r,int spc)
{
    if(l==r)
    {
        tr[spc].h=0;
        tr[spc].to=l;
        if(v[l].size())tr[spc].minval=v[l][0];
        else tr[spc].minval=0x3f3f3f3f;
        return ;
    }
    int mid=(l+r)>>1;
    build(l,mid,lll);
    build(mid+1,r,rrr);
    pushup(spc);
    return ;
}
void update(int l,int r,int ll,int rr,int spc)
{
    if(l>rr||ll>r)return ;
    if(ll<=l&&r<=rr)
    {
        add(spc,-1);
        return ;
    }
    pushdown(spc);
    int mid=(l+r)>>1;
    update(l,mid,ll,rr,lll);
    update(mid+1,r,ll,rr,rrr);
    pushup(spc);
    return ;
}
void check(int spc,int V)
{
    if(tr[spc].minval>0)return ;
    if(tr[spc].to)
    {
        ans=(ans+V)%Mod;
        tr[spc].h++;
        int i=tr[spc].to;
        if(tr[spc].h>=v[i].size())tr[spc].minval=0x3f3f3f3f;
        else{
            tr[spc].minval=v[i][tr[spc].h]+tr[spc].lzt;
        }
        return ;
    }
    pushdown(spc);
    check(lll,V);
    check(rrr,V);
    pushup(spc);
    return ;
}
int main()
{
    scanf("%d",&n);
    P.Insert(),C.Insert(),L.Insert(),R.Insert();
    p[1].sta();
    for(int i=2;i<=n;i++)
    {
        p[i].no=i;
        p[i].p=P.S(p[i-1].p,i);
        p[i].l=L.S(p[i-1].l,i);
        p[i].r=R.S(p[i-1].r,i);
        p[i].c=C.S(p[i-1].c,i);
        a[++cnt]=p[i].l;
        a[++cnt]=p[i].r;
        a[++cnt]=p[i].p;
    }
    a[++cnt]=p[1].l;
    a[++cnt]=p[1].r;
    a[++cnt]=p[1].p;
    std::sort(a+1,a+cnt+1);a[0]=-1;
    for(int i=1;i<=cnt;i++)
    {
        if(M.find(a[i])==M.end())M[a[i]]=++tot;
    }
    for(int i=1;i<=n;i++)
    {
        p[i].l=M[p[i].l];
        p[i].r=M[p[i].r];
        p[i].p=M[p[i].p];
        v[p[i].p].push_back(i);
    }
    for(int i=1;i<=n;i++)p[i].ch();
    build(1,tot,1);
    std::sort(p+1,p+n+1,cmp);
    for(int i=1;i<=n;i++)
    {
        update(1,tot,p[i].l,p[i].r,1);
        check(1,p[i].c);
    }
    printf("%lld\n",(ans%Mod+Mod)%Mod);
    return 0;
}