HDU 4630 No Pain No Game(树状数组)

题目链接

看的别人的题解，离线之后，按r排序，枚举1-n，利用pre[j]，存上次j的倍数出现的位置，树状数组里统计的当前位置到最后的最大值，树状数组是求区间最值其实应该很麻烦的，但是此题用法只是求到最后的最大值，插入的时候，往前更新就好了，类似求和。

 1 #include <cstdio>

 2 #include <cstring>

 3 #include <algorithm>

 4 using namespace std;

 5 int num[50100];

 6 int p[50100];

 7 int n;

 8 struct node

 9 {

10     int l,r,id;

11 } que[50001];

12 bool cmp(node a,node b)

13 {

14     return a.r < b.r;

15 }

16 bool cmp1(node a,node b)

17 {

18     return a.id < b.id;

19 }

20 int pre[50101];

21 int ans[50101];

22 int lowbit(int t)

23 {

24     return t&(-t);

25 }

26 void insert(int t,int d)

27 {

28     while(t)

29     {

30         p[t] = max(p[t],d);

31         t -= lowbit(t);

32     }

33 }

34 int getmax(int t)

35 {

36     int maxz = 0;

37     while(t <= n)

38     {

39         maxz = max(maxz,p[t]);

40         t += lowbit(t);

41     }

42     return maxz;

43 }

44 int main()

45 {

46     int i,j,t,m,tnum;

47     scanf("%d",&t);

48     while(t--)

49     {

50         memset(p,0,sizeof(p));

51         memset(pre,0,sizeof(pre));

52         scanf("%d",&n);

53         for(i = 1; i <= n; i ++)

54             scanf("%d",&num[i]);

55         scanf("%d",&m);

56         for(i = 0; i < m; i ++)

57         {

58             scanf("%d%d",&que[i].l,&que[i].r);

59             que[i].id = i;

60         }

61         sort(que,que+m,cmp);

62         tnum = 0;

63         for(i = 1; i <= n; i ++)

64         {

65             if(tnum == m) break;

66             for(j = 1; j*j <= num[i]; j ++)

67             {

68                 if(num[i]%j == 0)

69                 {

70                     if(pre[j] != 0)

71                         insert(pre[j],j);

72                     pre[j] = i;

73                     if (j * j== num[i]) continue ;

74                     if(pre[num[i]/j] != 0)

75                         insert(pre[num[i]/j],num[i]/j);

76                     pre[num[i]/j] = i;

77                 }

78             }

79             while(tnum < m&&que[tnum].r == i)

80             {

81                 ans[que[tnum].id] = getmax(que[tnum].l);

82                 tnum ++;

83             }

84         }

85         sort(que,que+m,cmp1);

86         for(i = 0;i < m;i ++)

87         {

88             if(que[i].l == que[i].r)

89             printf("0\n");

90             else

91             printf("%d\n",ans[i]);

92         }

93     }

94     return 0;

95 }