4.28 KMP（重点next数组）+线段树+map二分

KMP算法水题
String Freshman
KMP算法考察，会理解过KMP算法就是水题。

KMP算法的主要有两个方面
1.第一个方面是对于子串和主串匹配不成功的情况下子串如何根据next数组变化
2.如何获得子串的next数组

期中较为难理解的就是第二点。
如何理解next数组的由来，

我们来用一个例子来看一下
ABABAA

这样就能看出next数组的机制，也就是对于原本的数组的各个前缀来讲，如果对于这个前缀的长度一致前后缀相等，那么取最长的前后缀相等的长度作为next数组的值。
附KMP模板

void pre()
{
    p[1]=0;int j=0;
    for(int i=1;i<m;i++)
    {
        while(j>0&&s2[j+1]!=s2[i+1])
            j=p[j];
        if(s2[j+1]==s2[i+1])
            j++;
        p[i+1]=j;
    }
}
void kmp()
{
    int ans=0,j=0;
    for(int i=0;i<n;i++)
    {
        while(j>0&&s2[j+1]!=s1[i+1])
            j=p[j];
        if(s2[j+1]==s1[i+1]) j++;
        if(j==m)
        {
		    //ans++;匹配数量确认
            printf("%d\n",i-m);
            j=p[j];
            //j=0不可重叠确认
        }
    }
}

Bin Packing Problem
线段树的返回下标与map的自带二分
维护区间最大值和返回第一个合适的目标值的下标

#include <bits/stdc++.h>
using namespace std;
#define N 1000000+100
#define ls p<<1
#define rs p<<1|1
#define mid ((l+r)>>1)
int a[N];
int n,c;
int tree[N<<2];
void up(int p)
{
    tree[p]=tree[ls]>tree[rs]?tree[ls]:tree[rs];
}
void build(int l,int r,int p)
{
    if(l==r)
    {
        tree[p]=c;
        return;
    }
    build(l,mid,ls);
    build(mid+1,r,rs);
    up(p);
}
void update(int l,int r,int pos,int p,int k)
{
    if(l==r)
    {
        tree[p]-=k;
        return;
    }
    if(pos<=mid)
    update(l,mid,pos,ls,k);
    if(pos>mid)
    update(mid+1,r,pos,rs,k);
    up(p);
}
int query(int l,int r,int p,int k)
{
    if(l==r)
    {
        return l;
    }
    int ans=0;
    if(tree[ls]>=k)
    {
        ans=query(l,mid,ls,k);
    }
    else if(tree[rs]>=k)
    {
        ans=query(mid+1,r,rs,k);
    }
    return ans;
}
int main()
{
    int t;
    for(cin>>t;t;t--)
    {
        cin>>n>>c;
        for(register int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
        }
        build(1,n,1);
        int size1=0;
        for(int i=1;i<=n;i++)
        {
            int bag=query(1,n,1,a[i]);
            if(bag>size1)
            {
                size1++;
            }
            update(1,n,bag,1,a[i]);
        }
    return 0;
}

第二问用到了map和map自带的二分

#include <bits/stdc++.h>
using namespace std;
#define N 1000000+100
int a[N];
int n,c;
int main()
{
    int t;
    for(cin>>t;t;t--)
    {
        cin>>n>>c;
        for(register int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
        }
        map<int,int> mp;
        mp[c-a[1]]++;
        for(int i=2;i<=n;i++)
        {
            map<int,int>::iterator p=mp.lower_bound(a[i]);
            if(p==mp.end())
            {
                mp[c-a[i]]++;
            }
            else
            {
                mp[p->first]--;
                mp[(p->first)-a[i]]++;
                if(mp[p->first]==0)
                    mp.erase(p);
            }
        }
        int size2=0;
        for(auto x:mp)
        {
            size2+=x.second;
        }
        printf("%d %d\n",size1,size2);
    }
    return 0;
}

来了一个傻傻的精度题
A - New Year Table
注意精度控制一般在1e-8到1e-12之间

#include <bits/stdc++.h>
using namespace std;
double pi=acos(-1);
double eps=1e-10;
int main()
{
    double n,R,r;
    cin>>n>>R>>r;
    if(R<r)
    {
        cout<<"NO"<<endl;
        return 0;
    }
    if(r>R/2.0)
    {
        if(n==1)
            cout<<"YES"<<endl;
        else
            cout<<"NO"<<endl;
        return 0;
    }
    else
    {
        double t=asin(r/(R-r));
        if(n*2.0*t<=2*pi+eps)
            cout<<"YES"<<endl;
        else
            cout<<"NO"<<endl;
    }
    return 0;
}