2018 hdu 多校 4

这场的题感觉可作性高了……、
补题进度[5/12] 场上3题

Problem D. Nothing is Impossible
这道题大家众说纷纭，争议比较大吧。代码不是我写的，我个人的理解是鸽巢原理，既b+1个人里一定有1个人能作对，然后按b排序，一点一点取就是了

#include<bits/stdc++.h>
using namespace std;
int co[105];
int main()
{
    int T;
    cin>>T;
    int i,j,k;
    int p,n,a,b;
    while(T--)
    {
        cin>>p>>n;
        for(i=0;i<p;i++)
        {
            scanf("%d %d",&a,&b);
            co[i]=b;
        }
        sort(co,co+p);
        i=0;

        int cnt=0;
        while(n && i<p)
        {
            if(n>=co[i]+1)
            {
                cnt++;
            }
            n/=(co[i]+1);
            i++;
        }
        cout<<cnt<<endl;
    }
}

Problem K. Expression in Memories
无脑模拟就是

Problem L. Graph Theory Homework
题意是：一个完全图，两个点之间的边权是点权绝对值的开方向下取整，求1到n的最短路长度
这道题的正解是均值不等式，但是我猜了一个最短路经过的点最多三个，然后就过了，其实只要输出起点终点那条变得路径就行了。

Problem B. Harvest of Apples
为什么大家都能想到是莫队……自闭了
有1e5的查询无修改，就联想到mo算法，我们把$C_n^m$ 的m看作是l，n看作是r,然后我们利用组合数杨辉三角的性质推一下……大概就知道怎么转移了

#include <iostream>
#include <cstring>
#include <queue>
#include <cstdio>
#include <algorithm>
#include <set>
#include <map>
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int mod = 1000000007;
const int maxn = 1e5 + 5;
ll inv[maxn],fac[maxn];
void init()
{
    inv[0]=fac[0]=1;
    inv[1]=1;
    for(int i=1; i<maxn; i++)
        fac[i]=fac[i-1]*i%mod;
    inv[1]=1;
    for(int i=2; i<maxn; i++)
        inv[i]=(ll)(mod - mod/i)*inv[mod%i]%mod;
    inv[0]=1;
    for(int i=1; i<maxn; i++)
        inv[i]=inv[i-1]*inv[i]%mod;
}
ll com(int n,int m)
{
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int unit;
ll ans[maxn],val;
struct Query
{
    int n,k,id;
} q[maxn];
int T;
int cmp(Query a, Query b)
{
    return a.k/unit == b.k/unit ? a.n < b.n : a.k < b.k;
}
void solve()
{
    int L = 0,R = 1;
    val=1;
    for(int i = 1; i<=T; i++)
    {
        // printf("%d %d\n",q[i].k,q[i].n);
        while(R < q[i].n)
            val = (2*val- com(R++,L) + mod) % mod;
        while(R > q[i].n)
            val = ((val + com(--R,L))*inv[2]) % mod;
        while(L < q[i].k)
            val = (val + com(R,++L)) % mod;
        while(L > q[i].k)
            val = (val + mod - com(R,L--)) % mod;
        //cout<<L<<' '<<R<<endl;
        ans[q[i].id] = val;
        //cout<<L<<' '<<R<<endl;
    }
}
int main()
{
    // freopen("out.txt","w",stdout);
    scanf("%d",&T);
    init();
    unit = sqrt(10000);
    for(int i = 1; i<=T; i++)
    {
        scanf("%d %d",&q[i].n,&q[i].k);
        q[i].id = i;
    }
    sort(q+1,q+T+1,cmp);
    solve();
    for(int i = 1; i<=T; i++)
    {
        printf("%I64d\n",ans[i]);
    }
    return 0;
}

Problem J. Let Sudoku Rotate
不怎么会写旋转，就学习了dls的旋转操作，其他的操作还是比较好写的

#include <iostream>
#include <cstring>
#include <queue>
#include <cstdio>
#include <algorithm>
#include <set>
#include <map>
using namespace std;
int s[16][16],w[16][16];
char tmp[16][16];
void rot(int a,int b)
{
    for(int i = 0; i<4; i++)
    {
        for(int j = 0 ; j <4; j++)
        {
            w[j][3 - i] = s[a+i][b+j];
            //printf("%d ",w[i][j]);
        }
       // printf("\n");
    }
    for(int i = 0;i< 4;i++)
    {
        for(int j = 0 ;j <  4;j++)
        {
            s[a+i][b+j] = w[i][j];
        }
    }
}
int vis[16];
int ans,cnt;
bool check(int a,int b)
{
    for(int i = a;i < a+ 4;i++)
    {
        cnt++;
        for(int j = 0;j<b+4;j++)
        {
            if(vis[s[i][j]] == cnt) return 0;
            vis[s[i][j]] = cnt;
        }
    }
    for(int i = b;i<b+4;i++)
    {
        cnt++;
        for(int j = 0;j<a+4;j++)
        {
            if(vis[s[j][i]] == cnt) return 0;
            vis[s[j][i]] = cnt;
        }
    }

    return 1;
}
void dfs(int i,int j,int k)
{
    if(i == 4)
    {
        ans = min(ans,k);
        return ;
    }
    if(k >= ans) return;
    if(j == 4) dfs(i+1,0,k);
    for(int t = 0;t<4;t++)  // number of rotation
    {
        if(check(i*4,j*4))
        {
            dfs(i,j+1,k+t);
        }
        rot(i*4,j*4);
    }
}
int main()
{
    int ca;
    scanf("%d", &ca);
    while(ca--)
    {
        for(int i = 0;i<16;i++)
        {
            scanf(" %s",tmp[i]);
            for(int j = 0;j<16;j++)
            {
                if(tmp[i][j] >= '0' && tmp[i][j] <= '9') s[i][j] = tmp[i][j] - '0';
                else s[i][j] = tmp[i][j] - 'A' + 10;
            }
        }
        ans = 1000;
        dfs(0,0,0);
        printf("%d\n",ans);
    }
    return 0;
}
/*
2
5 2
1000 500
*/

Problem E. Matrix from Arrays
打了个表其实就会发现是由循环，按照循环的大小先打出来足够一部分，然后用二位前缀和解决这个问题。

#include <iostream>
#include <cstring>
#include <queue>
#include <cstdio>
#include <algorithm>
#include <set>
#include <map>
using namespace std;
int s[16][16],w[16][16];
char tmp[16][16];
void rot(int a,int b)
{
    for(int i = 0; i<4; i++)
    {
        for(int j = 0 ; j <4; j++)
        {
            w[j][3 - i] = s[a+i][b+j];
            //printf("%d ",w[i][j]);
        }
       // printf("\n");
    }
    for(int i = 0;i< 4;i++)
    {
        for(int j = 0 ;j <  4;j++)
        {
            s[a+i][b+j] = w[i][j];
        }
    }
}
int vis[16];
int ans,cnt;
bool check(int a,int b)
{
    for(int i = a;i < a+ 4;i++)
    {
        cnt++;
        for(int j = 0;j<b+4;j++)
        {
            if(vis[s[i][j]] == cnt) return 0;
            vis[s[i][j]] = cnt;
        }
    }
    for(int i = b;i<b+4;i++)
    {
        cnt++;
        for(int j = 0;j<a+4;j++)
        {
            if(vis[s[j][i]] == cnt) return 0;
            vis[s[j][i]] = cnt;
        }
    }

    return 1;
}
void dfs(int i,int j,int k)
{
    if(i == 4)
    {
        ans = min(ans,k);
        return ;
    }
    if(k >= ans) return;
    if(j == 4) dfs(i+1,0,k);
    for(int t = 0;t<4;t++)  // number of rotation
    {
        if(check(i*4,j*4))
        {
            dfs(i,j+1,k+t);
        }
        rot(i*4,j*4);
    }
}
int main()
{
    int ca;
    scanf("%d", &ca);
    while(ca--)
    {
        for(int i = 0;i<16;i++)
        {
            scanf(" %s",tmp[i]);
            for(int j = 0;j<16;j++)
            {
                if(tmp[i][j] >= '0' && tmp[i][j] <= '9') s[i][j] = tmp[i][j] - '0';
                else s[i][j] = tmp[i][j] - 'A' + 10;
            }
        }
        ans = 1000;
        dfs(0,0,0);
        printf("%d\n",ans);
    }
    return 0;
}