JZOJ2017.08.14 B组

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); 
    }
}