Flower

1.花
（flower.cpp/c/pas）
【问题描述】
商店里出售n种不同品种的花。为了装饰桌面，你打算买m支花回家。你觉得放两支一样的花很难看，因此每种品种的花最多买1支。求总共有几种不同的买花的方案？答案可能很大，输出答案mod p的值。
【输入格式】
一行3个整数n，m，p，意义如题所述。
【输出格式】
一个整数，表示买花的方案数。
【输入输出样例1】
flower.in   
4 2 5   
flower.out
1
见选手目录下的flower / flower1.in与flower / flower1.out
【输入输出样例1说明】
用数字1,2,3,4来表示花的种类的话，4种花里买各不相同的2支的方案有(1,2)、(1,3)、(1,4)、(2,3)、(2,4)、(3,4)，共6种方案，模5后余数是1。
【输入输出样例2】
见选手目录下的flower / flower2.in与flower / flower2.out
【数据范围】
对于30%的数据，n,m≤10
对于50%的数据，n,m≤1000
对于80%的数据，1≤m≤n≤50,000
对于100%的数据，1≤m≤n≤1,000,000，p≤1,000,000,000

cnm的暴力计算，有除法而且取模，刚开始想拓展欧几里得，但是x与p不一定互质。。。
只能用分解质因数
预处理出质数即可

#include <cstdio>
#include <iostream>
#include <cmath>
#define ll long long
using namespace std;
ll num[1100000],maxn,p;
ll prime[99999],cnt;
bool vis[1100000];
void Prime(int n)
{
    vis[1]=1;
    for(ll i=2;i<=n;i++)
     {
        if(!vis[i])
         prime[++cnt]=i;

        for(ll j=1;j<=cnt&&i*prime[j]<=n;j++)
         {
            vis[i*prime[j]]=1;
            if(!(i%prime[j]))
             break;
         }   
     }
}
ll fast_pow(ll x,ll y,ll z)
{
    ll tmp=1;
    while(y)
    {
        if(y&1)
        tmp=tmp*x%z;
        x=x*x%z;
        y>>=1;
    }
    return tmp;
}
int main()
{
    ll n,m;
    scanf("%lld%lld%lld",&m,&n,&p);
    Prime(m);
    for(ll i=n+1;i<=m;i++)
    {
        ll x=i;
        for(ll j=1;j<=cnt;j++)
         {
            if(!vis[x])
            {
                num[x]++;
                break;
            }
            while((x%prime[j])==0)
             x/=prime[j],num[prime[j]]++; 
            if(x==1) break;
         }
        //if(!vis[i]) num[i]++;
    }
    for(ll i=2;i<=(m-n);i++)
    {
        ll x=i;
        for(ll j=1;j<=cnt;j++)
         {
            if(!vis[x])
            {
                num[x]--;
                break;
            }
            while((x%prime[j])==0)
             x/=prime[j],num[prime[j]]--; 
            if(x==1) break;
         }

        //if(!vis[i]) num[i]--;
    }
    ll ans=1;
    //printf("%lld\n",cnt);
    for(ll i=1;i<=cnt;i++)
     ans=(ans%p*fast_pow(prime[i],num[prime[i]],p)%p)%p;
    printf("%lld",ans);
}