数学模板 - 组合数学

排列组合

生成1~n的排列：

void solve(int n, int *a, int num)///solve(n,mem,0),输入n,按字典序从小到大的顺序输出前n个数的所有排列
{
    if(num == n)///递归边界
    {
        for(int i = 0; i < n - 1; i++)
            printf("%d ",a[i]);
        printf("%d\n",a[n-1]);
    }
    else
    {
        for(int i = 1; i <= n; i++)///在a[num]中遍历一遍所有数
        {
            int ok = 0;
            for(int j = 0; j < num; j++)
            {
                if(a[j] == i)///判断a[0]~a[num-1]是否出现过该数字,出现标记为1
                    ok = 1;
            }
            if(!ok)
            {
                a[num] = i;
                solve(n,a,num+1);
            }
        }
    }
    return;
}

生成可重集的排列：

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

const int Max = 100000 + 10;
int mem[Max], q[Max];

void solve(int n, int *p, int *a, int num)///输入数组q,并按字典序从小到大的顺序输出数组a各元素的所有全排列
{
    if(num == n)///递归边界
    {
        for(int i = 0; i < n - 1; i++)
            printf("%d ",a[i]);
        printf("%d\n",a[n-1]);
    }
    else
    {
        for(int i = 0; i < n; i++)///在a[num]中遍历一遍p数组中所有的数
        {
            int bj1 = 0, bj2 = 0;
            if(!i || p[i] != p[i-1])///避免重复
            {
                for(int j = 0; j < num; j++)    if(a[j] == p[i])    bj1++;;
                for(int j = 0; j < n; j++)    if(p[i] == p[j])    bj2++;
                if(bj1 < bj2)///避免遗漏
                {
                    a[num] = p[i];
                    solve(n,p,a,num+1);
                }
            }
        }
    }
    return;
}

int main()
{
    int n;
    while (~scanf("%d",&n))
    {
        for(int i = 0; i < n; i++)
            scanf("%d",&q[i]);
        sort(q,q+n);
        solve(n,q,mem,0);
    }
    return 0;
}

next_permutation排列函数

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

const int Max = 100000 + 10;
int n, mem[Max];

void solve_int()///int 类型的 next_permutation 使用
{
    for(int i = 0; i < n; i++)
        scanf("%d",&mem[i]);
    sort(mem,mem+n);
    do
    {
        for(int i = 0; i < n - 1; i++)
            printf("%d ",mem[i]);
        printf("%d\n",mem[n-1]);
    }while(next_permutation(mem,mem+n));
}

void solve_char()///char 类型的 next_permutation 使用
{
    int len;
    char ch[Max];
    scanf("%s",ch);
    len = strlen(ch);
    sort(ch,ch+len);
    do
    {
        printf("%s\n",ch);
    }while(next_permutation(ch,ch+len));
}

void solve_string()///string 类型的 next_permutation 使用
{
    string str;
    while(cin >> str && str != "*")
    {
        sort(str.begin(),str.end());
        do
        {
            cout << str << endl;
        }
        while(next_permutation(str.begin(),str.end()));
    }
}

bool cmp(char a, char b)///自定义排序顺序,'A' < 'a' < 'B' < 'b' ... < 'Z' < 'z'
{
    if(tolower(a) == tolower(b))
        return a < b;
    else
        return tolower(a) < tolower(b);
}
void solve_mem()///自定义排序方法后 next_permutation 的使用
{

    int len;
    char ch[Max];
    scanf("%s",ch);
    len = strlen(ch);
    sort(ch,ch+len,cmp);
    do
    {
        printf("%s\n",ch);
    }while(next_permutation(ch,ch+len,cmp));
}

int main ()
{
    scanf("%d", &n);
    solve_int();
    solve_char();
    solve_string();
    solve_mem();
    return 0;
}

容斥原理

#include <bits/stdc++.h>
using namespace std;
int prime[40000], m;
bool f[40000];
vector<int>p;///存放质因数

int init()///用筛法40000以内的质数
{
    m = 0;
    for(int i = 2; i < 40000; i++)
    {
        if (f[i] == 0)
            prime[m++] = i;//质数

        for(int j = 0; j < m && i * prime[j] < 40000; j++)
        {
            f[i*prime[j]] = 1;
            if(i%prime[j] == 0)
                break;
        }
    }
}

void factor(int n)///对n分解质因数
{
    p.clear();
    for(int i = 0; i < m && prime[i] * prime[i] <= n; i++)
    {
        if(n % prime[i] == 0)
        {
            p.push_back(prime[i]);
            n /= prime[i];
            while(n % prime[i] == 0)
                n /= prime[i];
        }
    }
    if(n > 1)
        p.push_back(n);
}

int solve(int r)///用二进制实现容斥原理,求区间[1,r]内与n互素的数的个数
{
    int sum = 0;
    for (int i = 1; i < ( 1 << p.size()); ++i)
    {
        int mult = 1, bits = 0;
        for (int j = 0; j < p.size(); ++j)
            if(i & (1 << j))
            {
                bits++;///计算i中1的个数,也就是质因数的个数
                mult *= p[j];
            }
        int cur = r / mult;
        if(bits & 1)///若1的个数是奇数则进行加法，否则进行减法
            sum += cur;
        else sum -= cur;
    }
    return r - sum;
}

int main()
{
    init();
    int n, r;
    while(cin >> n >> r)
    {
        factor(n);
        cout<<solve(r)<<endl;
    }
    return 0;
}

Polya计数

int gcd(int a, int b)
{
    return b == 0 ? a : gcd(b,a % b);
}

ll pow(ll a, ll b)
{
    ll flag = 1;
    while(b--)
    {
        flag *= a;
    }
    return flag;
}

ll polya(ll n, ll m)///用m种颜色对n个珠子构成的环上色,旋转翻转后相同的只算一种
{
    ll ans = 0;
    for(ll i = 0; i < n; i++)
        ans += pow(m,gcd(n,i));
    if(n & 1)
        ans += n * pow(m,(n-1)/2+1);
    else
        ans += (n/2) * pow(m,n/2) + (n/2) * pow(m,n/2+1);
    return ans/(2*n);
}

Burnside定理

void solve(int z)///计算z个数中每个数的循环节长度
{  
    for(int i = 1; i <= z; i++)  
    {  
        int j = num[i];  
        int l = 1;  
        while(i != j)  
        {  
            l++;  
            j = num[j];  
        }  
        mem[i] = l;  
    }  
}  

const int MAX=1001;
int n,perm[MAX], visit[MAX];
void Polya()///sum求循环节个数,Perm用来存储置换,即一个排列
{
    int pos,sum=0;
    memset(visit,0,sizeof(visit));
    for(int i = 0; i < n; i++)
    {
        if(!visit[i])
        {
            sum++;
            pos = i;
            for(int j = 0; !visit[perm[pos]]; j++)
            {
                pos = perm[pos];
                visit[pos] = 1;
            }
        }
    }
    return sum;
}