[NOIP2011]刷水

　　前几天做了NOIP2011的题，感觉不是那么难。

　　这边先做了两天的前两题，T3还没打。

D1T1：顺次读入，分别判断是否覆盖即可，照例大水：

#include<cstdio>
int x,y,a[],b[],c[],d[],n,ans=-;
int main()
{
    scanf("%d",&n);
    for(int i=;i<=n;i++)scanf("%d%d%d%d",a+i,b+i,c+i,d+i);
    scanf("%d%d",&x,&y);
    for(int i=;i<=n;i++)
        if(a[i]<=x&&b[i]<=y&&c[i]+a[i]>=x&&d[i]+b[i]>=y)ans=i;
    printf("%d",ans);
    return ;
}

D1T2：一边扫过去，对于每家客栈，我们记下前一家相同色调的客栈，前一家能去的咖啡店，之前有多少家可以和这间客栈搭配和这是第几间该色调的客栈，然后如果前一家色调相同的客栈编号在前一家能去的咖啡店之前，说明这之前的客栈都能去，把可搭配客栈数更新后加入答案，否则直接加入答案。

代码如下：

#include<cstdio>
using namespace std;
int n,p,t[],h[],ans,l[],lst[];
int main()
{
    scanf("%d%*d%d",&n,&p);
    for(int i=;i<=;i++)lst[i]=;
    for(int i=;i<=n;i++)
    {
        int w,c;
        scanf("%d%d",&c,&w);
        t[i]=t[i-]+(w<=p);
        if(t[i]!=t[lst[c]-])h[c]=l[c];
        l[c]++;
        ans+=h[c];
        if(w<=p)h[c]=l[c];
        lst[c]=i;
    }
    printf("%d",ans);
    return ;
}

D2T1：用杨辉三角计算组合数然后用二项式定理直接出解：

#include<cstdio>
const int mod=;
int a,b,k,n,m,C[][],aa=,bb=;
void calC()
{
    C[][]=;
    for(int i=;i<=;i++)
        for(int j=;j<=i;j++)
            C[i][j]=(C[i-][j-]+C[i-][j])%mod;
}
int main()
{
    calC();
    scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);
    a%=mod;b%=mod;
    for(int i=;i<=n;i++)aa=(aa*a)%mod;
    for(int i=;i<=m;i++)bb=(bb*b)%mod;
    printf("%d",((aa*bb%mod)*C[k+][n+])%mod);
    return ;
}

D2T2：使和标准最近，而且是单调递减，容易想到用二分答案+check

二分就不说了。check中我们要O(m)计算m个区间的sigma，可以对于二分出的x计算从1到n能加入答案的计算前缀和然后轻松处理。

要注意防止爆long long，这边处理方法是加个最优上限maxn，否则一定不优。

代码如下：

#include<cstdio>
const long long maxn=10000000000000ll;
long long mn(long long x,long long y){return x<y?x:y;}
long long s[],ans=maxn,sa,v[],w[],sum[],l[],r[],ll,rr,n,m,mid;
bool check(long long x)
{
    long long tt=;
    for(int i=;i<=n;++i)
        sum[i]=sum[i-]+(w[i]>=x),
        s[i]=s[i-]+(w[i]>=x?v[i]:);
    for(int i=;i<=m;++i)
    {
        tt+=(sum[r[i]]-sum[l[i]-])*(s[r[i]]-s[l[i]-]);
        if(tt>maxn)return ;
    }
    ans=mn(-mn(tt-sa,sa-tt),ans);
    return tt<=sa;
}
void bina()
{
    ll=;rr=;
    while(ll<rr)
    {
        mid=(ll+rr)/;
        check(mid)?rr=mid:ll=mid+;
    }
}
int main()
{
    scanf("%lld%lld%lld",&n,&m,&sa);
    for(int i=;i<=n;++i)scanf("%lld%lld",w+i,v+i);
    for(int i=;i<=m;++i)scanf("%lld%lld",l+i,r+i);
    bina();
    check(mid-);check(mid);
    printf("%lld",ans);
    return ;
}