【扬中集训DAY1T1】 微信群

【题目链接】

          点击打开链接

【算法】

         对问题稍加分析后，发现其实要求的就是 ：

         C(N,K) + C(N,K+1) + C(N,K+2) + ... + C(N,N)

         因为N最大10^9,K最大10^5，很显然，我们不可能通过暴力或杨辉三角等来解决此题

         我们发现 ：

         C(N,K) + C(N,K+1) + C(N,K+2) + ... + C(N,N)

         = C(N,1) + C(N,2) + C(N,3) + ... + C(N,N) - C(N,1) - C(N,2) - ... - C(N,K-1) 

         =  2^N - 1 - C(N,1) - C(N,2) - ... - C(N,K-1)

         同时，我们可以推出 ： C(N,K) = C(N,K-1) * (N - K + 1) / K

         我们可以通过这个递推式在O(K)的时间求出C(N,1)..C(N,K-1) 

         但是问题还没有完，由于这个式子要模10^9+7

         因此在计算 C(N,K-1) * (N - K + 1) / K时，除法运算会变得十分棘手，我们怎么解决这个问题呢？

         其实很简单，(a / b) mod P = （a * b的逆元）mod P,有了这个小结论，问题就迎刃而解了！

         关于这个结论，具体的证明我就不给出了，详见这篇文章 :

         http://blog.csdn.net/star_weeper/article/details/52973879

       【代码】

        

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
  
const ll MOD = 1e9 + 7; 

ll i,N,K,a=1,b=1,ans;

template <typename T> inline void read(T &x) {
		ll f = 1; x = 0;
		char c = getchar();
		for (; !isdigit(c); c = getchar()) { if (c == '-') f = -f; }
		for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
		x *= f;
}

template <typename T> inline void write(T x) {
		if (x < 0) { putchar('-'); x = -x; }
		if (x > 9) write(x/10);
		putchar(x%10+'0');
}

template <typename T> inline void writeln(T x) {
		write(x);
		puts("");
}

inline ll pow_mod(ll a,ll n) {
		ll ans;
		if (!n) return 1;
		if (n == 1) return a;
		ans = pow_mod(a,n>>1);	
		ans = ans * ans % MOD;
		if (n & 1) ans = ans * a % MOD;
		return ans; 
}

inline void exgcd(ll a,ll b,ll &x,ll &y) {
		if (b == 0) {
				x = 1;
				y = 0;
		} else {
				exgcd(b,a%b,y,x);
				y -= a / b * x;
		}	
}

inline ll inv(ll a) {
		ll x,y;
		exgcd(a,MOD,x,y);	
		return (x % MOD + MOD) % MOD;
}

int main() {

		read(N); read(K);
	
		ans = pow_mod(2,N) - 1;
		 
		for (i = 1; i < K; i++) {
	   		    a = (a * i) % MOD;
				b = ((b % MOD) * ((N - i + 1) % MOD)) % MOD;
				ans = (ans - (b * inv(a)) % MOD) % MOD;
				if (ans < 0) ans += MOD;	
		}
	
		writeln(ans);
		
		return 0;

}