hdu--4407--一不留神就TLE了

相比于上次我做过的一个容斥题 和他很相似 就是多了一个modify操作

不过也不难 只要 通过 加加减减 操作就能完成了

这个操作 我想过很多种方法 最后觉得用Map迭代访问是最好的选择

如果遍历x->y这个区间 因为y<=400000可能会太大

如果我记录操作1之前执行了多少次的操作2  然后来访问这些操作2  但是可能同一个位置上的数据被多次修改 这样统计也会出错

这边 使用map还有一个原因 就是操作总共也就1000次 那么我最差的情况也就是修改了999次 询问了1次 那map里面最多也就是插入了999个元素 很快的

------

然后我的tle也是发生在了map身上  还是我自己的思维不好 这个加加减减操作被我搞复杂了 后来去看了别人的 瞬间 fml

void modify( LL x , LL y , LL p )
{
    LL a , b c , d;
    map<LL,LL>::iterator it;
    for( it=mp.begin();it!=mp.end();it++ )
    {
        c = it->first;
        d = it-second;
        a = gcd( c , p );
        b = gcd( d , p );
        if( it->first == it->second || it->first<=x-1 || it->first >y )
        {
            continue;
        }
        else if( a==1 )
        {
            if( b==1 )
                ans = ans - c + d;
            else
                ans = ans - c;
        }
        else
        {
            if( b==1 )
                ans = ans + d;
        }
    }
}

void modify( LL x , LL y , LL p )
{
    map<LL,LL>::iterator it;
    for( it=mp.begin();it!=mp.end();it++ )
    {
        if( it->first == it->second || it->first<=x-1 || it->first >y )
        {
            continue;
        }
        if( gcd(it->first,p)==1 )
            ans-=it->first;
        if( gcd(it->second,p)==1 )
            ans+=it->second;
    }
}

第一个 罗里吧嗦的版本就是我写的 真TM 又臭又长啊  害我TLE了

#include <iostream>
#include <map>
#include <vector>
using namespace std;

typedef long long LL;
LL ans;
vector<int> ve;
map<LL,LL>mp;
void getPrime( int n )
{
    for( int i = 2 ; i<=n/i ; i++ )
    {
        if( n%i==0 )
        {
            ve.push_back( i );
            while( n%i==0 )
                n /= i;
        }
    }
    if( n>1 )
        ve.push_back( n );
}

LL gcd( LL x , LL y )
{
    return x % y == 0 ? y : gcd( y , x%y );
}

LL solve( LL x )
{
    LL ans = 0 , var , num;
    int veSize = ve.size() , cnt;
    for( int i = 1 ; i<( 1<<veSize ) ; i++ )
    {
        var = 1;
        cnt = 0;
        for( int j = 0 ; j<veSize ; j++ )
        {
            if( i&(1<<j) )
            {
                var = (LL) (var * ve[j]);
                ++ cnt;
            }
        }
        num = x / var;
        if( cnt&1 )
            ans += num * ( var + num * var ) / 2;
        else
            ans -= num * ( var + num * var ) / 2;
    }
    return ans;
}

void modify( LL x , LL y , LL p )
{
    map<LL,LL>::iterator it;
    for( it=mp.begin();it!=mp.end();it++ )
    {
        if( it->first == it->second || it->first<=x-1 || it->first >y )
        {
            continue;
        }
        if( gcd(it->first,p)==1 )
            ans-=it->first;
        if( gcd(it->second,p)==1 )
            ans+=it->second;
    }
}

int main()
{
    cin.sync_with_stdio(false);
    int t , n , m , op;
    LL x , y , p , c , ans1 , ans2 , sum1 , sum2;
    cin >> t;
    while( t-- )
    {
        cin >> n >> m;
        mp.clear();
        while( m-- )
        {
            cin >> op;
            if( op==1 )
            {
                cin >> x >> y >> p;
                ve.clear( );
                getPrime( p );
                sum1 = (x-1) * x / 2;
                ans1 = solve( x-1 );
                sum2 = y * (y+1) / 2;
                ans2 = solve( y );
                ans = sum2 - ans2 - ( sum1 - ans1 );
                modify( x , y , p );
                cout << ans << endl;
            }
            else
            {
                cin >> x >> c;
                mp[x] = c;
            }
        }
    }
    return 0;
}

这种题目 还有一个很烦的地方就是 要用到long long或者__int64来处理

today:

　　对于 性格 外向的女生 毫无 招架之力