T1
思路:
容斥原理+快速幂
设t[i]表示以k*i为最小公约数的选数方案
可以设上下界l1=a/(k*i),r1=b/(k*i)
其实不算最大公约数大于k*i的情况下,就等于(r1-l1+1)^n-(r1-l1+1)
代码:
#include<cstdio>
#include<iostream>
using namespace std;
const int m=1000000007;
const int maxn=100000;
int n,k,a,b,l,r,t[200020];
int ksm(int x,int k)
{
int r=1;
while (k)
{
if (k&1) r=(1LL*r*x)% m;
k>>=1;
x=(1LL*x*x)% m;
}
return r;
}
int main()
{
freopen("number.in","r",stdin);
freopen("number.out","w",stdout);
scanf("%d%d%d%d",&n,&k,&a,&b);
l=a/k; r=b/k;
if (a%k) l++;
for (int i=maxn;i>=1;i--)
{
int l1=l/i,r1=r/i;
if (l%i) l1++;
if (l1<=r1)
{
t[i]=ksm(r1-l1+1,n);
t[i]=(t[i]-(r1-l1+1)+m)%m;
for (int j=i*2;j<=maxn;j+=i) t[i]=(t[i]-t[j]+m)% m;
}
}
if (l==1) t[1]=(t[1]+1)%m;
printf("%d",t[1]);
return 0;
}
T2
Description
Alice总是会提出很多奇怪的问题,一天他让他的朋友Bob跟他一起研究一个奇怪的问题。问题是:[A,B]中有多少个数满足组成这个数的数字之和为S,另一个问题是[A,B]内满足这一要求最小的数是哪个?
编程帮Bob解决这个问题。
Input
输入包含三个整数A,B和S(1<=A<=B<10^15,1<=S<=135)。
Output
第一行输出一个整数表示满足条件的数的个数;
第二行输出最小的满足条件的数。
输入文件保证至少存在一个数。
思路:
dp
设f[i][j][0/1][0/1]表示前i位,和为j,第一个0/1表示是否等于(0)或大于a(1),第二个0/1表示是否等于(0)或小于b(1)
那么求最小值,我们还可以定义一个与上面差不多的方程t[i][j][0/1][0/1]
这样可以用同样的时间做两件事
代码:
uses math;
var a,b:int64;
i,j,k,l,s,w,len,l1,l2,k1,r1:longint;
x,y:array[0..16]of longint;
f,t:array[0..16,0..136,0..1,0..1]of int64;
procedure main;
begin
for i:=1 to len do for j:=0 to s do for k:=0 to 1 do for l:=0 to 1 do
if f[i,j,k,l]>0 then
begin
if k=0 then l1:=x[i] else l1:=0;
if l=0 then r1:=y[i] else r1:=9;
for w:=l1 to r1 do
if (j+w<=s) then
begin
if w>x[i] then k1:=1 else k1:=k;
if w<y[i] then l2:=1 else l2:=l;
inc(f[i+1,j+w,k1,l2],f[i,j,k,l]);
if (t[i,j,k,l]*10+w)<t[i+1,j+w,k1,l2] then t[i+1,j+w,k1,l2]:=t[i,j,k,l]*10+w;
end;
end;
writeln(f[len+1,s,0,0]+f[len+1,s,0,1]+f[len+1,s,1,0]+f[len+1,s,1,1]);
writeln(min(min(min(t[len+1,s,0,0],t[len+1,s,0,1]),t[len+1,s,1,0]),t[len+1,s,1,1]));
end;
procedure init;
var i:longint;
s1:string;
begin
read(a,b,s);
str(b,s1);
len:=length(s1);
for i:=1 to len do y[i]:=ord(s1[i])-48;
str(a,s1);
while (length(s1)<len) do s1:='0'+s1;
for i:=1 to len do x[i]:=ord(s1[i])-48;
fillchar(t,sizeof(t),127);
f[1,0,0,0]:=1;
t[1,0,0,0]:=0;
end;
begin
init;
main;
end.
T3
Description
Alice又想到一个游戏:N个数每个数都在0到9之间,可以对每一个数进行加1操作,但这个加1比较特别,0-8加1后会相应变成1-9,但9加1后会变成0,给出N个数,进行M次操作,每次操作都会给出两个整数A和B(1<=A<=B<=N),要求输出第A个数到第B个数的和,然后把A到B之间的每一个数进行一次加1操作。
你能帮助Alice吗?
Input
输入文件第一行包含两个整数N和M(1<=N<=250000,1<=M<=100000);
第二行N个数字,每个数字都在0到9之间,中间没有空格,表示初始值。
接下来M行,每行包含两个整数A和B(1<=A<=B<=N)
Output
输出M行,每行表示对应区间的和,注意是先计算后操作。
思路
线段树
查询,区间和。
代码
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
#define N 250010
using namespace std;
char c[N];
int n,m;
int t[N<<2],tree[N<<2][10];
void rev(int* x,int k)
{
if(k>10) k%=10;
static int v[10];
memcpy(v,x+(10-k),k<<2);
memmove(x+k,x,(10-k)<<2);
memcpy(x,v,k<<2);
}
void pd(int m,int x)
{
if(t[x])
{
t[x<<1]+=t[x];
t[x<<1|1]+=t[x];
rev(tree[x<<1],t[x]);
rev(tree[x<<1|1],t[x]);
t[x]=0;
}
}
int sum(int x)
{
int ans=0;
for(int i=1;i<10;++i) ans+=i*tree[x][i];
return ans;
}
int query(int l,int r,int x,int head,int tail)
{
if(head<=l&&r<=tail) return sum(x);
int mid=l+r>>1,ans=0;
pd(r-l+1,x);
if(head<=mid) ans+=query(l,mid,x<<1,head,tail);
if(mid<tail) ans+=query(mid+1,r,x<<1|1,head,tail);
return ans;
}
void ps(int x)
{
for(int i=0;i<10;++i) tree[x][i]=tree[x<<1][i]+tree[x<<1|1][i];
}
void build(int l,int r,int x)
{
if(l==r){ tree[x][c[l]-'0']++; return; }
int mid=l+r>>1;
build(l,mid,x<<1);
build(mid+1,r,x<<1|1);
ps(x);
}
void plu(int l,int r,int x,int head,int tail)
{
if(head<=l&&r<=tail)
{
t[x]++; rev(tree[x],1); return;
}
int mid=l+r>>1; pd(r-l+1,x);
if(head<=mid) plu(l,mid,x<<1,head,tail);
if(mid<tail) plu(mid+1,r,x<<1|1,head,tail);
ps(x);
}
int main()
{
scanf("%d%d%s",&n,&m,c+1);
build(1,n,1);
for(int l,r,i=0;i<m;++i)
{
scanf("%d%d",&l,&r);
printf("%d\n",query(1,n,1,l,r));
plu(1,n,1,l,r);
}
}