The 2017 ACM - ICPC Asia Ho Chi Minh City Regional Contest

---

大一暑假第一场组队赛
2019/7/23 NCC_2590 Valkyrie dyx,jxt,cj

---

题 B - Barcode Kattis - barcode

题目：为了准备2017年的acm赛，主办方HCMUE决定在参赛者的T恤上打一个标签。这个标签要求是越简单越好，来减少成本，但是又体现科技风格。
HCMUE决定让每个标签都有红蓝条，至少满足下面的条件之一:
1.蓝色条纹数与红色条纹数相等
2.蓝色条纹不相邻
现在有k种方式来创建要求的这种包含N个条纹的标签。给出2个数N和M,M是一个指数，你的任务是帮他找到K对M取模后的余数。

题目涉及：排列组合数学知识，斐波那契规律，用lucus求解二项式，费马引理，快速幂，细节取模处理

解题思路：
一.确定求解模型
求蓝色条纹不相连的情况有几种的方法：
如有7个条纹，从3R4B开始，3r创造4个空，则c4里选4.
3R4B开始，则 $C_4^4$ ,
4R3B开始，则 $C_5^3$ ,
5R2B，则 $C_6^2$ ,
6R1B，则 $C_7^1$ ,
7R0B，则 $C_8^0$ ,
这些加起来之后就是7个条纹情况下，蓝色条纹不相邻的排列组合个数。
而打表后发现，这符合斐波那契数列规律。

当n是个奇数时，不存在红蓝个数相等，因此必须满足蓝色条纹不相邻，答案即为fl（n）；

当n是个偶数时，存在满足条件1的，满足条件2的，2个条件都满足的，求并集。
满足条件1的，即当红蓝个数相等时（在n个位置里，挑n/2的位置，那么剩下的就是另一个颜色的位置了），即$C_n^{n/2}$，得到数A；
满足条件2的，得到数B；
2个条件都满足的，即满足红白个数相等又满足蓝色条纹不相连的情况，即$C_{n/2+1}^{n/2}$,得到C。
故n为偶数时的答案为A+B-C=fl（n）+$C_n^{n/2}$-$C_{n/2+1}^{n/2}$

二.寻找算法。
而现在的$C_b^a$很大，
如果用 for循环，利用杨辉三角形计算，for(int i=1;i<=a;i++)c[i]=c[i-1]*(b-i+1)/i;，循环太多，并且循环一次取一次模是错的，如下

for(int i=1;i<=a;i++)c[i]=c[i-1]*(b-i+1)/i%mod;   //因为每次（b-i+1）/i 不一定是一个整数，故这样得出的答案不是原来的那个二项式取模

如果用 $C_b^{a}mod(g)=\frac{b!}{a!\ast (a-b)!}mod(g)$,那么a和b又太大，
所以就用了 lucas定理 https://www.cnblogs.com/Liuz8848/p/11019347.html

$C_b^{a}mod(g)=C_{b\%g}^{a\%g}\ast C_{b/p}^{a/p}\mathrm{m}\mathrm{o}\mathrm{d}(g)$,

$C_{bmod(g)}^{amod(g)}=C_p^q=\frac{p!}{q!\ast (p-q)!}\mathrm{m}\mathrm{o}\mathrm{d}(g)$,

此时不可以写成 $\frac{p!}{q!\ast (p-q)!}\mathrm{m}\mathrm{o}\mathrm{d}(p)=\frac{p!\mathrm{m}\mathrm{o}\mathrm{d}(p)}{q!\ast (p-q)!\mathrm{m}\mathrm{o}\mathrm{d}(p)}$,故要用逆元 ${(q!\ast (p-q)!)}^{p-2}$

———————————————————————————————————————————————
逆元
对于(ab)%mod,那么可以((a%mod)(b%mod))%mod，但是若是（a/b）%mod，则不可以如此，
则问有没有办法转化为（a*c）%mod,使之可以用那个方法，故c称为b的逆元，c=$a^{mod-2}$
证明：
$\because$由费马定理得 如果p是素数，则$a^{p-1}\equiv$ 1(mod p).
$\therefore a\ast a^{p-2}\equiv$ 1(mod p).
$\therefore a^{p-2}是a的逆元$
———————————————————————————————————————————————

解题坑点：
1、要用long long int，用int 只能过到第4个样例
2、对a*=a%mod 与a=a*a%mod 结果不一样
3、 a%mod 要用 （a%mod+mod）%mod 保证一下正数。
4.计算$C_b^a$时，要注意保证a<=b。

解题代码
如下，但还有改进的空间，快速幂可以用非递归形式，lucas也可以用非递归。

#include<cstdio>
#include<algorithm>
#define ll long long int
using namespace std;
const int maxn=1e6+10;
ll fibans[maxn];
ll fel(ll a,ll b)
{
    fibans[1]=2;
    fibans[2]=3;
    for(int i=3;i<=a;i++)fibans[i]=(fibans[i-1]+fibans[i-2])%b;
    return fibans[a]%b;
}
ll quick_pow(ll a,ll b,ll mod)
{
  if(b==0)return 1;
  ll get=quick_pow(a,b/2,mod);
  get=get*get%mod;
  if(b&1)return get*a%mod;
  else return get%mod;
}
ll com(ll up,ll down,ll mod)
{
    if(up>down)return 0;
    if(up==down)return 1;
    up=min(down-up,up);
    ll coa=1,cob=1;
    for(int i=1;i<=up;i++)
    {
        coa=coa*(down-i+1)%mod;
        cob=cob*i%mod;
    }
    ll niyuan=quick_pow(cob,mod-2,mod);
    return ((coa*niyuan)%mod+mod)%mod;
}
ll lucas(ll up,ll down,ll mod)
{
   if(!up||!down)return 1;
  return com(up%mod,down%mod,mod)*lucas(up/mod,down/mod,mod)%mod;
}
int main(void)
{
    int cas;
    scanf("%d",&cas);
    ll n,m,ans;
    for(int ca=1;ca<=cas;ca++)
    {
        scanf("%lld %lld",&n,&m);
        if(n&1)ans=fel(n,m);
        else ans=((lucas(n/2,n,m)+fel(n,m)-lucas(n/2,n/2+1,m))%m+m)%m;
        printf("%lld\n",ans);
    }
    return 0;
}

题H - Height Preservation Kattis - heightpreservation

题目：VR是最近的一个新应用趋势。穿戴VR头盔，你就会有一种虚拟现实的环境来教育、合作、娱乐等等。为了增强用户体验，ICPC建立了一间VR房间让用户走动和探索虚拟环境。用户可以通过头戴式设备来看到虚拟环境，感受到虚拟地形的坡度当在这间VR房间行走时。
这个VR房间的地板由mn个单元格组成，m行n列。每个单元格的真实高度可以被调整所以我们模拟了一个虚拟地势。
第i行和第j列的单元格分界线上有一个真实高度$S_{i,j}$,回应了虚拟画面中的虚拟高度$H_{i,j}$。这个现实高度对虚拟高度的维护，只需维护了同一行与同一列中的单元格的高度关系。
如第i行，如果$S_{i,j}$=$S_{i,k}$,那么$H_{i,j}$=$H_{i,k}$, 大小比较关系一致。在每一列中，也是这个情况。
现在给出每个单元格的虚拟高度，你的任务是决定最少种不同种类的真实高度。
输入：m和n，mn个单元格虚拟高度$H_{i,j}$
输出：最少种类的单元格
思路：
对点按大小排完序了之后，就要填补点了，如果有大小差距肯定要新增加一个点，那么种类数就得加一，每一次就是更新答案，
row_num和list_num表示更新到那一行，那一列的时候的答案。

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
struct re
{
    int x,y,height;
    bool operator <(const re &a)const
    {
        if(height!=a.height)return height<a.height;
    if(x!=a.x)return x<a.x;return y<a.y; }
};
const int maxn=1e6+10;
re all[maxn];
int row_num[maxn],list_num[maxn];
int row_max[maxn],list_max[maxn];
int main(void)
{
    int r,d,cnt=0;
    scanf("%d %d",&r,&d);
    for(int i=1;i<=r;i++)
        for(int k=1;k<=d;k++)
    {
        all[++cnt].x=i;
        all[cnt].y=k;
        scanf("%d",&all[cnt].height);
    }
    int ans=0;
    memset(row_num,0,sizeof(row_num));
    memset(list_num,0,sizeof(list_num));
    memset(row_max,0,sizeof(row_max));
    memset(list_max,0,sizeof(list_max));
    sort(all+1,all+1+cnt);
    for(int i=1;i<=cnt;i++)
    {
        int nx=all[i].x,ny=all[i].y,hei=all[i].height;
        if(hei>row_max[nx])row_num[nx]++;
        if(hei>list_max[ny])list_num[ny]++;
        int temans=max(row_num[nx],list_num[ny]);
        row_num[nx]=list_num[ny]=temans;
        ans=max(ans,temans);
        row_max[nx]=list_max[ny]=hei;//因为是已经排序的，所以当前高度肯定是最大值。
    }
    printf("%d\n",ans);
    return 0;
}

E - Engaging with Loyal Customers

Kattis - engaging
类型：KM算法的模板
题解：
1.首先要知道这是个KM算法的模型，二分匹配，求最优匹配。
https://www.cnblogs.com/wenruo/p/5264235.html
2.其次要知道如何对KM算法提速，因为女生越少，那么递归的次数越少，所以要让个数少的那个做选择方，做女生。
如下，这是在开始的时候做的处理，但是在输出答案的时候也要注意处理。

flag=0;
    scanf("%d %d %d",&m,&n,&k);
    if(n>m){swap(m,n);flag=1;}

    memset(pipei,0,sizeof(pipei));
    for(int i=1;i<=k;i++)
    {
        int a,b,c;
        scanf("%d %d %d",&a,&b,&c);
        if(flag)swap(a,b);
        pipei[a][b]=c;
    }