NOIP2011提高组解题报告day2.docVIP

计算系数 【问题描述】给定一个多项式(ax + by)^k，请求出多项式展开后(x^n)*(y^m)项的系数。【输入】输入文件名为 factor.in。共一行，包含 5 个整数，分别为a，b，k，n，m，每两个整数之间用一个空格隔开。【输出】输出文件名为 factor.out。输出共 1 行，包含一个整数，表示所求的系数，这个系数可能很大，输出对10007 取模后的结果。【输入输出样例】 factor.in factor.out 1 1 3 1 2 3 【数据范围】 对于 30%的数据，有0≤k≤10；对于 50%的数据，有a = 1，b = 1；对于 100%的数据，有0≤k≤1,000，0≤n, m≤k，且n + m = k，0≤a，b≤1,000,000。 【考察知识点】数论/二项式系数/组合【思路】学OI的应该知道二项式系数/组合数/杨辉三角，再不行手推几个也就能知道ans=C(k,n)*(a^n)*(b^m)这里由于是取模，同余只能用于加和乘。但又要算组合数，而组合数公式里有除法，我因此纠结了很久很久，甚至最后都有写高精加减乘除模的想法。瞬间灵光一闪，C(k,n)=C(k-1,n)+C(k-1,n-1)，行了！注意要边加边mod，否则跟没想到一样。【时间复杂度】O(k^2) program factor; //uses sysutils; const ????proname=factor; ????yu=10007; var ????fin,fout:text; ????i,j,k,m,n,l,r,s,t,x,y:longint; ????a,b:int64; ????ans:int64; ????c:array[0..1100,0..1100]of longint; ?? procedure pin; var ????i,j:longint; begin ????readln(fin,a,b,k,n,m); end; ?? function getc(n,r:int64):int64; var ????i,j:longint; begin ????fillchar(c,sizeof(c),0); ????c[0,0]:=1; ????c[1,0]:=1; ????c[1,1]:=1; ????for i:=2 to n do ????????begin ????????c[i,0]:=1; ????????for j:=1 to i do ????????????c[i,j]:=(c[i-1,j]+c[i-1,j-1]) mod yu; ????????end; ????exit(c[n,r] mod yu); end; ?? procedure main; var ????i,j:longint; begin ????ans:=1; ????a:=a mod yu; ????b:=b mod yu; ????for i:=1 to n do ????????ans:=ans*a mod yu; ????for i:=1 to m do ????????ans:=ans*b mod yu; ????if nm then ????????t:=m ????else ????????t:=n; ????ans:=ans*getc(k,t) mod yu; end; ?? procedure pout; var ????i,j:longint; begin ????writeln(fout,ans); end; ?? begin ????assign(fin,proname+.in); ????assign(fout,proname+.out); ????reset(fin); ????rewrite(fout); ????//time:=now; ????pin; ????main; ????pout; ????//writeln(fout,(now-time)*24*3600*1000:0:0); ????close(fin); ????close(fout); end. 聪明的质监员 【问题描述】小 T 是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有n 个矿石，从1到n 逐一编号，每个矿石都有自己的重量wi 以及价值vi。检验矿产的流程是：1、给定m 个区间[Li，Ri]；2、选出一个参数W；3、对于一个区间[Li，Ri]，计算矿石在这个区间上的检验值Yi ：Yi=sum*sumv。sum为[Li，Ri]之间Wi=W的矿石数量，sumv为其价值和。这批矿产的检