2018年4月10日训练日记

这两天补了省赛训练2的题目，以及做了手速练习四（全是原题）

省赛训练2：（完整） （zoj3936-zoj3947部分题目）

A：签到题，输出C，B+D，A，B+D

B：树分治，在网上好像只有浅蓝写了题解，然后看不懂，好像是cf上的Bear and Bowling 4的树上版。待以后补。

C：根据题意模拟即可

D：打表找规律，不难发现周期是400年（天数正好是7的倍数，为2058天），然后取模搞一下就行。我终于发现我比赛的时候为什么蜜汁WA了，我是记录在前，算这一天是不是幸运周在后，所以最后求的时候要加一。。。然后一加1一下就过了。。。

附代码：

#include <bits/stdc++.h>
#define ll long long
#define inf 0x3f3f3f3f
#define MAXN 200005
using namespace std;
const int maxn=50010;
const int mo=1e9+7;
int u,n,k,h,f;
int ans,tmp;
int mp[10000][13][32];
int y,m,d;
struct node
{
    int n,y,r;
}fmp[maxn];
bool jd(int year)
{
    int run=1;
    if (year%4 != 0) run = 0;
    else {
        if (year%100==0 && year%400!=0) run = 0;
        //if (year%3200 == 0) run = 0;
    }
    return run?1:0;
}
int day[]={0,31,28,31,30,31,30,31,31,30,31,30,31};
int main() {
    int T,t,cas;
    int dn=1753,dy=1,dr=0;
    int num=0;
    int sum=0;
    for(int i=0;i<=3009790;i++)
    {
        //if(dn>=1753+400) break;
        dr++;
        if(dy==2&&jd(dn))
        {
            if(dr==30) {dy++;dr=1;}
        }
        else if(dr>day[dy])
        {
            dy++;dr=1;
            if(dy>12) {dy=1;dn++;}
        }
        mp[dn][dy][dr]=num;//保存在前，所以最后需要+1
        if((dr==1||dr==11||dr==21)&&(i%7==0)){
            num++;fmp[num].n=dn;fmp[num].y=dy;fmp[num].r=dr;//计数
            }
    }
    while(cin>>T)
    {
        while(T--){
            int x,y,z;
            cin>>x>>y>>z>>n;
            n--;
            int tmp=n/2058;//400年正好是2058个幸运周
            n%=(2058);
            int mm=0;
            while(x>1753+400)
            {
                mm++;
                x-=400;
            }
            node ans=fmp[mp[x][y][z]+n+1];//比赛的时候这里少了+1
            cout<<(ll)((ll)(tmp+mm)*400+(ll)ans.n)<<" "<<(ll)ans.y<<" "<<(ll)ans.r<<endl;
    }
    }
    return 0;
}

E：优先队列，考虑一个数对一个区间取模，如果区间右端点大于等于这个数，就需要拆区间合并相同右端点区间（注意合并区间的细节，具体见代码），剩余区间入队，用优先队列维护最大右端点区间，记录合并后的区间数目，把结构体按右端点大小排序，然后保存每个区间的右端点和大于等于该右端点的区间总数，对于询问用lower_bound直接找出大于等于询问端点的下标，最后按题意相加即可。第一次交是TLE，然后加个注释再怎么交都是AC，迷，这道题仔细想想其实也不是很难

附代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=100010;
const int mo=1e9+7;
int ans[maxn];
int f,n,m,h;
int c[maxn],k,sum[maxn];
int b[maxn],flag,tmp;
int dp[maxn];
struct node
{
    int r,num;
    node(){r=0,num=0;}
    node(int x,int y){r=x,num=y;}
    bool operator<(node a)const{return r<a.r;}
}a[maxn];
int main()
{
    int T;
    int cas=1;
    while(cin>>T)
    {
        while(T--)
        {
            int n,m,s;
            cin>>n>>m;
            priority_queue<node> q;
            q.push(node(m,1));
            while(n--)
            {
                cin>>s;
                while(q.top().r>=s)
                {
                    node tmp=q.top();q.pop();
                    while(!q.empty()&&tmp.r==q.top().r)
                    {
                        tmp.num+=q.top().num;
                        q.pop();
                    }
                    q.push(node(s-1,(tmp.r+1)/s*tmp.num));
                    if((tmp.r+1)%s) q.push(node(tmp.r%s,tmp.num));
                }
            }
            int cnt=1;
            while(!q.empty())
            {
                node tmp=q.top(); q.pop();
                while(!q.empty()&&tmp.r==q.top().r)
                {
                    tmp.num+=q.top().num;
                    q.pop();
                }
                a[cnt++]=tmp;
            }
            sort(a+1,a+cnt);
            memset(sum,0,sizeof(sum));
            for(int i=1;i<cnt;i++)
            {
                ans[i]=a[i].r;
                sum[cnt-i]=sum[cnt-i+1]+a[cnt-i].num;
            }
            int qq;
            cin>>qq;
            ll anss=0;
            for(int i=1;i<=qq;i++)
            {
                cin>>s;
                int p=lower_bound(ans+1,ans+cnt,s)-ans;
                anss=(anss+(ll)i*(ll)(sum[p]%mo))%mo;
            }
            cout<<anss<<endl;
        }
    }
return 0;
}

F：DFS+贪心，最小的区间左端点是肯定要取的，然后考虑右端点是否取。我比赛的时候没什么思路，队友做的，漂亮的一血（然后就跑了），他的码风还真是简洁。。。连个const都没有

附代码：

#include<bits/stdc++.h>
using namespace std;
struct qj
{
    long long l,r;
};
int cmp(qj a,qj b)
{
    if(a.l==b.l)
        return a.r<b.r;
    return a.l<b.l;
}
long long m,k,sum,o;
qj a[50],b[50];
void dfs(long long val,long long i,long long ans,long long n)
{
    long long len,num,e,f,g;
	if(n>m)
		return;
    if((i>=o+1)||n==m)
    {
        sum=max(sum,ans);
		return;
    }
    if(val<b[i].r)
    {
        len=min(b[i].r-b[i].l+1,b[i].r-val);
        g=len%k;
        if(g==0)
            f=len/k;
        else
            f=len/k+1;
        if(f+n>m)
            sum=max(sum,ans+(m-n)*k);
        e=max(b[i].l+f*k-1,val+f*k);
        num=f*k;
        dfs(e,i+1,ans+num,n+f);
        dfs(b[i].r+k-1,i+1,ans+num+b[i].r+k-1-e,n+f+1);
    }
    else if(val>=b[i].r)
    {
        dfs(val,i+1,ans,n);
        dfs(b[i].r+k-1,i+1,ans+b[i].r+k-1-val,n+1);
    }
}
int main()
{
     long long t,x;
     int i,j;
     cin>>t;
     while(t--)
    {
         cin>>x>>k>>m;
         for(int i=1;i<=x;i++)
             cin>>a[i].l>>a[i].r;
         sort(a+1,a+x+1,cmp);
         o=0;
         for(i=1;i<=x;)
         {
             for(j=i+1;j<=x;j++)
             {
                if(a[j].l<=a[i].r+1)
                a[i].r=max(a[i].r,a[j].r);
                else break;
             }
             o++;
			 b[o]=a[i];
             i=j;
        }
        sum=0;
        dfs(0,1,0,0);
        cout<<sum<<endl;
     }
}

G：贪心。。。直接按照题目的形状判断即可。。。我当时写了8000+B，直接判断每个格子的每一种情况，也是够暴力的，人家几行就出来了，迷。。。（这就是错失一血的理由？代码又长又笨就不贴了）

H：多简单的一道spfa题目，求最短路的时候顺便记录最少价格就好了，比赛的时候居然没思路。。。傻

附代码：

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=200010;
int ans,f,n,m,h,k;
int a[maxn],c[maxn],t;
struct node
{
    int v;
    ll d,c;
    int nex;
}e[maxn*2];
int he[maxn],vis[maxn];
ll dis[maxn],di[maxn],sum,cas;
void add(int u,int v,int d,int c)
{
    e[t].v=v;
    e[t].d=d;
    e[t].c=c;
    e[t].nex=he[u];
    he[u]=t++;
}
void spfa(int st)
{
    int u,v;
    queue<int> q;
    q.push(st);
    dis[st]=0;
    vis[st]=1;
    while(!q.empty())
    {
        u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=he[u];i!=-1;i=e[i].nex)
        {
            v=e[i].v;
            if(dis[u]+e[i].d==dis[v])
            {
                if(di[v]>e[i].c)
                {
                    dis[v]=dis[u]+e[i].d;
                    di[v]=e[i].c;//保留最小价格
                }
                if(!vis[v]) vis[v]=1,q.push(v);
            }
            if(dis[u]+e[i].d<dis[v])
            {
                dis[v]=dis[u]+e[i].d;
                di[v]=e[i].c;//最短路是前提
                if(!vis[v]) vis[v]=1,q.push(v);
            }
        }
    }
}
void init()
{
    t=0;
    memset(vis,0,sizeof(vis));
    memset(dis,0x3f3f3f3f,sizeof(dis));
    memset(di,0x3f3f3f3f,sizeof(di));
    memset(he,-1,sizeof(he));
}
int main()
{
    int T;
    int cas=1;
    scanf("%d",&T);
    while(T--)
    {
       init();
       scanf("%d %d",&n,&m);
        for(int i=0;i<m;i++)
        {
           int x,y;
           ll d,c;
           scanf("%d%d%lld%lld",&x,&y,&d,&c);
           add(x,y,d,c);
           add(y,x,d,c);
        }
        spfa(0);
        ll ans=0,anss=0;
        for(int i=1;i<n;i++)
        ans+=dis[i];
        for(int i=1;i<n;i++)
        anss+=di[i];
        printf("%lld %lld\n",ans,anss);
    }
return 0;
}

I：n加上输入的每一个数的总和

对于昨晚的手速练习赛，都是原题，还真是比手速哇。。。有两道题一两秒的差距错失一血，论手速的重要性。。。C题忘了题意了，WA了三次，迷。。。后面有几道当时没做出来的题，幸好赛后补了，一个小时成功AK。。。

大部分是贪心暴力模拟就不说了，值得提一下的是F题，在降价允许的范围内使结尾的9最多，直接+1判断0最多，枚举一下余数就好了，十分简单，最后别忘了-1。还有J题，只需要交换行标和列标就可以了。。。

博客资料终于来了，也是时候总结和准备模板了（之前准备了一些常用的小技巧和一些知识点的模板）。

还有三个周省赛，不要让这一年多的时间白费，拼出自己的实力！！！