暑假集训日记——8.4( 点分治+康拓展开+组合数)

最新推荐文章于 2021-09-23 20:21:37 发布
最新推荐文章于 2021-09-23 20:21:37 发布
阅读量98 收藏
点赞数
分类专栏： 点分治 康拓展开 树状数组 二分查找法 思维题型/数论相关 容斥原理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43870114/article/details/98463648
版权

点分治详解：大佬的博客
P3806 【模板】点分治1
模板代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define inf 999999999
 
int n,m,len=0,Size;
struct node{int x,y,z,next;};
node e[20010];
int first[10010];
int root,ms,size[10010],mson[10010],sum[10000010];
bool v[10010];

void buildroad(int x,int y,int z)
{
    len++;
    e[len].x=x; e[len].y=y; e[len].z=z;
    e[len].next=first[x];
    first[x]=len;
}

//初始化以x为节点的子树
void getroot(int x,int fa)//fa表示x的父亲，防止往回搜
{
    size[x]=1;mson[x]=0;
    //size[i]记录以i为根的子树的大小，mson[i]表示i节点的最大子树的大小 
    for(int i=first[x];i;i=e[i].next)
    {
        int y=e[i].y;
        if(v[y]||y==fa)continue;
        //v[y]表示当前节点是否被分治过
        getroot(y,x); //往下继续搜索 
        size[x]+=size[y]; //加上子树大小 
        if(size[y]>mson[x]) mson[x]=size[y]; //更新最大的子树 
    }
    if(Size-size[x]>mson[x]) mson[x]=Size-size[x];
    //Size表示当前这整棵树的大小，那么Size-size[x]就表示不在x的子树内的节点数量，下面详解 
    if(ms>mson[x]) ms=mson[x],root=x;
    //ms表示树的重心的最大子树的大小（相当于mson[root]），这一步是用来更新树的重心的，root用来记录重心 
}

int t; //临时变量，记录dis的编号
int dis[10010];
int ask[110],ans[110];

//初始化每个节点到根节点x的距离
void getdis(int x,int fa,int z)//fa表示x的父亲，z表示x到目标点的距离
{
    dis[++t]=z;
    //这里写dis[++t]=z而不是dis[x]=z是因为只需要求得匹配数，所以我们只需要知道每一个点到目标点的距离，而不需要知道那个点是谁，
    for(int i=first[x];i;i=e[i].next)
    {
        int y=e[i].y;
        if(y==fa||v[y])continue;
        //这里加个v[y]是因为不能往被分治过的点走，分治过的点就没了
        getdis(y,x,z+e[i].z);
    }
}
struct asd{int x,y;}arr[10010];
int tt;
bool cmp(int x,int y){return x<y;}

//x表示getdis的起始点，y表示x到目标点的距离，id表示这一次统计出来的答案是合理的还是不合理的
void solve(int x,int y,int id)
{
    t=0;
    getdis(x,0,y);//算出树中的点到目标点的距离 
    tt=0;
    sort(dis+1,dis+t+1,cmp);
    dis[0]=-233;
    //合并相同距离的点，并记录数目
    for(int i=1;i<=t;i++)
    if(dis[i]!=dis[i-1])arr[++tt].x=dis[i],arr[tt].y=1;
    else arr[tt].y++;
    for(int i=1;i<=m;i++)
    {
    	//特判距离为k/2的点
    	if(ask[i]%2==0)for(int j=1;j<=tt;j++)
    		if(arr[j].x==ask[i]/2)ans[i]+=(arr[j].y-1)*arr[j].y/2*id;
    	//二分寻找其他满足距离和为k的点对，暴力求解复杂度为O(n^2),所以优化一下，枚举小于k/2的点，然后二分大于k/2的点，复杂度O(nlogn)
    	for(int j=1;j<=tt&&arr[j].x<ask[i]/2;j++)
    	{
    		int l=j+1,r=tt;
    		while(l<=r)
    		{
    			int mid=l+r>>1;
    			if(arr[j].x+arr[mid].x==ask[i])
    			{
    				ans[i]+=arr[j].y*arr[mid].y*id;
    				break;
                }
                if(arr[j].x+arr[mid].x>ask[i])r=mid-1;
                else l=mid+1;
            }
        }
    }
}

//分治操作 以x为重心的子树
void fenzhi(int x,int ssize)//ssize是当前这棵子树的大小
{
    v[x]=true;
    //代码保证每次进来的x都必定是当前这棵树的重心，我们将v[x]标记为true，表示x点被分治过
    solve(x,0,1);//计算这棵树以x为重心的所有组合，其中包含着不成立的情况
    for(int i=first[x];i;i=e[i].next)
    {
        int y=e[i].y;
        if(v[y])continue;//分治过，不能继续运算
        solve(y,e[i].z,-1);//计算不合法组合，用所有组合减去不合法组合
        ms=inf;root=0;//记得要初始化
        //分治的子树的大小
        Size=size[y]<size[x]?size[y]:(ssize-size[x]);
        getroot(y,0);//求出以y为根的子树
        fenzhi(root,(size[y]<size[x]?size[y]:(ssize-size[x])));
    }
}
 
int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<n;i++)
    {
        int x,y,z;
        scanf("%d %d %d",&x,&y,&z);
        buildroad(x,y,z);
        buildroad(y,x,z);
    }
    for(int i=1;i<=m;i++)
    scanf("%d",&ask[i]);
    root=0;ms=inf;Size=n;
    getroot(1,0);
    fenzhi(root,n);
    for(int i=1;i<=m;i++)
    if(ans[i]>0)printf("AYE\n");
    else printf("NAY\n");
}

A - Misha and Permutations Summation
题意：
给出两个排列，求出每个排列在全排列的排行，相加，模上n!（全排列个数)得出一个数k，求出排行为k的排列。

题解：康拓展开+树状数组+二分
按题意模拟即可

#include<bits/stdc++.h>
#define mp make_pair
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef long double ld;

const int N=1e6+10;
const int MAXN=20010;
const int INF=0x3f3f3f3f;
const double eps=0.0000001;
const ll mod=998244353 ;
int n,m,x,y,k,cnt;

int a[N];///维护康拓系数
int c[N];///用前缀和计算小于x的数的个数

int lowbit(int x){
    return x&-x;
}
int add(int x,int a)
{
    while(x<=n) c[x]+=a,x+=lowbit(x);
}

int sum(int x)
{
    int ans=0;
    while (x) ans+=c[x],x-=lowbit(x);
    return ans;
    //log2n;
}

int finds(int x)
{
    int l=1,r=n;
    while(l < r)///寻找x的上界
    {
        int mid=(l+r)>>1;
        int ans=sum(mid);
        if(ans< x) l=mid+1;
        else r=mid;
    }
    return l;
}

int main()
{
    scanf("%d",&n);
    memset(c,0,sizeof(c));
    ///初始化区间
    for(int i=1;i<=n;i++) add(i,1);
    int x;
    for(int i=1;i <= n; i++)
    {
        scanf("%d", &x);
        x++;///因为有数字0，所以整体后移一位
        ///找出剩余比x小的数的个数
        int temp=sum(x) - 1;
        ///删去已经询问过的数字
        add(x,-1);
        a[i]+=temp;
    }
    for(int i=1;i<=n;i++) add(i,1);
    for(int i=1;i <= n; i++)
    {
        scanf("%d", &x); x++;
        int temp=sum(x) - 1;
        add(x,-1);
        a[i]+=temp;
    }
    for(int i=n;i>=1;i--)
    {
        ///逆康拓展开，求余数
        a[i-1]+=a[i]/(n - i + 1);
        ///因为上式a[i]不可能大于(n-i+1)
        a[i] %=(n - i + 1);///对n!取余
    }
    ///逆康拓展开
    for(int i=1;i<=n;i++) add(i,1);
    for(int i=1;i<=n;i++)
    {
        int tmp=finds(a[i]+1);
        printf("%d ",tmp -1);
        add(tmp,-1);
    }
    printf("\n");

}

C - Color
组合数计算公式：
在这里插入图片描述
题解：（容斥+组合数+逆元）
容斥有 i个颜色没有选，i个颜色没有选的概率是Ai
在这里插入图片描述

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#define mp make_pair
using namespace std;

typedef long long LL;
typedef pair<int, int> pii;

typedef long double ld;


const int N=1e6+10;
const int MAXN=20010;
const int INF=0x3f3f3f3f;
const double eps=0.0000001;
const LL mod=1e9+7;

LL mul(LL a,LL b){return (a%mod*b%mod)%mod;}
LL add(LL a,LL b){return (a%mod+b%mod+mod)%mod;}
LL del(LL a,LL b){return (a%mod-b%mod+mod)%mod;}

LL inv[N];
LL pre[N];
int t,k,n,m;

LL qpow(LL a,LL b,LL M){
    LL r=1;
    while(b){
        if(b&1) r=r*a%M;
        b>>=1;
        a=a*a%M;
        }
    return r;

}


int cm[N],ck[N];
/// C(m,i) = C(m,i-1)*(m-i+1)/i;
void get_c()///计算组合数
{
    cm[0]=1,ck[0]=1;
    for(int i=1;i<=k;i++)
    {
        cm[i]=mul(mul(cm[i-1],m-i+1),inv[i]);
        ck[i]=mul(mul(ck[i-1],k-i+1),inv[i]);
    }
}

int main()
{
    inv[1] = 1;
    for(int i=2;i<=N;++i) inv[i]=qpow(i,mod-2,mod);///逆元打表
    scanf("%d",&t);
    LL s=t;
    while(t--)
    {
        scanf("%d%d%d",&n,&m,&k);
        get_c();
        LL fla=1,ans=0;
        for(int i=k;i>=1;i--)///容斥
        {
            ans=add(ans,mul(fla*ck[i]*i,qpow(i-1,n-1,mod)));
            fla=-fla;
            //cout<<ans<<" ";
        }

        ans=mul(ans,cm[k]);
        printf("Case #%d: ",s-t);
        printf("%lld\n",ans);
    }
    return 0;
}

D - Number Sequence
题解：（隔板法）+容斥
把 m个小球放到 n个盒子里,盒子可以为空的组合= C n + m − 1 n − 1 C_{n+m-1}^{n-1} Cn+m−1n−1​

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<map>
#define mp make_pair
using namespace std;

typedef long long LL;
typedef pair<int, int> pii;

typedef long double ld;


const int N=1e6+10;
const int MAXN=20010;
const int INF=0x3f3f3f3f;
const double eps=0.0000001;
const LL mod=1e9+7;

LL mul(LL a,LL b){return (a%mod*b%mod)%mod;}
LL add(LL a,LL b){return (a%mod+b%mod+mod)%mod;}
LL del(LL a,LL b){return (a%mod-b%mod+mod)%mod;}

LL C[505][505];
LL b[25] ;
map<int,int>m ;

int n ;

void init(){///杨辉三角推组合数
    memset(C,0,sizeof(C));
    C[0][0] = 1 ;
    for( int i = 1 ; i <= 500 ; i++ ){
        C[i][0] = 1 ;
        for( int j = 1 ; j <= i ; j++ ){
            C[i][j] = C[i-1][j] + C[i-1][j-1] ;
            if( C[i][j] > mod ) C[i][j] -= mod ;
        }
    }
}

LL work(){
    ///分解质因数
    for( int i =1 ;i<=n;i++ ){
        for( int j=2;j*j<=b[i];j++ ){
            while( b[i]%j==0){
                m[j]++;
                b[i]/=j ;
            }
        }
        if(b[i]> 1) m[b[i]]++ ;
    }
    map<int,int>::iterator it ;

    LL ans =1;///全部情况数
    for( it=m.begin();it!=m.end();it++ ){
        int num =it->second ;/// M_i 质因数的个数
        ans = (ans*C[num+n-1][n-1])%mod ;
    }
    for( int i =1;i<n;i++ ){///容斥
        LL val =C[n][i] ;
        for( it = m.begin();it!=m.end();it++){
            int num = it->second ;
            val=( val*C[num+n-i-1][n-i-1])%mod ;
        }
        if(i&1)
            ans=(ans-val+mod)%mod ;
        else
            ans=(ans+val)%mod ;
    }
    return ans ;

}
int main(){
    init() ;
    while( cin >> n ){
        m.clear();
        for( int i = 1 ; i <= n ; i++ )
            cin >> b[i] ;
       cout<<work()<<endl;
    }
    return 0;
}
THE END GAME
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
暑假集训日记——7.24（牛客+一堆算法）
weixin_43870114的博客
07-28 95
筱玛爱游戏 题解： 这题需要一些线性代数的知识 每个数可以看做一个 向量(即每一维都是 0 或 1 的向量) 这时数的异或就相当于向量的加法 那么集合存在一个非空子集异或和为0即为这个向量组线性相关 那么两个人在博弈过程中每一步都需保证向量组线性无关 那么这个向量组最大的大小即为所有向量的秩 而由线性代数基本结论，若当前选出的向量线性空间维数小于所有向量的秩，一定能加入一个另外的向量，使得向量组仍...
HDU - 1043 Eight （八数码）逆向BFS+康拓展开hash判重+打表
laaahu的博客
07-19 800
HDU - 1043 Eight （八数码）逆向BFS+康拓展开hash判重+打表 康拓展开 不理解康拓展开的可以学习一下啊：https://blog.csdn.net/laaahu/article/details/96425860 题意： 在九个方格中每个方格都有一个数字，其中一个是x（后面我们会用0表示这个点），所以就是0到8这九个数字在一个3X3的方格里运动，每次x与其相邻的上下左右处的数字...
康拓展开+康拓展开的逆运算
applehth的博客
05-16 367
康拓展开： 是一个全排列到一个自然数的双射，常用于构建哈希表时的空间压缩。康托展开的实质是计算当前排列在所有由小到大全排列中的顺序(名次)，因此是可逆的。 公式 X=a[n]*(n-1)!+a[n-1]*(n-2)!+…+a[i]*(i-1)!+…+a[1]*0! 其中,a[i]为整数,并且0 举例 例如，一个数组s={1,2,3,4,5,6,7,8,9}，它的一个排列s1={3 ,5
POJ - 1077 Eight bfs+康拓展开
a302549450的博客
06-26 157
题目链接题意：就是一个八数码问题，求从初始状态得到解的路径。思路：bfs，对于每个状态通过康拓展开来记录判重。康拓展开：设一个状态a[1...n],其康拓展开X=b[n]*(n-1)!+b[n-1]*(n-2)!+...+b[i]*(i-1)!+....+b[1]*0!;其中b[i]表示在a[i+1]到a[n]中&lt;a[i]的个数。逆康拓展开：因为b[i]&lt;i，所以(n-1)!&gt;b...
ACdream 1076 XXX的机器人 （dp+置换+康拓展开）
LzyRapX
07-20 475
题目链接： ACdream 1076题意： 中文题面。。。题解： dp dp 题。 设dp[i][j] dp[i][j]表示指令为i  i 时，全排列状态为j  j 。因为每条指令可以选择执行或不执行。那么 dp dp 的转移方程就是： dp[i+1][x]=min(dp[i+1][x],dp[i][j]+|t−s|+1) dp[i+1][x] = min(dp[i+1][x],dp[i
[LeetCode] 使字符串有序的最少操作次数 思维+康拓展开+组合数学
Garett的博客
04-19 354
题目链接：使字符串有序的最小操作次数 题面： 思路： 首先要想到题面所述的操作其实是将s转化为字典序恰好比当前排列小的一个排列。 那么问题就转化为当前排列在s的所有字符组成的全排列中的次序问题。 自然我们会想到用康拓展开，但由于s中有重复出现的字符，所以我们还需要用到多重集的全排列公式。此外由于存在取模操作，又存在除法所以还需要用到逆元。 代码： class Solution { public static final long MOD=1000000007; publ
HDU 1043 八数码 反向BFS+康拓展开
zhenlingcn的博客
07-16 283
题意八数码，末状态”12345678x”题解POJ八数码用普通BFS+康拓展开过了。但是HDU1043却过不了，看了题解才知道要反向BFS打表。有关康拓展开，可以参考相关百度文库文档。反向BFS的思路就是从末状态开始BFS搜索到达所有状态的路径，八数码问题总共也就三十多万个状态，每个状态路径可以用string进行存储，并使用康拓展开进行编码。针对每一组数据，进行编码后直接到string数组中查找即可
HDU Eight （BFS+康拓展开）
weixin_58110654的博客
09-23 128
HDU1043 Eight 题目 给出数字八数码游戏的一种棋局状态，求完成游戏的步骤，输出 任意一种即可 思路 每种对弈状态可以抽象成一个19的数字排列，所以所有的可能就是19的全排列，总共是10610^6106级，因此可以穷举。 考虑搜索，从初始态穷举所有状态。首先，如何搜索，如果使用BFS，显然会造成重复，并且无法通过简单的打标记方法解决。 但是，我们可以通过对他们的全排列进行编号，这样就可以采用标记解决。 那么，解决如何快速获取一对弈状态的全排列编号的问题。 我们可以使用康拓展开,排列的编号是倒
【 HDU1043-经典BFS+康拓展开 八数码】 （待更）
sxy201658506207的博客
12-01 237
给定一个序列，由1~8数字和字母x组成，表示的是一个3*3的矩形。每次操作x都能与相邻的数字交换，问如何操作才能使得序列为{1,2,3,4,5,6,7,8,x}。   //多组数据-需要计算全部路径后直接输出 //反向搜索+打表(离线) #include&lt;iostream&gt; #include&lt;cstring&gt; #include&lt;cstdio&gt; #inclu...
康拓展开的八数码广搜程序
11-28
康拓展开的八数码广搜程序 解决了搜索不能判断重复的难题 基于poj1077题目制作
康拓展开和逆康拓展开
09-22
ACM算法里的康拓展开和逆康拓展开，适合ACM竞赛的小白
康拓1数据延期发卡教程 康拓变种工具
06-13
康拓1数据延期发卡教程 康拓变种工具
康拓管理发卡器-园区码A1B2C3.exe
08-18
康拓梯控系统滚动码发卡延期详细教程，康拓梯控系统延期，突破电梯卡楼层限制，实现电梯卡通所有楼层，滚动码效验位计算教程，需要生成其他园区码的联系我
康拓系统分析
05-02
康拓系统分析
等比数列三角形
weixin_43870114的博客
05-15 1555
题解： （3|12）3能被12整除 #include<bits/stdc++.h> using namespace std; #define LL long long #define endl '\n' const int mod=1e9+7; LL n; double e=(1.0+sqrt(5.0))/2; LL b[5000]; int main() { whi...
普通递归关系
weixin_43870114的博客
05-29 545
注意：0 0的特判 纠结好久…哎 #include<algorithm> #include<cstdio> #include<cmath> #include<iostream> #include<cstring> #include<set> #include<map> #include<iterator&g...
求解前n个数的约数和
weixin_43870114的博客
04-03 463
#include<cstdio> int n,ans; int main() { scanf("%d",&n); for(int i=1,j;i<=n;i=j+1)//节约时间 { j=n/(n/i); ans+=(n/i)*(j-i+1); } printf("%d",ans); retur...
暑假集训日记——7.7（月赛+scc习题）
weixin_43870114的博客
07-08 413
A - 迷宫城堡 题解： 两两可达，就是判断一下这个图是否是一个强联通图，模板题。 #include<bits/stdc++.h> using namespace std; typedef long long ll; typedef pair<int, int> pii; typedef pair<ll, int> pli; typedef pair<...
康拓123发卡系统rar
最新发布
11-05
康拓123发卡系统rar是一个发卡系统的压缩文件，使用RAR格式进行压缩。康拓123发卡系统是一款可以用于发行各类卡片的系统软件，它提供了卡片的生成、管理、销售等功能。利用这个系统可以方便地制作各种类型的卡片，如...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值