Fence Building UVALive - 8471 2017乌鲁木齐区域赛D题

计蒜客题目

         一道组队赛题目，当时猜着是有什么规律，但是我们都不太擅长找规律，gg的一道题，今天补题时第一次了解到欧拉定理，这道题目就是用欧拉定理做的。

欧拉公式：V-E+F=2.其中V是顶点(即所有线段的断点数加上交点数),E是边数（即n段椭圆弧加上这些线段被切成的段数），F是面数（即土地块数加上椭圆外那个无穷大的面）。

这道题目让我们求最多能分割几个平面，因此这里 +1 就可以了（圆以外的平面不算），点数就是C（n,4）,边数是C（n，2），因此答案就是 C（n,4）+ C（n，2）+ 1。这里n的范围比较大，因此在求C（n，4）和C（n，2）的时候就要用到卢卡斯定理了。

#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
#include <cstdio>
#include <string>
#include <cstring>
#include <sstream>
#include <cmath>
#include <stack>
#include <map>
#define ll long long
#define ull unsigned long long
#define INF 0x3f3f3f3f
#define mod 1000000007
using namespace std;
ll _pow(ll a,ll b)
{
    ll sum = 1;
    while(b)
    {
        if(b & 1) sum = sum * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return sum;
}
ll c(ll n,ll m)
{
    if(m > n) return 0;
    ll sum = 1,a,b;
    m = min(m,n - m);
    for(int i=1;i<=m;i++)
    {
        a = (n - m + i) % mod;
        b = i % mod;
        sum = sum * (a * _pow(b,mod-2) % mod) % mod;
    }
    return sum;
}

ll lucas(ll n,ll m)
{
    if(m == 0) return 1;
    return c(n % mod,m % mod) * lucas(n / mod,m / mod) % mod;
}
int main()
{
    int t;
    ll n;
    scanf("%d",&t);
    for(int k=1;k<=t;k++)
    {
        scanf("%lld",&n);
        ll ans = (lucas(n,2) + lucas(n,4) + 1) % mod;
        printf("Case #%d: %lld\n",k,ans);
    }
    return 0;
}